timriffe / adultcoverage Goto Github PK

When I looked at the slopeint function the formula for the intercept is wrong:

(1) intercept <- with(codi,mean(leftterm[age %in% agesfit]) * (1/slope) - mean(rightterm[age %in% agesfit]))

I think it should be:

(2) intercept <- with(codi,mean(leftterm[age %in% agesfit]) - mean(rightterm[age %in% agesfit])) * slope

• add RMSE of age selection to output (mig branch)
• add r2 of line fit to ggb() output (mig branch)
• implement migration in ggb() (mig branch)
• implement migration in seg() (mig branch)
• ensure valid use of mig in ggbseg() (both stages should use mig and in the same way) (mig branch)
• use mig.summed flag (mig branch)
• warn, if total deaths < 25% of population change then either the deaths.summed flag is being used wrong or the user should consider not using the method. Use either message() or warning() (mig branch)
• add more example data (PJ will send spreadsheets)
• implement seg delta method, but think of a better name for it see RD spreadsheet AM_SEG_South Africa_males_3_1.xlsx, calibrate to spirit not values. To iterate or not to iterate?

Why are there points on the graph for ages 0 and 85? The method does not have information to estimate the X and Y variables for these ages.

The item called "coverage" returned by ggb is the coverage relative to census 2 or 1/b. The TDE workbook (http://demographicestimation.iussp.org/content/generalized-growth-balance-method) gives the completeness relative to the average (geometric) coverage of the two censuses. This is the completeness of the intercensal death rates.

The completeness listed in the graph title of ggbChooseAge does not match either of these values. I haven't looked at the code for that.

Hi Tim,

Thank you so much for developing this tool.

I was wondering if it were possible to make an adjustment for net migration before estimating completeness?

Thanks,
Nick