Nothing
nhppSpike <-
function(smoothRates, nSpike=25, cptLenR=4, cptLenMean=10, minGain=1.5, maxGain=10, minLoss=0.01, maxLoss=0.5, pGain=0.6) {
grid.mid = smoothRates$x
spikeRate = smoothRates$y
gridSize = grid.mid[2] - grid.mid[1]
grid.fix = grid.mid - gridSize/2
nGrid = length(grid.fix)
gridL = sample(1:nGrid, nSpike)
cptGridLen = rnbinom(nSpike, size=cptLenR, mu=cptLenMean)
gridR = gridL+cptGridLen
gridR[gridR>nGrid] = nGrid
relCN = sample(0:1, nSpike, replace=TRUE, prob=c(1-pGain, pGain))
relCN[1] = 0
relCN[2] = 1
relCN[relCN==0] = runif(sum(relCN==0), min=minLoss, max=maxLoss)
relCN[relCN==1] = runif(sum(relCN==1), min=minGain, max=maxGain)
spikeMat = cbind(gridL, gridR, grid.fix[gridL], grid.fix[gridR]+gridSize, relCN)
colnames(spikeMat) = c("gridL", "gridR", "readL", "readR", "relCN")
for(i in 1:nSpike) {
spikeRate[gridL[i]:gridR[i]] = spikeRate[gridL[i]:gridR[i]]*relCN[i]
}
spikeRate = spikeRate/mean(spikeRate)*mean(smoothRates$y)
caseRates = list(x=grid.mid, y=spikeRate)
return(list(spikeMat = spikeMat, caseRates=caseRates))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.