getDistrParam_DNA {MPRAnalyze} | R Documentation |
Get model distribution parameters from an MpraObject of a given candidate enhancer
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), 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)