[R] beta-binomial

Sat Jan 11 14:45:32 CET 2003

On 11 Jan 2003 at 20:20, Zhen Pang wrote:

You can finf the beta-binomila distribution in one of Jim Lindsay's 
libraries, in rmutil (dbetabinom      Density of Beta Binomial 
Distribution) and in repeated (biv.betab       Bivariate Beta-
binomial Regression). If you don't want to fit a regression model, 
just a distribution, you could combine with fitdistr in package MASS. 
An example:

> x <- rbetabinom(100, 10, 0.3,3)
> fitdistr(x, dbetabinom, list(m=0.5, s=4), size=10, method="L-BFGS-B",
     lower=c(0,0), upper=c(1,Inf) )
       m            s     
  0.22454182   4.39593104 
 (0.02145044) (1.01358666)

Kjetil Halvorsen

> 
> 
> Does anyone have R functions or library to fit a beta-binomial distribution 
> with glm? Thanks.
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> http://www.stat.math.ethz.ch/mailman/listinfo/r-help