Hi, I am trying to fit a SEM model (the ram specification of the model, along with the estimated polycoric VCmatrix (CC32), is attacched), but I get the following error message: "Error in 1:m : NA/NaN argument"I did not find any help on internet, so I tried to understend what is going wrong calling: getAnywhere(sem.default) and reproducing the sem.default function step by step.I found that the problem arises from:
    n.fix <- length(fixed.x)

    if (!is.null(fixed.x)) {

        for (i in 1:n.fix) {

            for (j in 1:i) {

                ram <- rbind(ram, c(2, fixed.x[i], fixed.x[j], 

                  0, S[fixed.x[i],
fixed.x[j]]))

            }

        }

    }

    m <- max(ram[, c(2, 3)])
Indeed, if I run just:    
for (i in 1:n.fix) {for (j in 1:i) {G<-c(2, fixed.x[i], fixed.x[j],                  0, S[fixed.x[i], fixed.x[j]]) }
I get only this (G=):  "2"       "Autunno" "Autunno" "0"       "1"
But I still cannot understand what is wrong.The code I used is:
library(sem)CLB<-classifyVariables(Mod.Bus)

fisse<-as.character(CLB$exogenous)Business32<-sem(Mod.Bus,
S= as.matrix(CC32), N<- 16029, fixed.x=fisse,objective= objectiveGLS,
maxiter=100000)

Can anybody help me please?Any suggestion will be truly appreciated.Annalisa 
 		 	   		  

Model ram specific.docx (21K) <http://r.789695.n4.nabble.com/attachment/4685052/0/Model%20ram%20specific.docx>
CC32.Rdata (20K) <http://r.789695.n4.nabble.com/attachment/4685052/1/CC32.Rdata>
C32.Rdata (17K) <http://r.789695.n4.nabble.com/attachment/4685052/2/C32.Rdata>




--
View this message in context: http://r.789695.n4.nabble.com/Help-with-SEM-Error-in-1-m-NA-NaN-argument-tp4685052.html
Sent from the R help mailing list archive at Nabble.com.
	[[alternative HTML version deleted]]

