plot_gev.m

极值理论中各种函数及图像的程序。matlab实现。

function c=plot_gev(gevres),
%Menu driven plotting facility for GEV distribution fit
%
% USAGE: plot_gev(gevres)
%
%gevres: Structure obtained from a former GEV fit i.e. gevres=gev(...)
%
residuals=(1 + (gevres.par_ests(1) * (gevres.data - gevres.par_ests(3)))/gevres.par_ests(2)).^(-1/gevres.par_ests(1));

r=1;
while r~=0, 
clc
r=input('\n\n\n\n\n\n\n\n Enter your choice for plotting !!! \n\n Enter (1) for scatterplot of residuals \n Enter (2) for QQ plot of residuals \n Enter (0) to exit \n                 >');
while (r~=1) & (r~=2) & (r~=0),
    r=input('!!!Enter either 0,1,2 >');
end
if r==1,
    plot(residuals,'.');
    sm=csaps(1:length(residuals),residuals(1:end),0.00005,1:length(residuals));
    hold on
    plot(1:length(residuals),sm,'k')
    hold off
    out=[residuals,sm'];    
    xlabel('Ordering');
    ylabel('Residuals');
    
 end

if r==2,
    out=qplot(residuals,0);
end

end