extractAIC.rglmb: Extract DIC from a Fitted Bayesian Model

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

Description

Computes the Deviance Information Criterion for a fitted Bayesian Generalized Linear Model

Usage

## S3 method for class 'rglmb'
extractAIC(fit, ...)

Arguments

fit	an object of class "rglmb", typically the result of a call to rglmb
...	further arguments to or from other methods

Value

A list with the Estimated effective number of parameters pD and the DIC from the object fit of class "glmb".

Examples

## ----dobson-------------------------------------------------------------------
## 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)

## Prior mean vector 
mu<-matrix(0,5)           
mu[1,1]=log(mean(counts)) 
## Prior standard deviation and Variance
mysd<-1           
V=((mysd)^2)*diag(5)  
## Call to glmb
glmb.D93<-glmb(n=1000,counts ~ outcome + treatment,
               family = poisson(),pfamily=dNormal(mu=mu,Sigma=V))
## ----glmb extractAIC-------------------------------------------------------------
extractAIC(glmb.D93)

knygren/glmbayes documentation built on Sept. 4, 2020, 4:39 p.m.