You are here

The OpenMx website will be down for maintenance from 9 AM EDT on Tuesday, September 17th, and is expected to return by the end of the day on Wednesday, September 18th. During this period, the backend will be updated and the website will get a refreshed look.

provide warning when data includes variables which are not modeled?

1 post / 0 new
tbates's picture
Offline
Joined: 07/31/2009 - 14:25
provide warning when data includes variables which are not modeled?

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: