provide warning when data includes variables which are not modeled?

No replies
tbates's picture
Offline
Joined: 07/31/2009

The code below grabs the first 6 variables from myFADataRaw as manifest, but then provides the full n=9-variable covariance matrix to the mxModel

It runs, and gives no warning that half the variables are not modelled. Would be very helpful to say
"variables x-y are included in the covariance matrix, but do not appear in your model!"

or something similar?

########## reduction ##########

require(OpenMx)
data(myFADataRaw); myData = myFADataRaw;
manifests = names(myData)[1:6] # Just the first 6
#[1] "x1" "x2" "x3" "x4" "x5" "x6"
latents = c("G")

factorModel = mxModel(
"One Factor", type="RAM",
manifestVars = manifests,
latentVars = latents,
mxPath(from = latents, to=manifests),
mxPath(from = manifests, arrows=2),
mxPath(from = latents , arrows=2, free=F, values=1.0),
mxData(cov(myData), type="cov",numObs=500)
)
summary(mxRun(factorModel))
...
name parameter estimate error estimate
1 1 1
2 1 1
3 1 1
4 1 1
5 1 1
6 1 1
7 1 1
8 1 1
9 1 1
10 1 1
11 1 1
12 1 1

Observed statistics: 45
Estimated parameters: 12
Degrees of freedom: 33
AIC:
BIC:
adjusted BIC:
RMSEA: