data(bcrp)
formula1<- I(cesdt1-cesdt3)~cond | nationality+marital+wcht1+
age+trext+comorbid+disopt1+uncomt1+negsoct1
control1<-quint.control(maxl=5,Bootstrap=FALSE)
quint1<-quint(formula1, data= subset(bcrp,cond<3),control=control1 )
summary(quint1)
Run the code above in your browser using DataLab