## create a PH with specific parameters
(phdist <- ph(alpha=c(1,0,0),
Q=rbind(c(-4,2,0),c(2,-5,1),c(1,0,-1)),
xi=c(2,2,0)))
## c.d.f. for 0, 0.1, ..., 1
pphase(q=seq(0, 1, 0.1), ph=phdist)
Run the code above in your browser using DataLab