rnaSeqMap (version 2.30.0)

generators: Generators for synt data and

Description

Various generators for experiments.

Usage

generatorAddSquare(nd, deg, length.prop=0.5)
generatorAdd(nd, deg, length.prop=0.5)
generatorMultiply(nd, deg, length.prop=0.5)
generatorTrunc(nd,deg)
generatorSynth(nd, deg, length.prop=0.5)
generatorPeak(nd, deg, sr=10, mult=10)

Arguments

nd

nucleotide distribution object

deg

degeneration level for the output profile

length.prop

a fraction of the genome region to be degenerated - (0,1)

sr

distance from the 5' end for the peak

mult

multiplier - how many times the peak is supposed to be higher than the maximum of the distribution

Examples

   if (1==0)
   {
   rs <- newSeqReads('chr2', 220238268, 220254744, -1)
   f <- c("test1.bam", "test2.bam", "test3.bam", "test4.bam", "test5.bam")
   ff <- sapply(f, function(x) system.file("extdata", x, package = "rnaSeqMap"))
   rs <- getBamData(rs, 1:5)    	 
   nd <- getCoverageFromRS(rs, 1:5) 
   generatorTrunc(nd,0.5)
   }