getDistrParam: Get model distribution parameters from an MpraObject of a...

In MPRAnalyze: Statistical Analysis of MPRA data

Description

Get model distribution parameters from an MpraObject of a given candidate enhancer

Usage

getDistrParam_DNA(obj, enhancer, full = TRUE)

getDistrParam_RNA(obj, enhancer = NULL, full = TRUE)

Arguments

obj	MpraObject to extract from
enhancer	enhancer to extract
full	whether to extract from full model

Value

fit parameters (numeric, samples x parameters)

Examples

data <- simulateMPRA(tr = rep(2,5), da=NULL, 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 <- analyzeQuantification(obj, dnaDesign = ~ batch + barcode, 
                              rnaDesign = ~1)
## get distributional parameters of the first enhancer:
dist.params.dna <- getDistrParam_DNA(obj, 1)
dist.params.rna <- getDistrParam_RNA(obj, 1)

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