R: Estimation of lower and upper bounds for missing values.

estim.bound {imp4p} | R Documentation

## Estimation of lower and upper bounds for missing values.

### Description

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`.

### Usage

    estim.bound(tab, conditions, q=0.95)

### Arguments

`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).

### 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.

### Value

A list composed of:

`tab.lower` | A matrix with the lower bounds for each missing value in `tab`.
`tab.upper` | A matrix with the upper bounds for each missing value in `tab`.

### Author(s)

Quentin Giai Gianetto <quentin2g@yahoo.fr>

### See Also

`prob.mcar`,`prob.mcar.tab`

### Examples

    #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);

[Package imp4p version 1.2 Index]