r - 测试 r 调查包中两种方法之间的差异
问题描述
尝试使用 at test 测试调查中两种均值之间的差异时遇到问题。基本上我有不同年份的重复健康调查,我想测试它们之间的平均值差异是否有意义。使用包调查我能够获得方法,以及此代码的标准错误
svyby(~CADL2,by=~ageco+year,design = wcd,FUN= svymean,na.rm=TRUE)
ageco year CADL2None CADL2Limited se.CADL2None se.CADL2Limited
50-54.2007 50-54 2007 0.8717507 0.12824934 0.02652932 0.02652932
55-59.2007 55-59 2007 0.8919843 0.10801569 0.01662038 0.01662038
60-64.2007 60-64 2007 0.9056865 0.09431355 0.01556955 0.01556955
65-69.2007 65-69 2007 0.8525438 0.14745624 0.02376984 0.02376984
70-74.2007 70-74 2007 0.7534787 0.24652131 0.03419399 0.03419399
75-79.2007 75-79 2007 0.7466576 0.25334237 0.04010796 0.04010796
80-85.2007 80-85 2007 0.5690972 0.43090276 0.06083682 0.06083682
85+.2007 85+ 2007 0.3853919 0.61460811 0.08913058 0.08913058
50-54.2017 50-54 2017 0.7150962 0.28490379 0.13929132 0.13929132
55-59.2017 55-59 2017 0.8720697 0.12793025 0.04088908 0.04088908
60-64.2017 60-64 2017 0.8503688 0.14963123 0.01783197 0.01783197
65-69.2017 65-69 2017 0.7931459 0.20685411 0.01829031 0.01829031
70-74.2017 70-74 2017 0.7764609 0.22353912 0.01895070 0.01895070
75-79.2017 75-79 2017 0.6666032 0.33339681 0.02428625 0.02428625
80-85.2017 80-85 2017 0.5462507 0.45374929 0.03324155 0.03324155
85+.2017 85+ 2017 0.3223467 0.67765331 0.03956227 0.03956227
现在我想测试一下,例如,2007 年 50-54 岁年龄组的平均值(0.8717)是否与 2017 年的平均值(0.7151)显着不同。我尝试了 svyttest,但由于我对这个包很陌生,我无法让它工作。
svyby(~CADL2,by=~ageco+year,design = wcd,FUN= svyttest,na.rm=TRUE)
Error in formula[[3]] : subscript out of bounds
我尝试在公式中交换变量,但每次我得到这个错误。你知道我怎样才能让它工作或以任何其他方式实现两种方式之间的测试吗?
解决方案
你想要svyby( CADL2 ~ year , ~ ageco , wcd , svyttest , keep.var = F )
吗?
数据来自的示例?svyttest
library(survey)
data(api)
dclus2<-svydesign(id=~dnum+snum, fpc=~fpc1+fpc2, data=apiclus2)
svyttest(enroll~comp.imp, dclus2)
# single subset
svyttest(enroll~comp.imp, subset(dclus2,stype=='E'))
# with svyby
svyby(enroll~comp.imp,~stype,dclus2,svyttest,keep.var=F)
推荐阅读
- angular - 使用 DomSanitizer 后拒绝加载图像
- html - 确定不同元素的位置
- java - 尽管没有抛出任何异常,但 ObjectOutputStream 不会写入文本文件
- wordpress - 仅在首页 Wordpress 中加载缓慢
- python - 无法在 pyspark 中应用 pandas_udf
- c# - Swashbuckle xml 注释中带有与号的查询字符串
- r - 在R中读取多个文件
- python - 如何使 [] 返回列表派生的子类而不是通常的列表类
- jupyter-notebook - 设置“保存此笔记本时省略代码单元格输出”是什么意思?
- ionic-framework - 离子日期时间在离子 4 中给整个月 31 天?