Nothing
#####################################################################################
# Update functions for the parameters which for an explicit solution does not exist #
# Solutions are calculated inside EM using root-finding algorithms #
#####################################################################################
# Update function for r.nb
# ------------------------
r.nb.update <- function(r, j, M, eta){
y <- sum(eta[1:M, j])
z <- sum(eta[1:M, j] * (1:M - 1 + r))
x <- sum(eta[1:M, j] * (digamma(1:M + r - 1) - digamma(r) + log(r * y / z)))
return(x)
}
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.