I am trying to fit a generalized extreme value distribution usinf the functioin of "gevfit" as well as the "fitdist" to a rainfall data set which is as follows
1028.5
860.3
1088.2
711.6
874.3
776.2
973.7
457.9
pdGEV=fitdist(sample,'Generalizedextremevalue'); % probability distribution object
paraGEV= gevfit(sample); % paramter estimation
Parameters are calculated are as follows 'shape' k = -1.03801443465120 'scale' sigma = 273.595539473463 'location' mu = 824.624147545310 but there is a warning messageWarning: Maximum likelihood estimation did not converge. Function evaluation limit exceeded. > In gevfit (line 115) In prob.GeneralizedExtremeValueDistribution.fit (line 170) In fitdist>localfit (line 238) In fitdist (line 185) In newcode (line 271) Warning: Maximum likelihood has converged to a boundary point of the parameter space.Confidence intervals and standard errors cannot be computed reliably. > In gevfit (line 127) In prob.GeneralizedExtremeValueDistribution.fit (line 170) In fitdist>localfit (line 238) In fitdist (line 185) In newcode (line 271) Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 7.552473e-21. > In gevlike (line 94) In prob.GeneralizedExtremeValueDistribution.fit (line 171) In fitdist>localfit (line 238) In fitdist (line 185) In newcode (line 271) Warning: Maximum likelihood estimation did not converge. Function evaluation limit exceeded. > In gevfit (line 115) In newcode (line 272) Warning: Maximum likelihood has converged to a boundary point of the parameter space.Confidence intervals and standard errors cannot be computed reliably. > In gevfit (line 127) In newcode (line 272)
How do I solve this problem?
Best Answer