1

surveyR のパッケージを使用して複数の設計を組み合わせる方法を理解しようとしています。

設計用語を次の順序で並べてみました。

n<- nrow(my_data)
d<- rep(N/n, n)
f<- rep(n/N, n)

#specifying SRS survey design
srs.design<- svydesign(ids= ~0, strata= NULL,
                       data= my_data, 
                       weight= ~d, fpc= ~f)

#specifying post-stratification survey design
ps.design<- postStratify(design= srs.design,
                         strata= ~postsurvey_strata, 
                         population= N.ps) #where N.ps is the poststrata population distribution


#specifying raking survey design
rake.design<- calibrate(design= ps.design, 
                        formula= ~as.factor(gender_age)+ 
                          as.factor(education)+
                          as.factor(race)+
                          as.factor(income),
                        calfun= "raking", 
                        population= pop.P_sam) #where pop.P_sam is the demographic distributions in the population

#specifying trimming survey design 
trim.design<- trimWeights(design= rake.design, lower= 0.2, upper= 6)

#specifying re-calibration
rerake.design<- calibrate(design= trim.design, 
                          formula= ~as.factor(gender_age)+ 
                            as.factor(education)+
                            as.factor(race)+
                            as.factor(income),
                          calfun= "raking", 
                          population= pop.P_sam)


(私が事後階層化を別のステップとして行う主な理由は、13 の事後階層があるためです。事後階層のいくつかは非常に小さいため、事後階層化変数をレーキングに折りたたむと収束につながるのではないかと懸念しています。問題。)

これが正しいアプローチであるかどうか、さらに簡潔な方法があるかどうかを知りたいです。ありがとう!

4

0 に答える 0