Skip to content

library/function that estimates parameters of well known distributions from empirical data?

3 messages · Ben Bolker, Ted

Ted
#
I found this a few months ago, but for the life of me I can't remember what
the function or package was, and I have had no luck finding it this week.

I have found, again, the functions for working with distributions like
Cauchy, F, normal, &c., and ks.test, but I have not found the functions for
estimating the distribution parameters given a vector of values.

What I need to do is estimate the distribution parameters for each candidate
distribution, and then test to see which gives the best fit to the data.

I want to examine the question, given this dataset (which may have thousands
of records), does the normal or cauchy distribution fit the data best, and
which what parameters.  It will not be known a priori whether or not the
most appropriate distribution is non-central, though we do know that often
(not always) values of medium size in absolute value are more often positive
than negative and that very large values are more often negative than
positive.

Could someone please give me a gentle reminder of the package and
function(s) I ought to be examining?

Thanks

Ted
#
Ted Byers <r.ted.byers <at> gmail.com> writes:
Look at the fitdistr function in the MASS package.  Consider
AIC comparisons for ranking the fits to these non-nested models.

  good luck
   Ben Bolker
Ted
#
Thanks Ben

That was the one I'd remembered but couldn't find.

Mark Leeds also told me about DistributionFits(fBasics), which I hadn't
seen.  There seems to be only a little overlap between the two.

Could I trouble you to expand on AIC (esp. what the function name and
package is to apply it to the output from these two functions)?  I just read
the help provided for each and neither mentions AIC.

Thanks again Ben

Ted
Ben Bolker wrote: