Closed joaoabr closed 3 years ago
Assmblerr version 0.0.0.9000. When generating a 2-compartment model with first-order absorption, numbers assigned to VC and VP do no seem to be correct, as there a depot compartment is present.
For instance:
pkm <- pk_model() + pk_absorption_fo() + pk_distribution_2cmp() + pk_elimination_linear() + obs_additive(conc~C["central"])
Leads to the following, where V1 (=VC) and V2 (=VP) should be V2 and V3, respectively:
$PROBLEM $INPUT ID TIME DV AMT $DATA data.csv IGNORE=@ $SUBROUTINES ADVAN4 TRANS4 $MODEL COMP=(DEPOT) COMP=(CENTRAL) COMP=(PERIPHERAL) $PK MAT = THETA(1) * EXP(ETA(1)) VC = THETA(2) * EXP(ETA(2)) VP = THETA(3) * EXP(ETA(3)) Q = THETA(4) * EXP(ETA(4)) CL = THETA(5) * EXP(ETA(5)) KA = 1/MAT V1 = VC V2 = VP $ERROR CONC = A(2)/VC Y = CONC + EPS(1) $ESTIMATION METHOD=COND INTERACTION $THETA (0, 1, Inf) ; POP_MAT $THETA (0, 1, Inf) ; POP_VC $THETA (0, 1, Inf) ; POP_VP $THETA (0, 1, Inf) ; POP_Q $THETA (0, 1, Inf) ; POP_CL $OMEGA 0.1; IIV_MAT $OMEGA 0.1; IIV_VC $OMEGA 0.1; IIV_VP $OMEGA 0.1; IIV_Q $OMEGA 0.1; IIV_CL $SIGMA 0.1; RUV_ADD
Assmblerr version 0.0.0.9000. When generating a 2-compartment model with first-order absorption, numbers assigned to VC and VP do no seem to be correct, as there a depot compartment is present.
For instance:
Leads to the following, where V1 (=VC) and V2 (=VP) should be V2 and V3, respectively: