nhppSpikeConstWindow: Spike NHPP rate with constant window width

In seqCBS: Copy Number Profiling using Sequencing and CBS

Description

Randomly spike the smoothed control rate of an NHPP according to the parameters, with constant window width.

Usage

nhppSpikeConstWindow(smoothRates, nSpike = 25, cptLen = 5, minGain = 1.5, 
  maxGain = 10, minLoss = 0.01, maxLoss = 0.5, pGain = 0.6)

Arguments

smoothRates	The smoothed rate estimate of the control process
nSpike	Number of signal spikes
cptLen	Window width of each signal region
minGain	Minimal Gain relative CN
maxGain	Maximal Gain relative CN
minLoss	Minimal Loss relative CN
maxLoss	Maximal Loss relative CN
pGain	Proportion of signal regions that are CN gain

Details

The signal strength are uniformly drawn between the two limits.

Value

spikeMat	A matrix containing the actual signal spike information
caseRates	The rate of the case PP to be simulated after signal spike

Author(s)

Jeremy J. Shen

See Also

nhppSimulate

seqCBS documentation built on May 2, 2019, 9:15 a.m.