Pessoal, alguém poderia me ajudar def.lcarandom <- randomLCA(definicao[,1:4], freq=definicao$x, initmodel=def.lca, random=TRUE, blocksize=1, quadpoints=190, probit=TRUE)
summary(def.lcarandom) Classes AIC BIC logLik penlogLik Link 2 3842.735 3896.336 -1911.368 -1911.373 Probit Class probabilities Class 1 Class 2 0.8213 0.1787 Conditional outcome probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.0065 0 0.0011 Class 2 0.9612 1.0000 1 1.0000 Marginal Outcome Probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.1987 0 0.1479 Class 2 0.7253 0.9545 1 1.0000 Loadings 2.76629> def.lcarandom$outcomep [,1] [,2] [,3] [,4] [1,] 5.653499e-50 0.006483299 5.024727e-33 0.001059331 [2,] 9.612131e-01 0.999999659 1.000000e+00 1.000000000
A primeira parte postada acima funciona normalmente, sendo que quando vou utilizar a função outcome.probs ela não funciona, ela fica "rodando" o tempo inteiro e não sai mais disso. def.lcarandom.probs<-outcome.probs(def.lcarandom, boot=T, R=100) Abraços!