Output PRED to table

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

Hi Thomas,

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

Hope, this helps.

Bernd

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=TVLS
exp(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=KA
A1-(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)/LS
AA*(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_noeta
AA*(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

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

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

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)