philmikejones / rakeR

Tools for easy spatial microsimulation (raking) in R
http://philmikejones.github.io/rakeR/
11 stars 1 forks source link

Handle zeros gracefully #64

Open philmikejones opened 6 years ago

philmikejones commented 6 years ago

Handle zeros in constraint tables gracefully

philmikejones commented 6 years ago

0s in cons done; what about 0s in aggregaged inds?

philmikejones commented 6 years ago

Error reported by Derrick Atherton by email October 2017:

Error in ipfp::ipfp(x, t(ind_cat), x0 = rep(1, nrow(ind_cat)), maxit = iterations) : 
  BLAS/LAPACK routine 'DGEMV ' gave error code -6

Able to replicate the error when an entire zone has an empty population (i.e. all zeros)