Я анализирую Ассоциацию по оценке глобального гендерного разрыва и антропометрических нарушений у детей (с Rstudio). Для этой цели я рассчитал z-баллы для антропометрических ошибок, например, задержки в росте, а затем вычислил распространенность задержки роста в этой соответствующей стране. Данные для расчета детской антропометрической недостаточности я получил от DHS, и z-показатели были рассчитаны в R с помощью пакета z-показателей. Сейчас я борюсь с расчетом доверительных интервалов для распространенности задержки роста.
Я использовал пакет «опрос» для сложных планов обследования.
библиотека (обзор)
Сложные параметры расчета образца
DHSdesign<-svydesign(id=CI_sample$v021, strata=CI_sample$v023, weights=CI_sample$wt, data=CI_sample, nest=TRUE)
табличный показатель по регионам
svyby(~stuntprev, ~CountryID, DHSdesign, svymean, vartype=c("se","ci"))
То, что R дал мне тогда, было:
```CountryID stuntprev se ci_l ci_u
NM07 NM07 0.308 0 0.308 0.308
NM13 NM13 0.225 0 0.225 0.225```
Например, распространенность задержки роста в Намибии в 2007 году составляла 30,8%, нижняя и верхняя границы также составляли 30,8%, так как мне изменить команду, чтобы она работала?