inst/models/values-1.1.1.3-U046_solve_saem-unix.R

In nlmixrdevelopment/nlmixr.examples: Nonlinear Mixed Effects Models in Population Pharmacokinetics and Pharmacodynamics

expected_values[[runno]] <- list(lik = c(-Inf, Inf, Inf), param = c(1.587, 4.471, 1.6025, 
4.0155, NaN), stdev_param = c(0.39167, 0.40172, 0.37913, 0.39653, 
NA), sigma = c(prop.err = NaN), parFixedDf = structure(list(Estimate = c(lCl = 1.58695532203411, 
lVc = 4.47097852854617, lQ = 1.60249299439807, lVp = 4.01547659835727, 
prop.err = NaN), SE = c(lCl = 0.0394241542239158, lVc = 0.0396078530596976, 
lQ = 0.0403902440147608, lVp = 0.039850526304257, prop.err = NA
), "%RSE" = c(2.4842636510638, 0.885887794065898, 2.52046306323681, 
0.992423323312603, NA), "Back-transformed" = c(4.8888412979973, 
87.4422458970768, 4.96539571088139, 55.4497163496269, NaN), "CI Lower" = c(4.52530690714031, 
80.910901117849, 4.58747411474492, 51.2835994399591, NA), "CI Upper" = c(5.28157973093313, 
94.5008183308702, 5.37445093943821, 59.9542753790866, NA), "BSV(CV%)" = c(41.8484980140224, 
40.7185948528273, 41.2645311112213, 39.3170647033797, NA), "Shrink(SD)%" = c(43.0218348272827, 
46.1938784891807, 43.1054388412945, 42.0587490874686, NA)), class = "data.frame", row.names = c("lCl", 
"lVc", "lQ", "lVp", "prop.err")), omega = structure(c(0.153407886892481, 
0, 0, 0, 0, 0.161378506274859, 0, 0, 0, 0, 0.143739376399012, 
0, 0, 0, 0, 0.157239748979002), .Dim = c(4L, 4L), .Dimnames = list(
    c("eta.Vc", "eta.Cl", "eta.Vp", "eta.Q"), c("eta.Vc", "eta.Cl", 
    "eta.Vp", "eta.Q"))), time = structure(list(saem = 20.989, 
    setup = 4.04465400000003, table = 0.0139999999999532, covariance = 0.00999999999999091, 
    other = 0.561346000000025, logLik = 0.286999999999978), class = "data.frame", row.names = "elapsed"), 
    objDf = structure(list(OBJF = Inf, AIC = Inf, BIC = Inf, 
        "Log-likelihood" = -Inf, "Condition Number" = 1.68010206242144), row.names = "gauss3_1.6", class = "data.frame"))