我想在svyby函数中使用covmat来比较我拥有的数据集的总卡路里的平均值。
目的是比较所有配对（不同糖尿病组）的平均差异：1.未意识到PreDM，2。PreDM意识，3。DM接受治疗，4。DM不接受治疗和5.无疾病的公众）。
我已经使用svyby来获取不同状态组的平均卡路里。但是，当我通过初始化“ covmat”来模仿帮助文档（“ http://r-survey.r-forge.r-project.org/survey/html/svyby.html”）提供的相同代码时。
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）