LikelihoodRatioTest {grandR}R Documentation

Compute a likelihood ratio test.

Description

The test is computed on any of total/old/new counts using DESeq2 based on two nested models specified using formulas.

Usage

LikelihoodRatioTest(
  data,
  name = "LRT",
  mode = "total",
  normalization = mode,
  target = ~Condition,
  background = ~1,
  no4sU = FALSE,
  columns = NULL,
  verbose = FALSE
)

Arguments

data

A grandR object

name

the user defined analysis name to store the results

mode

either "total", "new" or "old"

normalization

normalize on "total", "new", or "old" (see details)

target

formula specifying the target model (you can use any column name from the Coldata(data))

background

formula specifying the background model (you can use any column name from the Coldata(data))

no4sU

Use no4sU columns (TRUE) or not (FALSE)

columns

logical vector of which columns (samples or cells) to use (or NULL: use all)

verbose

Print status updates

Details

This is a convenience wrapper around the likelihood ratio test implemented in DESeq2.

DESeq2 by default performs size factor normalization. When computing differential expression of new RNA, it might be sensible to normalize w.r.t. to total RNA, i.e. use the size factors computed from total RNA instead of computed from new RNA. This can be accomplished by setting mode to "new", and normalization to "total"!

Value

a new grandR object including a new analysis table. The columns of the new analysis table are


[Package grandR version 0.2.0 Index]