inst/examples/Ex_dummy.coef.glmb.R

In knygren/glmbayes: Bayesian Generalized Linear Models (iid Samples)

set.seed(333)
## Dobson (1990) Page 93: Randomized Controlled Trial :
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
mysd<-1
mu<-matrix(0,5)
mu[1]=log(mean(counts))
V0<-((mysd)^2)*diag(5)
glmb.D93<-glmb(n=1000,counts ~ outcome + treatment, 
family = poisson(),pfamily=dNormal(mu=mu,Sigma=V0))
summary(glmb.D93)


myd=dummy.coef(glmb.D93)
print(myd)