R调查库均值差异检验

R调查库均值差异检验,r,statistics,survey,R,Statistics,Survey,我目前正在使用R的调查库分析调查数据。我有两个不同时期的样品。我的目标是测试两个加权样本均值之间的差值是否等于0问题:我如何使用R的调查库实现这一点? 我尝试了两种方法: 方法1:创建两个不同的后分层对象。玩具示例: 这种方法的问题是,我不知道如何对“第一”和“第二”对象(这两个svymean对象)进行均值差异测试 方法2:仅创建一个后分层对象。玩具示例: 这种方法的问题是,当我运行svyby只是为了检查创建的平均值时,svyby的输出不是我所期望的。当时间=0的理论加权平均值为0.55时,

我目前正在使用R的
调查
库分析调查数据。我有两个不同时期的样品。我的目标是测试两个加权样本均值之间的差值是否等于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)