R: Estimation of a vector of probabilities that missing values...

prob.mcar {imp4p}

R Documentation

Estimation of a vector of probabilities that missing values are MCAR.

Description

This function returns a vector of probabilities that each missing value is MCAR from specified confidence intervals.

Usage

prob.mcar(b.l,b.u,absc,pi.mcar,F.tot,F.na)

Arguments

`b.l`	A numeric vector of lower bounds for missing values.
`b.u`	A numeric vector of upper bounds for missing values.
`absc`	The interval on which is estimated the MCAR data mechanism.
`pi.mcar`	The estimated proportion of MCAR values.
`F.tot`	An estimation of the cumulative distribution function of the complete values on the interval `absc`.
`F.na`	An estimation of the cumulative distribution function of the missing values on the interval `absc`.

Value

A numeric vector of estimated probabilities to be MCAR for missing values in the confidence intervals defined by b.l and b.u. The input arguments absc, pi.mcar, F.tot and F.na can be estimated thanks to the function estim.mix.

Author(s)

Quentin Giai Gianetto <quentin2g@yahoo.fr>

Examples


#Simulating data
#Simulating data
res.sim=sim.data(nb.pept=2000,nb.miss=600,pi.mcar=0.2,para=10,nb.cond=2,nb.repbio=3,
nb.sample=5,m.c=25,sd.c=2,sd.rb=0.5,sd.r=0.2);

#Imputation of missing values with the slsa algorithm
dat.slsa=impute.slsa(tab=res.sim$dat.obs,conditions=res.sim$condition,repbio=res.sim$repbio);

#Estimation of the mixture model
res=estim.mix(tab=res.sim$dat.obs, tab.imp=dat.slsa, conditions=res.sim$condition);

#Computing probabilities to be MCAR
born=estim.bound(tab=res.sim$dat.obs,conditions=res.sim$condition);

#Computing probabilities to be MCAR in the first column of result$tab.mod
proba=prob.mcar(b.l=born$tab.lower[,1],b.u=born$tab.upper[,1],absc=res$abs.mod,
pi.mcar=res$pi.mcar[1], F.tot=res$F.tot[,1], F.na=res$F.na[,1]);

[Package imp4p version 0.7 Index]