Nothing

```
expsnps <-
function(x, m.rate, c.rate, tau) {
veclen <- length(x)
sumterm <- numeric(veclen)
for (k in 1:veclen) {
for (i in 0:x[k]) {
sumterm[k] <- sumterm[k] + dgeom(i, c.rate/(2*m.rate+c.rate))*dpois(x[k]-i, m.rate*tau)
}
}
return(sumterm)
}
```

**Any scripts or data that you put into this service are public.**

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.