There are specific functions for obtaining AIC, confidence intervals, and distribution parameter estimates related to each fitting functions respectively. Details are provided below:
1) The fitting functions ("evfit", "expfit", "gamfit", "lognfit", "normfit", "wblfit") contain an option to specify confidence intervals. Estimates under certain confidence intervals will be returned.
2) The likelihood functions ("evlike", "explike", "gamlike", "lognlike", "normlike", "wbllike") returns the negative log-likelihood given estimates from the previous step as input.
3) The "aicbic" function returns Akaike information criteria (AIC) given log-likelihood and number of estimated parameters.
Taking the Extreme value parameter estimates as an example...
1) "parmhat = evfit(data)" returns the estimates for 95% confidence intervals by default. It can also return estimates for defined confidence interval from the "alpha" value via the command "[parmhat,parmci] = evfit(data,alpha)".
2) "nlogL = evlike(params,data)" returns the negative of the log-likelihood for the extreme value distribution.
3) "aicbic(logL, numParam)" returns AIC given the log-likelihood and number of estimated parameters from previous two steps.
Best Answer