Output PRED to table

Phoenix Certara Forums Modelling and Simulation
1

Hi everyone,

Is there a way to refer to the population prediction in PHX modelling object…Run Options…Table01…Variables, so that I can get it into Table01 for indicated Times, just like the individual prediction (“C”).

“PRED” doesn’t seem to generate an output.

Thx

BR Thomas Klitgaard

2

Hi Thomas,

you can use pred() function to generate the output:

PRED.png
PRED.png696×475 5.44 KB

Hope, this helps.

Bernd

3

if you had an old version that does not support pred you can always edit your model to have PRED where etas are zeroed:

test(){
covariate(cov1,cov2,…)
fcovariate(cov3,…)
fixef(TVSM(freeze)=c(,xxx,))
stparm(SMAX=TVSMexp(ETASM))
fixef(TVSC(freeze) = c(,xxx,) )
stparm(SC50 = TVSC ((EPOCOR/6000)^coefepoSC1)exp(ETASC))
fixef(coefepoSC1 = c(,xxx,))
fixef(TVLS=c(,x,))
stparm(LS=TVLSexp(dLSdDialysis1(dialysis==1)) * exp(dLSdDialysis2(dialysis==2))exp(ETALS))
fixef(dLSdDialysis1= c(, -xxx, ))
fixef(dLSdDialysis2 = c(, -xxx, ))
fixef(TVHB0=c(,9.44287796378036,))
stparm(HB0=TVHB0 exp(dHBdDialysis1*(dialysis==1)) * exp(dHBdDialysis2*(dialysis==2))*(exp(ETAHB0)))
fixef(dHBdDialysis1= c(, -xxx, ))
fixef(dHBdDialysis2 = c(, -xxx, ))

ranef(block(ETASM,ETASC,ETALS,ETAHB0) = c(1.6371413,
2.242861,4.388878,
0.15415201,-0.15082374,0.14648225,
-0.20766072,-0.20145853,-0.040917645,0.040534749))

SMAX_noeta = TVSM
SC50_noeta = TVSC*((EPOCOR/6000)^coefepoSC1)
LS_noeta = TVLSexp(dLSdDialysis1(dialysis==1)) * exp(dLSdDialysis2*(dialysis==2))
HB0_noeta = TVHB0* exp(dHBdDialysis1*(dialysis==1)) * exp(dHBdDialysis2*(dialysis==2))

dosepoint(A1,bioavail= pedi == 0 ? xxx: yyy)
dosepoint(A2)
sequence { A3 = hbbas }
fcovariate(CL,VD,KA)

deriv(A1=-KAA1)
deriv(A2=KAA1-(CL/VD)*A2)
CONC1=A2/VD
A2del=delay(A2,LS)# drug amount for t<TR is 0
CONCdel=A2del/VD

AA=(t>LS?0:1)
sequence { A3 = hbbas }
sequence { A3_PRED = hbbas }

deriv(A3=STIM1-STIM2+STIM3)
deriv(A3_PRED=STIM1_PRED-STIM2_PRED+STIM3_PRED)
Hb_pred = A3_PRED
STIM1=HB0/LS*(1+SMAXCONC1/(SC50+CONC1))
STIM2=HB0/LS(1+SMAXCONCdel/(SC50+CONCdel))
STIM3 =(HB0-hbbas)/LSAA*(1-COR)

STIM1_PRED=HB0_noeta/LS_noeta*(1+SMAX_noetaCONC1/(SC50_noeta+CONC1))
STIM2_PRED=HB0_noeta/LS_noeta(1+SMAX_noetaCONCdel/(SC50_noeta+CONCdel))
STIM3_PRED =(HB0_noeta-hbbas)/LS_noetaAA*(1-COR)

error(EEps=0.76242779)
observe(EObs=A3+EEps)
}

so above got the normal HB diff eq and the one for the pred version where I zero the etas ( I had to remove the parameter values

4

Hi again,

Thank you for your input, smouksassi1 and Bernd. I wanted to dump IPRED, DV, PRED in my Table01 in one run - I guess the suggest approach with setting ETAs to zero will not help here?

Im using 8.1.1.279 - the pred(CObs) approach doesnt seem to work on that version no. Further, I couldn’t see it referenced in the PHX ML documentation for 8.3 ?

BR

5

Samer might know better when this function was introduced, but here is the link to the part of online help from Phoenix 8.3:

https://onlinehelp.certara.com/phoenix/8.3/index.html#t=topics%2FSimple_run_mode.htm&rhsearch=population%20prediction&rhhlterm=population%20prediction&rhsyns=%20

Bernd

6

Yes I do it in one go here is an example of my table statement: the GUI builds this for you automatically depending on what options you input in the field while in RsNLME you need to type it in:

below I ask for the IPRED, PRED and the observed DV in one go if you can copy paste your model I can help you recode it to give you the PRED or if possible at your company upgrade to the latest version which had several enhancements from 8.1

table(file=“table01.csv”,time(seq(0,400,1)), covr(CL,VD,KA,study,dialysis), dose(A1,A2), obs(EObs),A3,A3_PRED,SMAX,SC50,LS,HB0,SMAX_noeta,SC50_noeta,LS_noeta,HB0_noeta, specvar(TAD))
table(file=“posthoc.csv”, time(0), SMAX,SC50,LS,HB0)