Nothing
pgp <- function(y, theta, lambda) {
n <- length(y)
prob <- numeric(n)
x <- 0:max(y)
a <- log( theta ) + (x - 1) * log( theta + lambda * (x) ) - theta - x * lambda - lgamma( x + 1 )
b <- cumsum( exp(a) )
for ( i in 1:n ) prob[i] <- sum( b * ( y[i] == x ) )
prob
}
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.