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