I'm running a factor analysis, which I have done before using EQS, with OpenMx and finding an error. The code and the output are below. I'd appreciate any diagnostics and directions!

Regards,

Carlos.

# load library

require(OpenMx)

# data

AttractivenessDataRaw <- read.table("/home/pos/carlosdenner/Desktop/logattractiveness.txt",header=TRUE)

names(AttractivenessDataRaw)

# model specification

AttractivenessFactorModel<-mxModel("Attractiveness Factor Model",

type="RAM",

mxData(

observed=AttractivenessDataRaw,

type="raw"),

manifestVars=c("pageviews","downloads","members"),

latentVars="attractiveness",

# residual variances

mxPath(from=c("pageviews","downloads","members"),

arrows=2,

free=TRUE,

values=c(1,1,1),

labels=c("e1","e2","e3")

),

# latent variance

mxPath(from="attractiveness",

arrows=2,

free=TRUE,

values=1,

labels ="varAttractiveness"

),

# factor loadings

mxPath(from="attractiveness",

to=c("pageviews","downloads","members"),

arrows=1,

free=c(FALSE,TRUE,TRUE),

values=c(1,1,1),

labels =c("l1","l2","l3")

),

# means

mxPath(from="one",

to=c("pageviews","downloads","members","attractiveness"),

arrows=1,

free=c(TRUE,TRUE,TRUE,FALSE),

values=c(1,1,1,0),

labels =c("meanpageviews","meandownloads","meanmembers",NA)

)

) # close model

AttractivenessFactorFit <- mxRun(AttractivenessFactorModel)

Output:

Running Attractiveness Factor Model

Warning message:

In model 'Attractiveness Factor Model' NPSOL returned a non-zero status code 6. model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)