copied scripts from passing
[openmx:openmx.git] / demo / UnivariateSaturated_PathRaw.R
1 require(OpenMx)
2
3 #Simulate Data
4 set.seed(100)
5 x <- rnorm (1000, 0, 1)
6 testData <- as.matrix(x)
7 selVars <- c("X")
8 dimnames(testData) <- list(NULL, selVars)
9 summary(testData)
10 mean(testData)
11 var(testData)
12
13 #example 2: Saturated Model with Raw Data and Path-Style Input
14 univSatModel2 <- mxModel("univSat2",
15     manifestVars= selVars,
16     mxPath(
17         from=c("X"), 
18         arrows=2, 
19         free=T, 
20         values=1, 
21         lbound=.01, 
22         labels="vX"
23     ),
24     mxData(
25         observed=testData, 
26         type="raw", 
27     ),
28     type="RAM"
29     )
30 univSatFit2 <- mxRun(univSatModel2)
31 EM2 <- mxEvaluate(M, univSatFit2)
32 EC2 <- mxEvaluate(S, univSatFit2)
33 LL2 <- mxEvaluate(objective,univSatFit2);
34
35
36 #Mx answers hard-coded
37 #example Mx..2: Saturated Model with Raw Data
38 Mx.EM2 <- 0.01680516
39 Mx.EC2 <- 1.061050
40 Mx.LL2 <- 2897.135
41
42
43 #Compare OpenMx results to Mx results (LL: likelihood; EC: expected covariance, EM: expected means)
44 #2:RawPat 
45 omxCheckCloseEnough(LL2,Mx.LL2,.001)
46 omxCheckCloseEnough(EC2,Mx.EC2,.001)
47 omxCheckCloseEnough(EM2,Mx.EM2,.001)