## create a PH with specific parameters
(param1 <- 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)))
## create a CF1 with specific parameters
(param2 <- cf1(alpha=c(1,0,0), rate=c(1.0,2.0,3.0)))
## create a hyper Erlang with specific parameters
(param3 <- herlang(shape=c(2,3), mixrate=c(0.3,0.7), rate=c(1.0,10.0)))
## variance
ph.var(param1)
ph.var(param2)
ph.var(param3)
Run the code above in your browser using DataLab