在poLAC模型中使用子集数据帧

时间:2016-03-10 06:12:41

标签: r

我正在使用具有权重,分层和psu的国家数据集。我利用调查数据指南中的代码来创建方差估计。然后,我为我感兴趣的子集创建了一个数据框。现在我想弄清楚如何在我的LCA模型中使用这个数据框。我尝试了一些不同的东西,我绝对卡住了。如何将此数据框应用于我的LCA模型?

顺便说一下,我的LCA在没有子集的情况下运行得非常好,但我必须将子集合并以解释方差。

#Variance Estimation
nh2 <- svydesign(id=~PSU_P, strata=~STRAT_P,
             nest = TRUE,
             weights=~WTFA_SA,
             data= nh)

# subset 
subgrp <- subset(nh2, Black < 2 & SEX > 1)

#LCA codes 
f5<-  cbind( hit_use, hivtest, worrhccost, talkgp, lastvisit=lastvisit +  1) ~ age18_24 + age25_35 + age36_45 + age46_55 +  inc0_34 + inc35_49 + inc50_74 + inc75_99 + inc0_49 + inc50_99

lc3<-poLCA(f5, nh, nclass=3, na.rm = FALSE, nrep=10, maxiter=3000)

0 个答案:

没有答案