# \donttest{
## Generate random normally distributed observations.
junk<-rnorm(1000,3,2)
## Fit the dataset using fun.data.ml
fit<-cbind(fun.RMFMKL.ml(junk),starship(junk)$lambda)
## Compare the resulting fits. It is usually the case the maximum likelihood
## provides better estimation of the moments than the starship method.
fun.comp.moments.ml.2(fit,junk)
# }
Run the code above in your browser using DataLab