请教一下大家,有份复杂抽样数据采用survey程序包进行分析:
library(foreign)
data1<-read.spss("F://复杂抽样数据分析/data1.sav")
library(survey)
data2 <- svydesign(id=~jcd,strata=~cx,weights =~Wf,data =data1)
如要想知道某一变量加权后的结果时,采用svytable(~F01,data2) 可得总数,svytable(~F01,Ntotal =100,data2) 可得到频率。
但做交叉表时,svytable(~F01+F02,data2) 可得总数,假设F01是两分类变量(1,2),我试用svytable(~F01+F02,Ntotal =200,data2),即定义总数是200。结果是:
F02
F01 NO YES
1 91.01270 12.83669
2 83.79347 12.35714
只能限定总数为200,但无法定义F01为1相加为100,F01为2相加为100,即为F01的分类频率。
R survey帮助中提及:svytable(formula, design, Ntotal = sum(weights(design, "sampling")), round = FALSE,...)
Ntoal后面应该怎么写才能得出F01变量分类的频率,或有survey中有什么其它函数能得出复杂抽样的频数。