plot_genewise_moments {BatchQC} | R Documentation |
Visualize gene-wise moments
plot_genewise_moments(data, batch, robust)
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 |
Gene-wise moments
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)