DRIMSeq (version 1.0.2)

plotTest: Plot p-values distribution

Description

Plot p-values distribution

Usage

plotTest(x, ...)
"plotTest"(x, out_dir = NULL)
"plotTest"(x, out_dir = NULL)

Arguments

x

dmDStest or dmSQTLtest object.

...

Other parameters that can be defined by methods using this generic.

out_dir

Character string that is used to save the plot in paste0(out_dir, plot_name, ".pdf") file. plot_name depends on type of a plot produced, for example, plot_name = "hist_features" for histogram with number of features per gene. If NULL, the plot is returned as ggplot object and can be further modified, for example, using theme().

Value

Plot a histogram of p-values.

See Also

data_dmDSdata, data_dmSQTLdata, plotData, plotDispersion, plotFit

Examples

###################################
### Differential splicing analysis
###################################
# If possible, use BPPARAM = BiocParallel::MulticoreParam() with more workers

d <- data_dmDSdata

### Filtering
# Check what is the minimal number of replicates per condition 
table(samples(d)$group)
d <- dmFilter(d, min_samps_gene_expr = 7, min_samps_feature_expr = 3, 
 min_samps_feature_prop = 0)

### Calculate dispersion
d <- dmDispersion(d, BPPARAM = BiocParallel::SerialParam())

### Fit full model proportions
d <- dmFit(d, BPPARAM = BiocParallel::SerialParam())

### Fit null model proportions and test for DS
d <- dmTest(d, BPPARAM = BiocParallel::SerialParam())

plotTest(d)