Learn R Programming

REMixed (version 1.1.2)

plotSAEM: Display the value of parameters at each iteration

Description

Display the value of parameters at each iteration

Usage

plotSAEM(fit, paramToPlot = "all", trueValue = NULL)

Value

For each parameters, the values at the end of each iteration of remix algorithm is drawn. Moreover, the SAEM steps of each iteration are displayed.

Arguments

fit

object of class remix, from remix or a certain build from cv.remix output.

paramToPlot

Population parameters to plot (which have been estimated by SAEM) ;

trueValue

(for simulation purpose) vector named of true values ;

See Also

remix, cv.remix.

Examples

Run this code
if (FALSE) {
project <- getMLXdir()

ObsModel.transfo = list(S=list(AB=log10),
                        linkS="yAB",
                        R=rep(list(S=function(x){x}),5),
                        linkR = paste0("yG",1:5))

alpha=list(alpha0=NULL,
           alpha1=setNames(paste0("alpha_1",1:5),paste0("yG",1:5)))

y = c(S=5,AB=1000)
lambda = 1440

res = remix(project = project,
            dynFUN = dynFUN_demo,
            y = y,
            ObsModel.transfo = ObsModel.transfo,
            alpha = alpha,
            selfInit = TRUE,
            eps1=10**(-2),
            eps2=1,
            lambda=lambda)

plotConvergence(res)

trueValue = read.csv(paste0(dirname(project),"/demoSMLX/Simulation/populationParameters.txt"))

plotSAEM(res,paramToPlot = c("delta_S_pop","phi_S_pop","delta_AB_pop"),trueValue=trueValue)
}

Run the code above in your browser using DataLab