Hello.
I'd like to fit single dose transdermal PK profile using Phoenix model.
Basically, I used "Transdermal example of nicotine patch" text model with some modification
test(){
deriv(A1 = Frate + Srate - Ke * A1 - K12 * A1 + K21 * A2)
deriv(A2 = K12 * A1 - K21 * A2)
Dose = 10000
Frate = (t <= Ftau ? ((Dose-SDose)/Ftau) : 0)
Srate = (t <= Stau ? (SDose/Stau) : 0)
C = A1 / V
error(CEps = 1)
observe(CObs = C * (1 + CEps))
stparm(Ftau = tvFtau * exp(nFtau))
stparm(Stau = tvStau * exp(nStau))
stparm(SDose = tvSDose * exp(nSDose))
stparm(V = tvV * exp(nV))
stparm(Ke = tvKe * exp(nKe))
stparm(K12 = tvK12 * exp(nK12))
stparm(K21 = tvK21 * exp(nK21))
stparm(Tlag = tvTlag * exp(nTlag))
fixef(tvV = c(0, 384, ))
fixef(tvKe = c(0, 0.15, ))
fixef(tvK12 = c(0, 0.01, ))
fixef(tvK21 = c(0, 0.9, ))
fixef(tvTlag = c(0, 0.75, ))
fixef(tvFtau = c(8, 12, 16))
fixef(tvStau = c(16, 24, 28))
fixef(tvSDose = c(0, 8000, ))
ranef(diag(nV, nKe, nK12, nK21, nTlag, nFtau, nStau, nSDose) = c(1, 1, 1, 1, 1, 1, 1, 1))
}
But, Our PK profile has "Tlag" and I'd like to enter the Tlag as a parameter in above model.
How can I put the "Tlag" in the model?
Can you modify this model to Phoenix graphical model?
and How can I simulate multiple dose profile (Tau = 24hr) using this model?
Please give me advices
Thanks [file name=TDS_modeling.phxproj size=365000]http://www.pharsight.com/extranet/media/kunena/attachments/legacy/files/TDS_modeling.phxproj[/file]