# NOT RUN {
ymodel= y~ t1+t3+t6 + group + t12 + t32+t62 + perf + sever
gmodel= as.factor(g)~ t3+t6+group+yp+perf+sever
qolef$t12<-qolef$t1*qolef$group
qolef$t32<-qolef$t3*qolef$group
qolef$t62<-qolef$t6*qolef$group
qolef.isni=isnigls(ymode=ymodel, gmodel=gmodel, alldata=qolef)
summary(qolef.isni)
# }
Run the code above in your browser using DataLab