cond.prob: Calculate conditional probabilities for observing records at...

In ELCIC: The Empirical Likelihood-Based Consistent Information Criterion

Calculate conditional probabilities for observing records at each time point

Description

A function calculates conditional probabilities for longitudinal missing data. The observing probability is at observation-level.

Usage

cond.prob(x_mis,gamma,id,time)

Arguments

x_mis	A matrix containing covariates for the missing data model. The first column should be all ones corresponding to the intercept.
gamma	coefficients calculated from missing data model
id	A vector indicating subject id.
time	The number of observations in total for each subject

Value

a vector containing conditional probabilities.

Examples

## tests
# load data
data(wgeesimdata)
library(wgeesel)
data_wgee<-data.frame(do.call(cbind,wgeesimdata))
corstr<-"exchangeable"
dist<-"binomial"
id<-data_wgee$id
# obtain the estimates
fit<-wgee(y~x1+x2+x3,data_wgee,id,family=dist,corstr =corstr,scale = NULL,
         mismodel =obs_ind~x_mis1)
beta<-as.vector(summary(fit)$beta)
rho<-summary(fit)$corr
phi<-summary(fit)$phi
#calculate observing probabilies for all observations
gamma<-as.vector(summary(fit$mis_fit)$coefficients[,1])
x_mis<-wgeesimdata$x_mis
pi<-cond.prob(x_mis,gamma,id,time=3)

ELCIC documentation built on Feb. 16, 2023, 7:18 p.m.