glmDM-methods: Method glmDM

Description Usage Arguments Value See Also Examples

Description

glmDM perform inference and estimation on RNA differential modification log2FC.

GLMs with interactive design between dummy variables of IP/input and Treatment/control are fitted for each peaks/sites:

log2(Q) = intercept + I(Treatment) + I(IP) + I(IP)*I(Treatment)

The log2FC and the associated statistics are based on the coefficient estimate of the interactive term: I(IP)*I(Treated).

Under default setting, the returned log2FC are the RR estimates with Couchey priors defined in apeglm.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
glmDM(
  sep,
  glm_type = c("DESeq2", "NB", "Poisson"),
  LFC_shrinkage = c("apeglm", "ashr"),
  ...
)

## S4 method for signature 'SummarizedExomePeak'
glmDM(
  sep,
  glm_type = c("DESeq2", "NB", "Poisson"),
  LFC_shrinkage = c("apeglm", "ashr", "none"),
  ...
)

Arguments

sep

a SummarizedExomePeak object.

glm_type

a character speciefies the type of Generalized Linear Model (GLM) fitted for the purpose of statistical inference during peak calling, which can be one of the c("DESeq2", "NB", "Poisson").

DESeq2

Fit the GLM defined in the function DESeq, which is the NB GLM with regulated estimation of the overdispersion parameters.

NB

Fit the Negative Binomial (NB) GLM.

Poisson

Fit the Poisson GLM.

By default, the DESeq2 GLMs are fitted on the data set with > 1 biological replicates for both the IP and input samples, the Poisson GLM will be fitted otherwise.

LFC_shrinkage

a character for the method of emperical bayes shrinkage on log2FC, could be one of c("apeglm", "ashr", "none"); Default = "apeglm".

see lfcShrink for details; if "none" is selected, only the MLE will be returned.

...

Optional arguments passed to DESeq

Value

a SummarizedExomPeak object.

See Also

glmM

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
### Load the example SummarizedExomPeak object
f1 = system.file("extdata", "sep_ex_dm.rds", package="exomePeak2")

sep <- readRDS(f1)

### Normalize the GC contents biases
sep <- normalizeGC(sep)

### Calculate GLM Statistics on the Modification Peaks
sep <- glmDM(sep)

exomePeak2 documentation built on Nov. 8, 2020, 5:27 p.m.