Я хочу использовать covmat в функции svyby, чтобы сравнить среднее значение общего количества калорий в наборе данных, который у меня есть.
Цель состоит в том, чтобы сравнить средние различия по всем парам (различные диабетические группы: 1. PreDM Unaware, 2. PreDM осознает, 3. DM с лечением, 4. DM без лечения, и 5. общественность без болезней ).
Я уже использовал svyby для получения средних калорий по различным группам статуса. Однако, когда я имитирую тот же код, который предлагается в справочном документе ("http://r -survey.r-forge.r-project.org / survey / html / svyby.html "), инициализируя "covmat" ».
R сообщил об ошибке как:
"covmat = TRUE не реализовано для этого типа дизайна"
TOTAL_SAMPLE_diet<-svydesign(ids = ~SDMVPSU, weights = ~MEC6YR, strata = ~SDMVSTRA, data=TOTAL_SAMPLE, nest=TRUE)
mean_kcal<-svyby(~M_KCAL, by=~as.numeric(STATUS), TOTAL_SAMPLE_diet, svymean, na.rm=TRUE, vartype = "ci", covmat=TRUE)
Результаты опроса:
as.numeric(STATUS) M_KCAL ci_l ci_u
1 1 2032.693 1924.894 2140.493
2 2 2090.285 2057.790 2122.781
3 3 2008.129 1943.930 2072.328
4 4 2008.131 1918.108 2098.155
5 5 2124.223 2094.359 2154.087
Я ожидал увидеть сравнение среднего потребления калорий среди всех статусных групп со значениями p для проверки статистической значимости.
(например, SAS: LSMEANS или MEANDIFF)