simulateMPRA: Simulate an MPRA dataset
In MPRAnalyze: Statistical Analysis of MPRA data

Description Usage Arguments Details Value Examples

Simulate an MPRA dataset

simulateMPRA(tr = rep(2, 100), da = c(rep(0, ceiling(length(tr)/2)),
  rep(0.5, floor(length(tr)/2))), dna.noise.sd = 0.2,
  rna.noise.sd = 0.3, dna.inter = 5, dna.inter.sd = 0.5, nbc = 100,
  coef.bc.sd = 0.5, nbatch = 3, coef.batch.sd = 0.5)

`tr`	a vector of the true transcription rates. The length of the vector determines the number of enhancers included in the dataset Default is 100 enhancers of identical transcription rate of 2.
`da`	a vector of the differential activity signal. Must be the same length as tr. if NULL (default) only one condition is simulated, with no differential activity. The values provided are used as the logFold Change between the conditions, treating the tr vector as the reference condition. For non-differentially active enhancers, this value should be 0.
`dna.noise.sd`	level of noise to add to the DNA library
`rna.noise.sd`	level of noise to add to the RNA library
`dna.inter`	the baseline DNA levels (intercept term), controlling the true mean abundance of plasmids
`dna.inter.sd`	the true variation of the plasmid levels
`nbc`	number of unique barcode to include per enhancer
`coef.bc.sd`	true variation between barcodes
`nbatch`	number of batches to simulate
`coef.batch.sd`	the level of true variation that distinguishes batches (the size of the batch effects)

the data is generated by using the same nested-GLM construct that MPRAnalyzes uses, with non-strandard log-normal noise models (whereas by default MPRAnalyze uses a Gamma-Poisson model). The data generated can have multiple batches, and either 1 or 2 conditions, and the simulated data is always paired (DNA and RNA extracted from the same library). User can control both true and observed variation levels (noise), the number of expected plasmids per barcode, the true transcription ratio, the size of the batch and barcode effects.

a list:

true.dna The true dna abundances
obs.dna the observed dna counts
true.rna the true rna abundances
obs.rna the observed rna counts
annot the annotations data.frame for each sample

data <- simulateMPRA()
# single condition
data <- simulateMPRA(da=NULL)
# more observed noise
data <- simulateMPRA(dna.noise.sd = 0.75, rna.noise.sd = 0.75)
# gradually increasing dataset
data <- simulateMPRA(tr = seq(2,3,0.01), da=NULL)