orig_pplr {puma}R Documentation

Probability of positive log-ratio

Description

This is the original version of the pplr function as found in the pplr package. This should give exactly the same results as the pplr function. This function is only included for testing purposes and is not intended to be used. It will not be available in future versions of puma.

This function calculates the probability of positive log-ratio (PPLR) between any two specified conditions in the input data, mean and standard deviation of gene expression level for each condition.

Usage

orig_pplr(e, control, experiment)

Arguments

e a data frame containing the mean and standard deviation of gene expression levels for each condition.
control an integer denoting the control condition.
experiment an integer denoting the experiment condition.

Details

The input of 'e' should be a data frame comprising of 2*n components, where n is the number of conditions. The first 1,2,...,n components include the mean of gene expression values for conditions 1,2,...,n, and the n+1, n+2,...,2*n components contain the standard deviation of expression levels for condition 1,2,...,n.

Value

The return is a data frame. The description of the components are below.

index The original row number of genes.
cM The mean expression levels under control condition.
sM The mean expression levels under experiment condition.
cStd The standard deviation of gene expression levels under control condition.
sStd The standard deviation of gene expression levels under experiment condition.
LRM The mean log-ratio between control and experiment genes.
LRStd The standard deviation of log-ratio between control and experiment genes.
stat A statistic value which is -mean/(sqrt(2)*standard deviation).
PPLR Probability of positive log-ratio.

Author(s)

Xuejun Liu, Marta Milo, Neil D. Lawrence, Magnus Rattray

References

Liu,X., Milo,M., Lawrence,N.D. and Rattray,M. (2005) Probe-level variances improve accuracy in detecting differential gene expression, technical report available upon request.

See Also

Related method bcomb

Examples

  data(exampleE)
  data(exampleStd)
  r<-bcomb(exampleE,exampleStd,replicates=c(1,1,1,2,2,2),method="map")
  p<-orig_pplr(r,1,2)

[Package puma version 1.8.1 Index]