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

This function is to peform mBeta t-test with rho=1 and omega=1 on simulated data. The result lists differentially expressed genes or isoforms marked by symbol="+" and their rho values. The rho values are used to calculate omega value for performance of mBeta t-tests on the real data.

1 | ```
smbetattest(X, na, nb, alpha)
``` |

`X` |
simulated count data with N genes or isoforms. |

`na` |
number of replicate libraries in condition A. |

`nb` |
number of replicate libraries in condition B. |

`alpha` |
statistical probabilistic threshold, default value is 0.05. |

Before performing mbeta t-test on real data, user needs omega (w) value for the threshold of rho(*ρ*). To determine omega value, user is requred to simulate null data having the same gene or isoform number and the same numbers of replicate libraries in two conditions and then performs mbeta t-test on the simulated null data by setting rho =1 and omega =1. To calculate accurately omega value, user needs such performance on 4-6 simulated null
datasets. Manual provides method for omega calculation.

Return results from multple beta t-tests on simulated data.

Yuan-De Tan [email protected]

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.

See Also as `mbetattest`

1 2 3 | ```
data(skjt)
mysim<-smbetattest(X=skjt[1:500,],na=3,nb=3,alpha=0.05)
``` |

Yuande/MBttest documentation built on April 29, 2018, 8:05 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.