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]]

