# NOT RUN {
params <- c(4.114,0.1333,0.0193,0.1588) #taken straight from some paper
gld.rv <- rgld(100,params)
hist(gld.rv)
fit.davies.q(gld.rv) #remember the Davies distn has 3 DF and the GLD 4...
# }
Run the code above in your browser using DataLab