Question

我是R的新手，有关于poLCA的问题。当你运行一个潜在的课程在poLCA中进行分析，它为每个受访者提供了一个价值 “属于”每个潜类的后验概率。这些都是存储的作为元素'后验'中的矩阵。

我想创建一个包含每个受访者唯一的数据框 ID号（存储为用于poLCA的数据帧中的变量）和他们的匹配来自'后验'矩阵的后验概率。一世然后，我希望将此数据帧写入csv文件，以便在另一个文件中使用程序

我知道这是可能的，但我似乎无法做到正确（责怪我的 R）无能。任何帮助都会受到非常热烈的赞赏。

祝福。罗伯特德弗里斯

编辑 - 示例

#Loading the poLCA package
library(poLCA)

#The ID variable
serial <-sample(1:1000,100,replace=F)

#The variables used in the latent class model
V1 <- sample(1:2,100,replace=T)
V2 <- sample(1:2,100,replace=T)
V3 <- sample(1:2,100,replace=T)
V4 <- sample(1:2,100,replace=T)

#the data given to the lca
lcadata <- data.frame(serial,V1,V2,V3,V4)

#the lca formula
f <- cbind(V1,V2,V3,V4) ~1

#A 2 class LCA model in poLCA
lca2 <- poLCA(f,lcadata,nclass=3,maxiter=5000,nrep=10)

Answer 1

像

这样简单

write.csv(cbind(lcadata$serial, lca2$posterior), 'new_data.csv', row.names=FALSE)

应该有效。

匹配poLCA R包中的后验概率

1 个答案: