Hi all, I am attempting to fit, and test the goodness of fit of, a mixed exponential distribution to my dataset which consists of 15minute rainfall intensity data. FYI, the dataset spanning approx.2 years and 7 rainfall stations consists of some three hundred thousand 15min data records, of which some 30 thousand are non-zero rainfall amounts.
Could anyone please tell me how i could do this/which package i could use in R? I am familiar with MASS functions fitdistr which i've used for the single exponential distributions, but have not yet successfully found out how to fit mixed distributions. The single gamma and exponential distributions did not fit the data well at all, hence my interest in fitting a mixed exponential distribution. Also, i thought of formulating the likelihood function and using mle() in package stats4-- but had trouble writing down the likelihood function for the case of a mixed exponential distribution... thanks, vishal Vishal Mehta PhD Candidate Crop and Soil Science Cornell University ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
