Thanks for the response. In my project, my PD response are in 6 ordered categories (2,1,0,-1,-2,-3), 2 is the lowest (correlates with low concentration) and -3 is the highest (correlates with high concentration). they are sedation scores, negative values represent sedated and calm while positive values are agitated and awake.
I have a PK model and I fixed the parameters, now I am adding the PD data. I attempted to write the code for the logistic model as below:
stparm(th1=tvth1+nth1)
stparm(th2 = tvth2)
stparm(th3 = tvth3)
stparm(th4 = tvth4)
stparm(th5 = tvth5)
stparm(th6 = tvth6)
fixef(tvth1= c(, 1, ))
fixef(tvth2 = c(0, 1, ))
fixef(tvth3=c(0,1,))
fixef(tvth4 = c(0,1,))
fixef(tvth5 = c(0,1,))
fixef(tvth6 = c(0,1,))
#ilogit
Prob(-3) = ilogit(th1+th2*C)
Prob(-3 or -2) = ilogit(th1+th2*C+th3)
Prob(-3 or -2 or -1) = ilogit(th1+th2*C+th3+th4)
Prob(-3 or -2 or -1 or 0) = ilogit(th1+th2*C+th3+th4+th5)
Prob(-3 or -2 or -1 or 0 or 1) = ilogit(th1+th2*C+th3+th4+th5+th6)
multi(EObs,ilogit,(th1+th2*C+th3+th4+th5+th6),(th1+th2*C+th3+th4+th5),(th1+th2*C+th3+th4),(th1+th2*C+th3),(th1+th2*C))
model wouldn't run, i'm not sure if my code is correct, could you help?
thank you!