testLrt: Calculate likelihood ratio test for the specific nested model

Description Usage Arguments Value Note Examples

View source: R/hypothesisTesting.R

Description

Calculate likelihood ratio test for the specific nested model

Usage

1
testLrt(obj)

Arguments

obj

the MpraObject containing the full and reduced

Value

results data frame

Note

Must be run after running an LRT-based analysis

Examples

1
2
3
4
5
6
7
8
9
data <- simulateMPRA(tr = rep(2,5), da=c(rep(0,2), rep(1,3)), 
                     nbatch=2, nbc=15)
obj <- MpraObject(dnaCounts = data$obs.dna, 
                  rnaCounts = data$obs.rna, 
                  colAnnot = data$annot)
obj <- estimateDepthFactors(obj, lib.factor = "batch", which.lib = "both")
obj <- analyzeComparative(obj, dnaDesign = ~ batch + barcode + condition, 
                              rnaDesign = ~ condition, reducedDesign = ~ 1)
results <- testLrt(obj)

YosefLab/MPRAnalyze documentation built on Nov. 14, 2020, 2:35 a.m.