R svytable“inf”结果;与Stata“svy”代码不符

时间:2014-10-26 00:42:02

标签: r statistics

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中不起作用?任何见解将不胜感激!

非常感谢!

0 个答案:

没有答案