Updated copyright to 2013 for R/ demo/ models/passing and src/ folders, and also...
[openmx:openmx.git] / models / passing / FIMLLikelihoodFrontend.R
1 #
2 #   Copyright 2007-2013 The OpenMx Project
3 #
4 #   Licensed under the Apache License, Version 2.0 (the "License");
5 #   you may not use this file except in compliance with the License.
6 #   You may obtain a copy of the License at
7
8 #        http://www.apache.org/licenses/LICENSE-2.0
9
10 #   Unless required by applicable law or agreed to in writing, software
11 #   distributed under the License is distributed on an "AS IS" BASIS,
12 #   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 #   See the License for the specific language governing permissions and
14 #   limitations under the License.
15
16 require(OpenMx)
17
18 model <- mxModel('model')
19 obsData <- matrix(c(1:10), 10, 1)
20 data <- mxData(obsData, "raw")
21 chol <- mxMatrix("Full", 2, 2, c(T,T,F,T), c(1,.2,0,1), name="Chol") 
22 expCov <- mxAlgebra(Chol %*% t(Chol), name="expCov") 
23 expMean <- mxMatrix("Full", 1, 2, T, c(0,0), "expMean")
24 objective <- mxFIMLObjective("expCov", "expMean")
25 foo <- mxMatrix('Full', 1, 1, name = 'foo')
26 bar <- mxMatrix('Full', 10, 1, name = 'bar')
27 model <- mxModel(model, data, chol, expCov, expMean, objective, foo, bar)
28 omxCheckIdentical(mxEval(objective + foo, model, TRUE), matrix(as.numeric(NA), 1, 1))
29 omxCheckError(mxEval(objective + bar, model, TRUE), "The following error occurred while evaluating the expression 'model.fitfunction + model.bar' in model 'model' : non-conformable arrays")
30 objective <- mxFIMLObjective("expCov", "expMean", vector = TRUE)
31 model <- mxModel(model, objective)
32 omxCheckError(mxEval(objective + foo, model, TRUE), "The following error occurred while evaluating the expression 'model.fitfunction + model.foo' in model 'model' : non-conformable arrays")
33 omxCheckIdentical(mxEval(objective + bar, model, TRUE), matrix(as.numeric(NA), 10, 1))