Hi, guys
I’m trying to use SAEM to evaluate my pre-clinical data. I tried to do my analysis using FOCE, but it did not work. I thought it might be because of the high variability in my profiles and the small number of samples.
So, I used this valuable topic of discussion as a reference to adapt my model (SAEM algorithm parameters - #8 by bobbrown) beyond the standard documentation of the Pumas of course.
I’m working with the following code:
mdl_2cmp_prop_saem = @ emmodel begin
@ param begin
tvcl ~ 1 | LogNormal #when FC == 1
tvv ~ 1 | LogNormal
tvq ~ 1 | LogNormal
tvvp ~ 1 | LogNormal
beta_Cl_FC_2 ~ 1 | Normal # when FC == 2
beta_V1_FC_2 ~ 1 | Normal
beta_Q_FC_2 ~ 1 | Normal
end
@ random begin # Previously I found IVV only in Cl, Vc, and Q, but not in Vp (Monolix)
η1 ~ 1 | Normal
η2 ~ 1 | Normal
η3 ~ 1 | Normal
end
@ covariance (3)
@ covariates FC
@ pre begin #here I'm evaluating my categorical COV (1 or 2, where 1 is the standard)
CL = tvcl * beta_Cl_FC_2 * (FC ==2)* exp(η1)
Vc = tvv * beta_V1_FC_2 * (FC ==2)* exp(η2)
Q = tvq * beta_Q_FC_2 * (FC ==2)* exp(η3)
Vp = tvvp
end
@dynamics
#Central1Periph1 - I've tried using it as a library but it was the same
Central' = -(CL+Q)/Vc*Central + Q/Vp*Peripheral
Peripheral' = Q/Vc*Central - Q/Vp*Peripheral
@post begin
cp = @. Central / Vc
end
@error begin
DV ~ Normal(cp)
end
end
params_2cmp_prop_saem = (
tvcl = 7119.54,
tvv = 1795.2,
tvq = 3471.5,
Ω = (0.1, 0.1, 0.1),
tvvp =7124.9,
σ = 1,
η1 = 0.1,
η2 = 0.1,
η3 = 0.1,
beta_Cl_FC_2 = 1,
beta_V1_FC_2 = 1,
beta_Q_FC_2 = 1,
)
SAEM_fit = fit(mdl_2cmp_prop_saem, pop1, params_2cmp_prop_saem, Pumas.SAEM())
And after countless attempts at correction, I have this error:
ERROR: TaskFailedException
nested task error: UndefVarError: `Central` not defined
Stacktrace:
[1] macro expansion
@ c:\Users\Ana\.julia\Curcumina\popPK IV\popPK iv 2cpt_monolix to pumas.jl:274 [inlined]
[2] (::var"#18#20")(t::Float64, param::NamedTuple{(:tvcl, :tvv, :tvq, :tvvp, :beta_Cl_FC_2, :beta_V1_FC_2, :beta_Q_FC_2, :η1, :η2, :η3, :FC, :CL, :Vc, :Q, :Vp), Tuple{Float64, Float64, Float64, Float64, Float64, Float64, Float64, Float64, Float64, Float64, Int64, Vararg{Float64, 4}}})
@ Main C:\a\PumasSystemImages\PumasSystemImages\julia_depot\packages\Pumas\Zl9W3\src\dsl\emmodel_macro.jl:38
[3] _calc_expectations_loop!
I would appreciate it if you could suggest something!