Question

我试图使用svytotal来计算每对的人口分布（PUMA，BDSP）。当我使用svyby时，它可以像我期望的那样工作：

> svyby(~NP, ~factor(PUMA)+factor(BDSP), ss12hsf, svymean)
        factor(PUMA) factor(BDSP)        NP         se
7501.0          7501            0 0.8328120 0.10738595
7502.0          7502            0 1.1790473 0.10128726
7503.0          7503            0 1.0707710 0.05241559
7504.0          7504            0 0.9982076 0.16726445
...

然而，当我尝试对svytotal做同样的事情时，它似乎只是做一个列表连接而不是交叉产品：

> svytotal(~factor(PUMA)+factor(BDSP), ss12hsf)
                  total      SE
factor(PUMA)7501  66732  514704
factor(PUMA)7502  66350  570194
factor(PUMA)7503  71129 1597630
factor(PUMA)7504  54108  595059
factor(PUMA)7505  44327  498771
factor(PUMA)7506  40461  424705
factor(PUMA)7507  33792  216696
factor(BDSP)0     50036  867466
factor(BDSP)1    111847 1049884
factor(BDSP)2    111976 1586937
factor(BDSP)3     69279 1112640
factor(BDSP)4     24386  496505
factor(BDSP)5      5841  384093
factor(BDSP)10     3534   88227

有任何想法或建议吗？

Answer 1

有一个如何在?svytotal

中执行此操作的示例

svytotal( ~interaction( var1 , var2 ) , design )

svyby与svytotal的多重因素

1 个答案: