R调查库均值差异检验
我目前正在使用R的R调查库均值差异检验,r,statistics,survey,R,Statistics,Survey,我目前正在使用R的调查库分析调查数据。我有两个不同时期的样品。我的目标是测试两个加权样本均值之间的差值是否等于0问题:我如何使用R的调查库实现这一点? 我尝试了两种方法: 方法1:创建两个不同的后分层对象。玩具示例: 这种方法的问题是,我不知道如何对“第一”和“第二”对象(这两个svymean对象)进行均值差异测试 方法2:仅创建一个后分层对象。玩具示例: 这种方法的问题是,当我运行svyby只是为了检查创建的平均值时,svyby的输出不是我所期望的。当时间=0的理论加权平均值为0.55时,
调查
库分析调查数据。我有两个不同时期的样品。我的目标是测试两个加权样本均值之间的差值是否等于0问题:我如何使用R的调查库实现这一点?
我尝试了两种方法:
方法1:创建两个不同的后分层对象。玩具示例:
这种方法的问题是,我不知道如何对“第一”和“第二”对象(这两个
svymean
对象)进行均值差异测试
方法2:仅创建一个后分层对象。玩具示例:
这种方法的问题是,当我运行
svyby
只是为了检查创建的平均值时,svyby
的输出不是我所期望的。当时间=0的理论加权平均值为0.55时,它得出平均值=0.5714
对于理论平均值与svyby的不同之处,我们将非常感谢。
非常感谢您抽出时间 你在找这个吗?谢谢
library(survey)
q1 = c(1,1,1,1,0, 1,1,0,0,1)
# edited #
group = c(0,0,0,1,1, 2,2,3,3,3)
time = c(0,0,0,0,0, 1,1,1,1,1) #Variable that distinguishes between the samples
df = data.frame(q1, group, time)
svy_design = svydesign(ids = ~1 , data = df)
# edited #
pop_data = data.frame(group = c(0,1,2,3), Freq = c(10,90,20,80))
ps_design = postStratify(svy_design, strata = ~group,pop_data)
svyby(~q1, ~time, ps_design, svymean)
svyttest(q1~time, ps_design)
你在找这个吗?谢谢
library(survey)
q1 = c(1,1,1,1,0, 1,1,0,0,1)
# edited #
group = c(0,0,0,1,1, 2,2,3,3,3)
time = c(0,0,0,0,0, 1,1,1,1,1) #Variable that distinguishes between the samples
df = data.frame(q1, group, time)
svy_design = svydesign(ids = ~1 , data = df)
# edited #
pop_data = data.frame(group = c(0,1,2,3), Freq = c(10,90,20,80))
ps_design = postStratify(svy_design, strata = ~group,pop_data)
svyby(~q1, ~time, ps_design, svymean)
svyttest(q1~time, ps_design)
library(survey)
q1 = c(1,1,1,1,0, 1,1,0,0,1)
# edited #
group = c(0,0,0,1,1, 2,2,3,3,3)
time = c(0,0,0,0,0, 1,1,1,1,1) #Variable that distinguishes between the samples
df = data.frame(q1, group, time)
svy_design = svydesign(ids = ~1 , data = df)
# edited #
pop_data = data.frame(group = c(0,1,2,3), Freq = c(10,90,20,80))
ps_design = postStratify(svy_design, strata = ~group,pop_data)
svyby(~q1, ~time, ps_design, svymean)
svyttest(q1~time, ps_design)