gettstat: Calculate the t-statistic

In cuiyingbeicheng/BEDFC: Log fold change detection based on an empirical bayesian method.

Description

Calculate the t-statistic to detect if the ratio of $mu^g_2$ and $mu^g_1$ (i.e., fold change) is within or outside a region of interest.

Usage

gettstat(alpha, beta, lfc, null.lfc)

Arguments

alpha	shape parameters
beta	scale parameters
lfc	The log ratio with base 2 of gene expression between treatment and control group
null.lfc	The log ratio with base 2 of gene expression between treatment and control group under the null hypothesis

Details

The t-statistic adopting the proposed bayesian approach is constructed as

$t_g = barY^g - d_0 /sqrttildesigma^2_g/n_2 = sum^n_2_j=1Y^g_j/n_2 - d_0 /sqrttildesigma^2_g/n_2$,

where $d_0$ can be replaced by $d_1$ and $d_2$ to test regions with two sided boundaries.

Value

naive	T-statistic
bayes	T-statistic adopting the proposed bayesian approach

cuiyingbeicheng/BEDFC documentation built on Nov. 8, 2019, 12:34 a.m.