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

Description Usage Arguments Value Examples

Description

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

Usage

1
2
3
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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.