R的新手,在使用"调查"编码时遇到了两个问题。包裹:
(1)使用"调查"来自svytable的代码包不符合Stata估计的Stata:tab。我打电话给
svyd.nation <- svydesign(ids = ~1, probs = ~wt_national, strata = ~stratum, data=nats.sub)
svytable(formula = ~wpev, design = svyd.nation, Ntotal = 100)
等效的Stata将是:
svyset [pw=wt_national], strata(stratum)
svy: tab wpev
这些为表格提供了不同的%s。关于我在R中做错了什么的任何想法?
(2)我在R中验证的替代权重与Stata中的分布相同 - 由于某种未知原因,似乎不能在R中工作。
svyd.st <- svydesign(ids = ~1, probs = ~wt_state, strata = ~stratum, data=nats.sub)
svytable(~wpev, design = svyd.st, Ntotal=100)
产生:
wpev
0 1
和 svytable(~wpev,svyd.st)
产生
wpev
0 1
Inf Inf
为什么这些替代权重 - 在重置svyset时在Stata中起作用 - 在R中不起作用?任何见解将不胜感激!
非常感谢!