egumbel.m

来自「一个非常实用的统计工具箱」· M 代码 · 共 13 行

function [a, u] = egumbel(x)
%EGUMBEL  Estimate parameters of the Gumbel distribution
%
%         [a, u] = egumbel(x)
%
%	  The parameter  a  is scaling and  u  is position. The
%	  method of moments is used.

%       GPL Copyright (c) Anders Holtsberg, 1999

a = mean(x) - .57721566490153286; % Euler's gamma
u = std(x) / 1.2825498301618641 % pi/sqrt(6)