getDistrParam: Get model distribution parameters from an MpraObject of a...
In YosefLab/MPRAnalyze: Statistical Analysis of MPRA data

Description Usage Arguments Value Examples

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

1
2
3

getDistrParam_DNA(obj, enhancer, full = TRUE)

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

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

fit parameters (numeric, samples x parameters)

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