1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
#
# Copyright 2007-2018 by the individuals mentioned in the source code history
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
# a simple one factor ordinal model
require(OpenMx)
data(myFADataRaw)
oneFactorOrd <- myFADataRaw[,c("z1", "z2", "z3")]
oneFactorOrd$z1 <- mxFactor(oneFactorOrd$z1, levels=c(0, 1))
oneFactorOrd$z2 <- mxFactor(oneFactorOrd$z2, levels=c(0, 1))
oneFactorOrd$z3 <- mxFactor(oneFactorOrd$z3, levels=c(0, 1, 2))
oneFactorModel <- mxModel("Common Factor Model Path Specification",
type="RAM",
mxData(
observed=oneFactorOrd,
type="raw"
),
manifestVars=c("z1","z2","z3"),
latentVars="F1",
# residual variances
mxPath(
from=c("z1","z2","z3"),
arrows=2,
free=FALSE,
values=c(1,1,1),
labels=c("e1","e2","e3")
),
# latent variance
mxPath(
from="F1",
arrows=2,
free=TRUE,
values=1,
labels ="varF1"
),
# factor loadings
mxPath(
from="F1",
to=c("z1","z2","z3"),
arrows=1,
free=c(FALSE,TRUE,TRUE),
values=c(1,1,1),
labels=c("l1","l2","l3")
),
# means
mxPath(
from="one",
to=c("z1","z2","z3","F1"),
arrows=1,
free=FALSE,
values=0,
labels=c("meanz1","meanz2","meanz3","meanF")
),
# thresholds
mxThreshold(vars=c("z1", "z2", "z3"),
nThresh=c(1,1,2),
free=TRUE,
values=c(-1, 0, -.5, 1.2)
)
) # close model
oneFactorCon <- omxConstrainMLThresholds(oneFactorModel)
oneFactorResults <- mxRun(oneFactorCon)
|
