## 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)))
## generate 10 samples
rphase(n=10, ph=phdist)
Run the code above in your browser using DataLab