[R-br] Ajuda com o pacote randomLCA
Victor Eduardo
victorduca08 em gmail.com
Quarta Janeiro 30 22:39:25 BRST 2013
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!
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130130/f6a9cb41/attachment.html>
Mais detalhes sobre a lista de discussão R-br