-
Notifications
You must be signed in to change notification settings - Fork 11
Description
Hi
Thank you for your repo and the guidance on calculations for the WHO standard - really very helpful
I am replicating this in python and am matching your values but am struggling with the overlap at 5 y where i am out consistently by a small amount. Anthroplus begins at 60 months (1826.25 days), but anthro finishes at 60.98 months (1856 days). So here there is an overlap of the datasets of not quite a month.
When I look at height in boys, for example, at 1826 days we have L, M and S values of:
anthro/data-raw/growthstandards/lenanthro.txt
Line 1828 in e715811
| 1 1826 1 109.9593 0.04214 H |
but looking at height in boys at 60 months in Anthroplus we have L, M and S values of:
sex age l m s
1 60 1 109.7265 0.04156
And in fact, if we look at the published WHO 2007 tables we find that the data begins at 61 months, not 60. So the values in your Anthroplus library for 60 weeks I figure must have been interpolated from the values in Anthro.
But, If I were to interpolate for the extra 0.25 days I would get
sex age l m s
1 1826.25 1 109.9279 0.0421325
And these values don't match what is in Anthroplus at 60 weeks either.
I have tried to work through your R code but am sorry am more of a pythonista than an R expert and cannot see how you are handling these values between WHO 2006 and 2007 at the 5 y cut off for that month. Any help would be gratefully received. Thank you :)