svyby 대 svytotal의 여러 요소

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

당신이 다음`svytable`를 사용하지 않아야 교차 분류를 원하는 경우 ?svytotal

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