dmTest: Likelihood ratio test

Description Usage Arguments Value Author(s) See Also Examples

Description

First, estimate the null Dirichlet-multinomial model proportions, i.e., feature ratios are estimated based on pooled (no grouping into conditions) counts. Use the likelihood ratio statistic to test for the difference between feature proportions in different groups to identify the differentially spliced genes (differential splicing analysis) or the sQTLs (sQTL analysis).

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
dmTest(x, ...)

## S4 method for signature 'dmDSfit'
dmTest(x, compared_groups = levels(samples(x)$group),
  prop_mode = "constrOptimG", prop_tol = 1e-12, verbose = 0,
  BPPARAM = BiocParallel::MulticoreParam(workers = 1))

## S4 method for signature 'dmSQTLfit'
dmTest(x, prop_mode = "constrOptimG",
  prop_tol = 1e-12, verbose = 0,
  BPPARAM = BiocParallel::MulticoreParam(workers = 1))

Arguments

x

dmDSfit or dmSQTLfit object.

...

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

compared_groups

Vector that defines which experimental conditions should be tested for differential splicing. By default, we test for a difference between any of the groups specified in samples(x)$group. Values in this vector should indicate levels or numbers of levels in samples(x)$group.

prop_mode

Optimization method used to estimate proportions. Possible values "constrOptim" and "constrOptimG".

prop_tol

The desired accuracy when estimating proportions.

verbose

Numeric. Definie the level of progress messages displayed. 0 - no messages, 1 - main messages, 2 - message for every gene fitting.

BPPARAM

Parallelization method used by bplapply.

Value

Returns a dmDStest or dmSQTLtest object.

Author(s)

Malgorzata Nowicka

See Also

data_dmDSdata, data_dmSQTLdata, plotTest, dmDispersion, dmFit

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
###################################
### 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)

head(proportions(d))
head(statistics(d))
head(results(d))


Search within the DRIMSeq package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.