Как получить правильную частотную таблицу, взвешенную и невзвешенную в комплексном обследовании?

Question

Я использую survey package для анализа сложного опроса

Проблема, с которой я сталкиваюсь, когда использую этот код

tab_w<- svytable( ~edu+ smok , design= dsub)

В результате была показана таблица частотвся популяция, а не только выборочная совокупность в моем опросе.

         smok
edu         1         2
       1   26466299  87948351
       2   20898455 102012829

Это таблица, которую я собираюсь сделать

----------------------------------------------
    |     Smok             |  Smok
    | Unweighted           | Weighted
    | Sample size          | population Size
    |   1     |  2         |  1     |  2         
----------------------------------------------
edu |-----------------------------------------                    
  1 |26466299| 87948351    |  200 | 204 
  2 |20898455|102012829    |  400 | 300
-----------------------------------------------

Это включает всю популяцию, а не только мою выборочную совокупность. Можно ли включить в таблицу частот только выборочную совокупность?

Answer 1

Вам не нужен пакет опроса для невзвешенной выборочной таблицы, вы можете просто использовать функцию table или xtabs.

Или, если у вас есть только данные, удобно доступные вобъект обследования

with(model.frame(dsub), table(edu,smok))

В качестве воспроизводимого примера

library(survey)
data(api)


with(apistrat, table(comp.imp, sch.wide))
#stratified sample
dstrat<-svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)
svytable(~comp.imp+sch.wide, design=dstrat)
with(model.frame(dstrat), table(comp.imp, sch.wide))