Hi all,
I am running a multistate model in RMark that runs fine and gives point estimates. However, it gives an error when I try to get the se.TransitionMatrix. Can you provide any advice on how to resolve this error? At first I got a warning message when I used the 'get.real' function:
>top=ms.results$S.stratum.p.stratum.Psi.stratum.t
> Psilist=get.real(top,"Psi",vcv=TRUE)
#In get.real(top, "Psi", vcv = TRUE) :
#Improper V-C matrix for beta estimates. Some variances non-positive.
And then it gives an error when I run se.TransitionMatix with vcv.real (but it runs fine without requesting vcv.real variance matrices):
>Psi.values=Psilist$estimates
> top.psi=TransitionMatrix(Psi.values[Psi.values$time==2014,],vcv.real=Psilist$vcv.real)
#Error in se.TransitionMat/(TransitionMat * (1 - TransitionMat)) :
#non-conformable arrays
#In addition: Warning message:
#In se.TransitionMat[!is.na(pin)] <- x$se[match(pin[!is.na(pin)], :
# number of items to replace is not a multiple of replacement length
Thanks for your help!
Ellen