# \donttest{
# Generate random normally distributed observations.
junk<-rnorm(1000,3,2)
# Fit the dataset using fun.data.ml
fit<-fun.data.fit.ml(junk)
# 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(fit,junk)
# }
Run the code above in your browser using DataLab