Estimation of lower and upper bounds for missing values.
This function allows estimating lower and upper bounds for missing values of an input matrix. It can be used before to use the functions prob.mcar and prob.mcar.tab.
estim.bound(tab, conditions, q=0.95)
tab |
A data matrix containing numeric and missing values. Each column of this matrix is assumed to correspond to the intensities measured in an experimental sample, and each row to the ones of an identified peptide. |
conditions |
A vector of factors indicating the biological condition to which each sample belongs. |
q |
A numeric value allowing to define confidence intervals for missing values (see Details). |
In each condition, this function estimates lower and upper bounds for missing values of row i by:
upper(i)=max(tab[i,]);
lower(i)=max(tab[i,])-quant_diff(q);
where quant_diff(q) corresponds to a quantile value of the differences between the maximum and the minimum of the observed values for all the peptides in the condition. As a result, if q is close to 1, quant_diff(q) represents an extrem value between the maximum and the minimum of the intensity values in a condition for a peptide.
A list composed of:
tab.lower |
A matrix with the lower bounds for each missing value in |
tab.upper |
A matrix with the upper bounds for each missing value in |
Quentin Giai Gianetto <quentin2g@yahoo.fr>
#Simulating data res.sim=sim.data(nb.pept=2000,nb.miss=600,pi.mcar=0.2,para=0.5,nb.cond=2,nb.repbio=3, nb.sample=5,m.c=25,sd.c=2,sd.rb=0.5,sd.r=0.2); data=res.sim$dat.obs; cond=res.sim$conditions; #Estimation of lower and upper bounds for each missing value res=estim.bound(data,conditions=cond);
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.