BisRNA-package {BisRNA}R Documentation

Analysis of RNA Cytosine-5 Methylation from bisulfite sequencing experiments


BisRNA analyses non-conversion of bisulfite-treated RNA of a set of samples, in order to distinguish actual cytosine-5 methylation marks from artifacts. Calculations are done as follows:

A small non-conversion p-value points to methylation of a RNA cytosine, or another event blocking bisulfite conversion.


 ## Load data

 ## Ratio (Poisson parameter / coverage), derived p-values, 
 #  adjusted for multiple testing using either BH (here), or IHW if available.
 lambda1 <- RNAmeth.poisson.par(Bisdata1)$estimate
 BisXP1  <- RNAmeth.poisson.test(Bisdata1,lambda1,method="BH")
 lambda2 <- RNAmeth.poisson.par(Bisdata2)$estimate
 BisXP2  <- RNAmeth.poisson.test(Bisdata2,lambda2,method="BH")
 lambda3 <- RNAmeth.poisson.par(Bisdata3)$estimate
 BisXP3  <- RNAmeth.poisson.test(Bisdata3,lambda3,method="BH")

 ## Combined p-values ; median and standard error of
 ## bisulfite non-conversion ratio
 BisXP.combined <- samples.combine(BisXP1,BisXP2,BisXP3)

[Package BisRNA version 0.2.2 Index]