logo标签列表

如何在R中计算调查数据的估计标准误差而不使用调查设计？

rdplyrdata.tablesurvey
4

我想知道在给定的例子中，您如何生成标准误差的估计值而不使用调查设计。

library(survey)

data(api)

dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)

例子 1

(total_enroll_mean=svymean(~enroll,by=~stype,design= dstrat, na.rm=TRUE))

        mean     SE
 enroll 595.28 18.509

example 2

(total_enroll_total=svytotal(~enroll,by=~stype,design= dstrat, na.rm=TRUE))

       total     SE
enroll 3687178 114642

示例3

(stype_enroll_mean=svyby(~enroll,by=~stype,design= dstrat, na.rm=TRUE,svymean))

stype  enroll       se
E      416.78    16.41740

H      1320.70   91.70781

M      832.48    54.52157

例子 4

(stype_enroll_total=svyby(~enroll,by=~stype,design= dstrat, na.rm=TRUE,svytotal)

stype     enroll       se
E        1842584.3     72581.33

H         997128.5     69239.40

M         847464.7     55502.96

基本上，我希望能够在R中手动为所有4个示例生成标准误差，而不是使用survey包、srvyr包或任何其他调查相关的包。手动生成svymean和svytotal的标准误差应与上述所有示例中使用survey包生成的标准误差匹配。

1
嗨，您可以键入“debug（survey :: svymean）”来查看每个计算步骤，逐步发生的情况，当您键入“svymean（...）”时。 - Anthony Damico
1
嗨，谢谢回复。但我想知道调查包使用的数学公式来计算标准误差，这样我就可以使用基本的R函数创建自己的函数来计算标准误差估计值。 - Sid Sath
1个回答
0
2
目前你的回答不够清晰，请编辑并添加更多细节，以帮助其他人理解它如何回答问题。你可以在帮助中心找到有关如何编写好答案的更多信息。 - Community
网页内容由stack overflow 提供, 点击上面的
可以查看英文原文，
原文链接