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!