LPEseq.predict.var: A function predicting per-gene variance using LPEseq.var()...

In JungsooGIM/LPEseq: LPEseq: local-pooled-error test for RNA-seq data with a small number of replicates

Description

A function predicting per-gene variance using LPEseq.var() output

Usage

LPEseq.predict.var(gene_expr, var.spline)

Arguments

gene_expr	gene expression whose variance to be estimated
var.spline	output of LPEseq.var() function

Details

See the reference paper

Value

numeric

estimated variance

Author(s)

Jungsoo Gim

References

J. Gim and T. Park, PLoS ONE under revision (http://bibs.snu.ac.kr/software/LPEseq)

Examples

## Generating a non-replicated simulation dataset
simData.norep <- generateData(n.deg=1000, n.rep=1)

## Normalization
normData.norep <- LPEseq.normalise(simData.norep[,-3])

## Generating LPE object
normData.var <- LPEseq.var(normData.norep, n.bin=100, df=10, d=1.2)

## Predicting variance
LPEseq.predict.var(8.423, normData.var)

JungsooGIM/LPEseq documentation built on May 16, 2019, 10:04 a.m.