0

我尝试使用expss软件包进行调查数据分析,但标准误差、方差和置信区间的结果与survey软件包结果不同。

survey

dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)
svyby(~api99, ~stype, dclus1, svymean)

  stype    api99       se
E     E 607.7917 22.81660
H     H 595.7143 41.76400
M     M 608.6000 32.56064

expss

apiclus1 %>% tab_cells(api99) %>% 
         tab_rows(stype) %>% tab_weight(pw) %>%
         tab_stat_fun(w_mean,w_se) %>% tab_pivot()

 |       |    |       |        | #Total |
 | ----- | -- | ----- | ------ | ------ |
 | stype |  E | api99 |   mean |  607.8 |
 |       |    |       |     se |    1.6 |
 |       |  H | api99 |   mean |  595.7 |
 |       |    |       |     se |    4.7 |
 |       |  M | api99 |   mean |  608.6 |
 |       |    |       |     se |    3.7 |

我如何使用surveypackage inside的功能expss

4

0 回答 0