Adding mxCI() to a model makes all SE's reported as NaN. I don't think it should.
data(demoOneFactor) Rff <- mxMatrix(type="Stand",nrow=1,ncol=1,free=F,name="Rff") L <- mxMatrix(type="Full",nrow=5,ncol=1,free=T,values=0.2,labels=paste("l",1:5,sep=""),name="L") I <- mxMatrix(type="Iden",nrow=5,ncol=5,name="I") C <- mxAlgebra(L %*% Rff %*% t(L),name="C") U <- mxAlgebra(I-(I*C),name="U") SD <- mxMatrix(type="Full",nrow=5,ncol=1,free=T,values=0.6, labels=c("sd1","sd2","sd3","sd4","sd5"),lbound=0,name="SD") Sigma <- mxAlgebra( (C+U) * (SD%*%t(SD)), name="Sigma", dimnames=list(colnames(demoOneFactor),colnames(demoOneFactor))) factorModelStand <- mxModel( "1FactorStd", Rff, L, I, C, U, SD, Sigma, mxMLObjective(covariance="Sigma"), # mxCI("L"), mxData(observed=cov(demoOneFactor), type="cov", numObs=500) ) run_fms <- mxRun(factorModelStand,intervals=T) summary(run_fms) factorModelStand <- mxModel( "1FactorStd", Rff, L, I, C, U, SD, Sigma, mxMLObjective(covariance="Sigma"), mxCI("L"), mxData(observed=cov(demoOneFactor), type="cov", numObs=500) ) run_fms <- mxRun(factorModelStand,intervals=T) summary(run_fms)
Comments
#1
This does seem wrong, and is still true with current beta
#2
Verify that SVN 3209 fixes it and close the bug.
#3
This works.
#4
Thanks, Joshua!