Jump to content


Photo

Output PRED to table


  • Please log in to reply
5 replies to this topic

#1 Thomas Klitgaard

Thomas Klitgaard

    Member

  • Members
  • PipPip
  • 15 posts

Posted 04 April 2022 - 01:31 PM

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 bwendt@certara.com

bwendt@certara.com

    Advanced Member

  • Administrators
  • 282 posts

Posted 04 April 2022 - 01:51 PM

Hi Thomas,

 

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

 

PRED.png ]

 

Hope, this helps.

 

Bernd



#3 smouksassi1

smouksassi1

    Advanced Member

  • Members
  • PipPipPip
  • 231 posts
  • LocationMontreal

Posted 04 April 2022 - 02:12 PM

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=TVSM*exp(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 =   TVLS*exp(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=-KA*A1)
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+SMAX*CONC1/(SC50+CONC1))
STIM2=HB0/LS*(1+SMAX*CONCdel/(SC50+CONCdel))
STIM3 =(HB0-hbbas)/LS*AA*(1-COR)
 
STIM1_PRED=HB0_noeta/LS_noeta*(1+SMAX_noeta*CONC1/(SC50_noeta+CONC1))
STIM2_PRED=HB0_noeta/LS_noeta*(1+SMAX_noeta*CONCdel/(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

  • bwendt@certara.com likes this

#4 Thomas Klitgaard

Thomas Klitgaard

    Member

  • Members
  • PipPip
  • 15 posts

Posted 04 April 2022 - 02:34 PM

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 bwendt@certara.com

bwendt@certara.com

    Advanced Member

  • Administrators
  • 282 posts

Posted 04 April 2022 - 02:45 PM

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.c...iction&rhsyns=

 

Bernd



#6 smouksassi1

smouksassi1

    Advanced Member

  • Members
  • PipPipPip
  • 231 posts
  • LocationMontreal

Posted 05 April 2022 - 07:33 AM

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)





0 user(s) are reading this topic

0 members, 0 guests, 0 anonymous users