testLrt: Calculate likelihood ratio test for the specific nested model

In MPRAnalyze: Statistical Analysis of MPRA data

Description

Calculate likelihood ratio test for the specific nested model

Usage

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

data <- simulateMPRA(tr = rep(2,5), da=c(rep(2,2), rep(2.5,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)

MPRAnalyze documentation built on Nov. 8, 2020, 8:22 p.m.