## make sample
wsample <- rweibull(n=100, shape=2, scale=1)
## PH fitting for general PH
(result1 <- phfit.point(ph=ph(2), x=wsample))
## PH fitting for CF1
(result2 <- phfit.point(ph=cf1(2), x=wsample))
## PH fitting for hyper Erlang
(result3 <- phfit.point(ph=herlang(3), x=wsample))
## mean
ph.mean(result1$model)
ph.mean(result2$model)
ph.mean(result3$model)
## variance
ph.var(result1$model)
ph.var(result2$model)
ph.var(result3$model)
## up to 5 moments
ph.moment(5, result1$model)
ph.moment(5, result2$model)
ph.moment(5, result3$model)
Run the code above in your browser using DataLab