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
|
#
# Copyright 2007-2021 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.
setClass(Class = "MxFitFunctionR",
representation = representation(
fitfun = "function",
units = "character",
model = "MxModel",
flatModel = "MxFlatModel",
state = "list"),
contains = "MxBaseFitFunction")
setMethod("initialize", "MxFitFunctionR",
function(.Object, ...) {
.Object <- callNextMethod()
.Object@fitfun <- ..1
.Object@state <- ..2
.Object@expectation <- as.integer(NA)
.Object@units <- ..3
.Object
})
setMethod("genericFitFunConvert", signature("MxFitFunctionR"),
function(.Object, flatModel, model, labelsData, dependencies) {
.Object <- callNextMethod()
.Object@model <- model
.Object@flatModel <- flatModel
.Object
})
mxFitFunctionR <- function(fitfun, ..., units="-2lnL") {
if (!is.function(fitfun)) {
stop("First argument 'fitfun' must be of type function")
}
if (length(formals(fitfun)) != 2) {
stop("The fit function must take exactly two arguments: a model and a persistent state")
}
state <- list(...)
return(new("MxFitFunctionR", fitfun, state, units))
}
displayRFitFun <- function(fitfunction) {
cat("MxFitFunctionR", omxQuotes(fitfunction@name), '\n')
cat("$fitfun (fitness function) \n")
cat("$units: ", omxQuotes(fitfunction@units), '\n')
print(fitfunction@fitfun)
if (length(fitfunction@result) == 0) {
cat("$result: (not yet computed) ")
} else {
cat("$result:\n")
}
print(fitfunction@result)
invisible(fitfunction)
}
setMethod("print", "MxFitFunctionR", function(x,...) {
displayRFitFun(x)
})
setMethod("show", "MxFitFunctionR", function(object) {
displayRFitFun(object)
})
|
