# NOT RUN {
library("mixor")
data("schizophrenia")
### Random intercept
SCHIZO1.fit<-mixor(imps79o ~ TxDrug + SqrtWeek + TxSWeek, data=schizophrenia,
id=id, link="probit")
print(SCHIZO1.fit)
### Random intercept and slope
SCHIZO2.fit<-mixor(imps79o ~ TxDrug + SqrtWeek + TxSWeek, data=schizophrenia,
id=id, which.random.slope=2, link="probit")
print(SCHIZO2.fit)
# }
Run the code above in your browser using DataLab