N<-10 #Number of extractions
beta<-c(0.5, -0.8, 1.0, 0, 0.4, -0.7) #Coefficient vector
p<-length(beta)
r<-2 #Number of success parameter
aux_cov<-rnorm((p-1)*N, 0,1)
Covariates<-data.frame(matrix(aux_cov, ncol=p-1, nrow=N)) #Generating the Covariates data.frame
colnames(Covariates)<-c("X1", "X2", "X3", "X4", "X5")
base<-gen_base_NegBinomial_reg(N, beta, r, Covariates) #Generating the data
base
Run the code above in your browser using DataLab