as.parframe.parlist {dMod} | R Documentation |
Coerce object to a parameter frame
Description
Coerce object to a parameter frame
Usage
## S3 method for class 'parlist'
as.parframe(x, sort.by = "value", ...)
as.parframe(x, ...)
Arguments
x |
object to be coerced |
sort.by |
character indicating by which colum the returned parameter frame
should be sorted. Defaults to |
... |
other arguments |
Value
object of class parframe.
Examples
## Generate a prediction function
regfn <- c(y = "sin(a*time)")
g <- Y(regfn, parameters = "a")
x <- Xt(condition = "C1")
## Generate data
data <- datalist(
C1 = data.frame(
name = "y",
time = 1:5,
value = sin(1:5) + rnorm(5, 0, .1),
sigma = .1
)
)
## Initialize parameters and time
pars <- c(a = 1)
times <- seq(0, 5, .1)
plot((g*x)(times, pars), data)
## Do many fits from random positions and store them into parlist
out <- as.parlist(lapply(1:50, function(i) {
trust(normL2(data, g*x), pars + rnorm(length(pars), 0, 1), rinit = 1, rmax = 10)
}))
summary(out)
## Reduce parlist to parframe
parframe <- as.parframe(out)
plotValues(parframe)
## Reduce parframe to best fit
bestfit <- as.parvec(parframe)
plot((g*x)(times, bestfit), data)
[Package dMod version 1.0.2 Index]