Beta-Binomial Regression in R

Hi All:

I am trying to generate Beta-Binomial data with regressors using R.  
I have
used the following code to generate Beta-Binomial data.  Now I want  
to add
a covariate to the equation. I would then like to use the simulated  
data to
run the Beta-Binomial model with covariates on it. Appreciate any  
help.

set.seed(111)
k<-20
n<-60
x<-NULL

p<-rbeta(k,3,3)# so that the mean nausea rate is alpha/(alpha+beta)
for(i in 1:k)
x<-cbind(x,rbinom(300,n,p[i]))

Thanks
Anamika

	[[alternative HTML version deleted]]