question on nbinom1

Thank you.
This is perfectly clear and enables further analysis.
I hope you can arrange to put it in some (even better, much) documentation.
I wish it were in wikipedia.

Ben Bolker writes:
 > 
 >    I believe the R d/p/q/r functions corresponding to glmmTMB's 
 > implementation of *nbinom1 would look like this:
 > 
 > 
 > rnbinom1 <- function(n, mu, phi) {
 >      ## var = mu*(1+phi) = mu*(1+mu/k) -> k = mu/phi
 >      rnbinom(n, mu=mu, size=mu/phi)
 > }
 > 
 > dnbinom1 <- function(x, mu, phi, ...) {
 >      dnbinom(n, mu=mu, size=mu/phi, ...)
 > }
 > 
 > pnbinom1 <- function(q, mu, phi, ...) {
 >      pnbinom(q, mu=mu, size=mu/phi, ...)
 > }
 > 
 > qnbinom1 <- function(p, mu, phi, log=FALSE) {
 >      pnbinom(p, mu=mu, size=mu/phi, ...)
 > }
 > 
 > 
 >    (there would be an even more clever/inscrutable way to do this by 
 > transforming the body of the code, without repeating oneself so much, 
 > but it would probably be a bad idea)
 >

question on nbinom1

Thread (6 messages)