DP-class: Dirichlet Process DP applied to censored data

Description Examples

Description

Dirichlet Process DP applied to censored data

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
## Not run: 
weights <- matrix(c(1,0,0,0,1,0,0,0,1), ncol=3)

data <- sim.data(n=100, J=10, weights)

G1 <- init.DP(DataStorage=data, prior=list(mu=0, n=0.1, v=3, vs2=1*3), L=35, thinning=50,
          burnin = 5000, max_iter = 55000 )
          
G1 <- MCMC.DP(G1, data, 55000)

validate.DP(G1, data)
 
plot.ICDF(G1@theta, G1@phi, G1@weights, G1@L, grid=0:500,
           distribution=data@presentation, xlim=500)

## End(Not run)

AlexPiche/DPsurv documentation built on May 5, 2019, 4:52 a.m.