mbetattest: Performance of multiple beta t-test on simulated data
In MBttest: Multiple Beta t-Tests

Description Usage Arguments Details Value Author(s) References See Also Examples

This function is to peform multiple beta t-test method on real data. The result lists geneid or isoformid, gene name, the other information, t-value, p-value, adjusted p-value, adjusted alpha value, rho (ρ), and symb. All these lists are ordered by absolution of t-values.

1	mbetattest(X, na, nb, W, alpha=0.05, file)

`X`	count data of RNA reads with na replicates in condition A ans nb replicates in condition B.
`na`	number of replicate libraries in condition A.
`nb`	number of replicate libraries in condition B.
`W`	numeric parameter, called omega (ω) that is a constant determined by null simulation.
`alpha`	the probabilistic threshold. User can set alpha (α)= 0.05 or 0.01 or the other values. Defalt value is 0.05
`file`	a csv file. User needs to give file name and specify direction path. But if user uses setwd function, drive is not necessarily specified in file.

t-statistic is defined as t-statistic multiplied by (rho/omega), that is,

T=t\times \frac{ρ}{ω}

where

t=\frac{(P_A-P_B)}{√(V_A+V_B)}

ρ=√{ψ ζ}

where

ψ =max(\frac{min(X_A)}{max(X_B)+1},\frac{min(X_B)}{max(X_A)+1})

ζ=log(1+\frac{\bar{X}σ^2+1}{\bar{X}_Aσ^2_A+\bar{X}_Bσ^2_B+1})

ω is a constant as threshold estimated from null data.

return a dat list: the data ordered by abs(t) contain information cloumns, data columns, t-values, rho and symb that are used to make heatmap and MAplot.

Yuan-De Tan tanyuande@gmail.com

Yuan-De Tan Anita M. Chandler, Arindam Chaudhury, and Joel R. Neilson(2015) A Powerful Statistical Approach for Large-scale Differential Transcription Analysis. Plos One, 10.1371/journal.pone.0123658.

smbetattest.

1
2
3

data(jkttcell) 

dat<-mbetattest(X=jkttcell[1:1000,],na=3,nb=3,W=1,alpha=0.05,file="jurkat_NS_48h_tag_mbetattest.csv")

Loading required package: gplots

Attaching package: 'gplots'

The following object is masked from 'package:stats':

    lowess

Loading required package: gtools