dMod (version 0.3.1)

parlist: Parameter list

Description

The special use of a parameter list is to save the outcome of multiple optimization runs provided by mstrust, into one list.

Parameter list

Fitlists carry an fit index which must be held unique on merging multiple fitlists.

Usage

parlist(...)
as.parlist(x = NULL)
"summary"(object, ...)
"c"(...)

Arguments

...

Objects to be coerced to parameter list.

x

list of lists, as returned by trust

object

a parlist

See Also

load.parlist, plot.parlist

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)

## Plot array with all solutions
plotArray(parframe, g*x, times, data)