Calculate the ISNI when the regression outcome is subject to missingness and follows linear mixed-effects models (LMMs)
isnilmm(ymfix, ymran, gmodel, predprobobs = NULL, misni = FALSE, alldata)
an object of class "formula": a two-sided linear formula description of the fixed-effects part of the model to be fitted for the outcome
an object of class "formula": an one-sided linear formula description of the random-effects part of the model to be fitted for the outcome
an object of class "formula": a symbolic description of the selection model to be fitted for the missingness indictor g
Null if using buil-in multinomial transitional logistic model to obtain predicted probabilities of being observed; otherwise user supply a vector of these probabilities for all the observations in alldata.
FALSE if using the default approach to computing ISNI with a scalar nonignorability parameter; TRUE when computing ISNI with multiple nonignorability parameters.
the name of data frame containing all the variables in the model and all the observations including those intended to be collected but became missing.
# NOT RUN {
data(qolef)
qolef$t12<-qolef$t1*qolef$group
qolef$t32<-qolef$t3*qolef$group
qolef$t62<-qolef$t6*qolef$group
ymfix=y~ t1+t3+t6 + group + t12 + t32+t62 + perf + sever
gmodel=factor(g)~ group+factor(time)+yp+perf+sever
ymran=~1|id
##Random intercept model
result=isnilmm(ymfix,ymran, gmodel, alldata=qolef)
summary(result)
# }
Run the code above in your browser using DataLab