disc_phase_type <- matrix(c(0.4, 0, 0.2,
0.5, 0.3, 0.2,
0, 0.7, 0.2), ncol = 3)
Y <- DPH(disc_phase_type)
dDPH(3:4, Y)
pDPH(5, Y)
qDPH(0.5, Y)
set.seed(0)
rDPH(6, Y)
rFullDPH(Y)
Run the code above in your browser using DataLab