plot_genewise_moments: Visualize gene-wise moments

Description Usage Arguments Value Examples

Description

Visualize gene-wise moments

Usage

1
plot_genewise_moments(data, batch, robust)

Arguments

data

Given data or simulated data from rnaseq_sim()

batch

Batch covariate

robust

Boolean indicator of using robust (TRUE) or non-robust test (FALSE) in visualization

Value

Gene-wise moments

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
nbatch <- 3
ncond <- 2
npercond <- 10
data.matrix <- rnaseq_sim(ngenes=50, nbatch=nbatch, ncond=ncond, npercond=
    npercond, basemean=10000, ggstep=50, bbstep=2000, ccstep=800, 
    basedisp=100, bdispstep=-10, swvar=1000, seed=1234)
batch <- rep(1:nbatch, each=ncond*npercond)
condition <- rep(rep(1:ncond, each=npercond), nbatch)
data_adjusted <- batchQC_condition_adjusted(data.matrix, batch, condition)
gene_moments <- plot_genewise_moments(data_adjusted, batch, robust=FALSE)

BatchQC documentation built on Nov. 8, 2020, 8:30 p.m.