DPpsBF: Computes Pseudo Bayes Factors from DPpackage output

Description Usage Arguments Author(s) Examples

Description

This function computes Pseudo Bayes Factors from DPpackage output.

Usage

1

Arguments

...

DPpackage output objects. These have to be of the same class.

Author(s)

Alejandro Jara <atjara@uc.cl>

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
## Not run: 
    # Respiratory Data Example

      data(indon)
      attach(indon)

      baseage2 <- baseage**2
      follow <- age-baseage
      follow2 <- follow**2 

    # Prior information

      beta0 <- rep(0,9)
      Sbeta0 <- diag(1000,9)
      tinv <- diag(1,1)
      prior <- list(a0=2,b0=0.1,nu0=4,tinv=tinv,
                    mub=rep(0,1),Sb=diag(1000,1),
                    beta0=beta0,Sbeta0=Sbeta0)

    # Initial state
      state <- NULL

    # MCMC parameters

      nburn <- 5
      nsave <- 100
      nskip <- 5
      ndisplay <- 100
      mcmc <- list(nburn=nburn,nsave=nsave,nskip=nskip,ndisplay=ndisplay)

    # Fit the Probit model
      fit1 <- DPglmm(fixed=infect~gender+height+cosv+sinv+xero+baseage+
                     baseage2+follow+follow2,random=~1|id,
                     family=binomial(probit),
                     prior=prior,mcmc=mcmc,state=state,status=TRUE)

    # Fit the Logit model
      fit2 <- DPglmm(fixed=infect~gender+height+cosv+sinv+xero+baseage+
                     baseage2+follow+follow2,random=~1|id,
                     family=binomial(logit),
                     prior=prior,mcmc=mcmc,state=state,status=TRUE)

    # Model comparison
      DPpsBF(fit1,fit2)



## End(Not run)

DPpackage documentation built on May 1, 2019, 10:23 p.m.

Related to DPpsBF in DPpackage...