Adding arguments to mxRun() for checkpointing and socket communication
[openmx:openmx.git] / man / mxCompare.Rd
1 \name{mxCompare}
2 \alias{mxCompare}
3
4 \title{Assign Model Parameters}
5
6 \description{
7     Compare the fit of one or more models to a base model.
8     The output is a table with one row per model comparison, and the following columns:
9 \describe{
10 \item{base}{Name of the base model}
11 \item{comparison}{Name of the comparison model}
12 \item{ep}{Estimated parameters of the comparison model}
13 \item{minus2LL}{Minus 2*log-likelihood of the comparison model}
14 \item{df}{Degrees in freedom of the comparison model}
15 \item{AIC}{Akaike's Information Criterion for the comparison model}
16 \item{diffLL}{Change in minus 2*log-likelihood}
17 \item{diffdf}{Change in degrees of freedom}
18 \item{p}{Significance level of the change in fitness function}}
19 }
20
21 \usage{
22 mxCompare(base, comparison, digits = 3)
23 }
24
25 \arguments{
26    \item{base}{an MxModel object or list of MxModel objects.}
27    \item{comparison}{a MxModel object or list of MxModel objects.}
28    \item{digits}{a numeric setting table decimal precision.} 
29 }
30
31 \seealso{
32 \code{\link{mxModel}} 
33 }
34
35 \examples{
36
37 data(demoOneFactor)
38 manifests <- names(demoOneFactor)
39 latents <- c("G1")
40 model1 <- mxModel("One Factor", type="RAM",
41       manifestVars = manifests,
42       latentVars = latents,
43       mxPath(from=latents, to=manifests),
44       mxPath(from=manifests, arrows=2),
45       mxPath(from=latents, arrows=2, free=FALSE, values=1.0),
46       mxData(cov(demoOneFactor), type="cov",numObs=500)
47 )
48
49 fit1 <- mxRun(model1)
50
51 latents <- c("G1","G2")
52 model2 <- mxModel("Two Factor", type="RAM",
53       manifestVars = manifests,
54       latentVars = latents,
55       mxPath(from=latents[1], to=manifests[1:3]),
56       mxPath(from=latents[2], to=manifests[4:5]),
57       mxPath(from=manifests, arrows=2),
58       mxPath(from=latents, arrows=2, free=FALSE, values=1.0),
59       mxData(cov(demoOneFactor), type="cov",numObs=500)
60 )
61 fit2 <- mxRun(model2)
62
63 mxCompare(fit1, c(fit1,fit2), digits=3)
64 }