two_conditions_analysis_Gene.Rd
This function is used to perform differential ASE detection for two conditions analysis of a given gene.
two_conditions_analysis_Gene( dat, phased = FALSE, varList = NULL, adaptive = TRUE, n_resample = 10^6, parallel = FALSE, n_core = 1 )
dat: | bulk RNA-seq data of a given gene. Must contain variables
|
---|---|
phased: | a logical value indicates whether the haplotype phase of the data is known or not. Default is FALSE |
varList: | a character string specifies fomula of covariates that users want to adjusted in the model. An example could be "`var1`+`var2`". Default is NULL |
n_resample: | a numeric value indicates the maximum number of resamplings performed to obtain estimated p-value. Default is 10^6 |
adaptive: | a logical value indicates whether the resampling is done through an adaptive procedure or not. Only applies when n_resample >= 10^3. Default is TRUE |
parallel: | a logical value indicates whether do parallel computing for the resampling precedure or not. Default is FALSE |
n_core: | a numeric value indicates number of clusters used for parallel computing when parameter "parallel" is set to TRUE. Default is 1 |
A vector with two elements:
`LRT statistic`: numeric, the likelihood ratio test (LRT) statistics of differential ASE effect;
`p-value`: the estimated p-value of the LRT statistic;