使用 MCMCpack 作为潜在变量并从 mcmc 对象中提取潜在分数

Question

我正在使用 R 包

MCMCpack

中的函数

MCMCordfactanal()

和 MCMCirtKd() 生成后验分布并创建潜在因子变量。

我可以使用

summary()

和

dimnames()

检查潜在分数，但我仍然不知道如何提取这些值并将它们变成一个潜在变量。

我运行了以下代码：

fact_posterior <- MCMCordfactanal(~ x1 + x2 + x3 + x4, data=vdem, 
    factors=1, lambda.constraints=list(x1=list(1, "+")), burnin=5000, 
    mcmc=50000, thin=100, verbose=500, L0= 0.5, store.lambda=TRUE, 
    store.scores=TRUE, tune=1.2)

irt_posterior <- MCMCirtKd(t(vdem), dimensions=1, burnin=5000, mcmc=50000, 
    thin=10, B0=.25, store.item=TRUE, item.constraints=list("1"=list(2, "-")))

我期望从

mcmc

和

MCMCordfactanal()

生成的

MCMCirtKd()

对象中得到一个潜在变量