Kurtoza

Przykład kurtozy - inf. przy legendzie podano w języku japońskim

Kurtoza (z gr. κυρτός, kyrtos, kurtos – wydęty) – jedna z miar kształtu rozkładu wartości cechy statystycznej. W praktyce definiuje się ją najczęściej w formie tzw. kurtozy nadwyżkowej (współczynnika ekscesu, ang. excess kurtosis):

${\mbox{Kurt}}_{E}={\frac {\mu _{4}}{\sigma ^{4}}}-3,$

gdzie:

$\mu _{4}$ – czwarty moment centralny,

$\sigma$ – odchylenie standardowe.

Interpretacja

Klasyczny współczynnik kurtozy definiowany jest jako standaryzowany moment centralny czwartego rzędu:

$\operatorname {Kurt} [X]=\operatorname {E} \left[\left({\frac {X-\mu }{\sigma }}\right)^{4}\right]={\frac {\operatorname {E} \left[(X-\mu )^{4}\right]}{\left(\operatorname {E} \left[(X-\mu )^{2}\right]\right)^{2}}}={\frac {\mu _{4}}{\sigma ^{4}}}.$

Kurtoza nadwyżkowa (inaczej współczynnik ekscesu lub po prostu eksces: ${\mbox{Kurt}}_{E}[X]={\mbox{Kurt}}[X]-3$) jest jednak wygodniejsza, gdyż:

• kurtoza nadwyżkowa rozkładu normalnego wynosi 0,
• jeśli $Y$ jest sumą $n$ niezależnych zmiennych losowych, każdej o rozkładzie identycznym z rozkładem zmiennej losowej $X,$ zachodzi własność: ${\mbox{Kurt}}_{E}[Y]={\mbox{Kurt}}_{E}[X]/n.$

Wbrew stwierdzeniom obecnym w niektórych podręcznikach, kurtoza nie mierzy „spłaszczenia”, „wysmukłości” ani „spiczastości” rozkładu. Na kurtozę ma wpływ intensywność występowania wartości skrajnych, mierzy więc ona, co się dzieje w „ogonach” rozkładu, natomiast kształt „czubka” rozkładu jest praktycznie bez znaczenia.

Rozkłady prawdopodobieństwa można podzielić ze względu na wartość kurtozy na rozkłady:

• mezokurtyczne (Kurt_E = 0) – wartość kurtozy nadwyżkowej wynosi 0, intensywność wartości skrajnych jest podobna do intensywności wartości skrajnych rozkładu normalnego (dla którego kurtoza nadwyżkowa wynosi dokładnie 0),
• leptokurtyczne (Kurt_E > 0) – kurtoza nadwyżkowa jest dodatnia, intensywność wartości skrajnych jest większa niż dla rozkładu normalnego („ogony” rozkładu są „grubsze”),
• platykurtyczne (Kurt_E < 0) – kurtoza nadwyżkowa jest ujemna, intensywność wartości skrajnych jest mniejsza niż w przypadku rozkładu normalnego („ogony” rozkładu są „węższe”).

Kurtoza nadwyżkowa z próby wyraża się wzorem:

${\mbox{Kurt}}_{E}={\frac {{\frac {1}{n}}{\sum _{i=1}^{n}(x_{i}-\mu )^{4}}}{\sigma ^{4}}}-3,$

gdzie:

$x_{i}$ – $i$-ta wartość cechy,

$\mu$ – wartość oczekiwana w populacji,

$\sigma$ – odchylenie standardowe w populacji,

$n$ – liczebność próby.

Powyższa statystyka jest obciążonym estymatorem kurtozy nadwyżkowej z populacji, estymator nieobciążony wyraża się wzorem:

${\mbox{Kurt}}_{E}={\frac {n(n+1)}{(n-1)(n-2)(n-3)}}\sum _{i=1}^{N}\left({\frac {x_{i}-{\bar {x}}}{s}}\right)^{4}-{\frac {3(n-1)^{2}}{(n-2)(n-3)}},$

gdzie:

${\bar {x}}$ – średnia z próby,

$s$ – odchylenie standardowe z próby,

$x_{i}$ – kolejne wartości cechy,

$n$ – liczebność próby.

Obliczenie kurtozy dla rozkładu normalnego

${\mbox{Kurt}}_{E}={\frac {\mu _{4}(X)}{\sigma ^{4}}}-3=0$

Dowód

Niech:

$X\colon \Omega \to \mathbb {R} ,$

${\mbox{Kurt}}_{E}={\frac {\mu _{4}(X)}{\sigma ^{4}}}-3,$

$\mu _{n}(X)=\operatorname {E} ((X-\operatorname {E} X)^{n})$ – moment centralny n–tego rzędu,

$m_{n}(X)=\operatorname {E} ((X)^{n})$ – moment zwykły n–tego rzędu,

Dystrybuanta rozkładu normalnego oraz gęstość prawdopodobieństwa rozkładu normalnego to odpowiednio:

$P_{x}(A)=\int _{A}f(x)dx,\qquad f(x)={\frac {1}{\sigma {\sqrt {2\pi }}}}e^{\frac {-(x-m)^{2}}{2\sigma ^{2}}},dla\qquad x\in \mathbb {R} .$

Wiadomo, że w rozkładzie normalnym:

$m=\operatorname {E} X=m_{1}(X),$

$\sigma ^{2}=\operatorname {D} ^{2}X=\operatorname {E} ((X-\operatorname {E} X)^{2})=\mu _{2}(X).$

Mamy:

a)

$m_{n}(X)=\int _{\Omega }X^{n}(w)P(dw)=\int _{\mathbb {R} }x^{n}f(x)dx,$

b)

${\begin{aligned}\mu _{4}(X)&=\operatorname {E} ((X-\operatorname {E} X)^{4})\\&=\operatorname {E} (X^{4}-4(\operatorname {E} X)X^{3}+6(\operatorname {E} X)^{2}X^{2}-4(\operatorname {E} X)^{3}X+(\operatorname {E} X)^{4})\\&=\operatorname {E} (X^{4})-4m\operatorname {E} (X^{3})+6\operatorname {E} (X^{2})(\operatorname {E} X)^{2}-4(\operatorname {E} X)^{3}m+m^{4}\\&=\operatorname {E} (X^{4})-4m\operatorname {E} (X^{3})+6m^{2}\operatorname {E} (X^{2})-3m^{4}.\end{aligned}}$

Obliczamy momenty zwykłe:

$\operatorname {E} (X^{2})=\int _{-\infty }^{+\infty }x^{2}{\frac {1}{\sigma {\sqrt {2\pi }}}}e^{\frac {-(x-m)^{2}}{2\sigma ^{2}}}dx={\frac {1}{\sigma {\sqrt {2\pi }}}}\int _{-\infty }^{+\infty }(z{\sqrt {2}}\sigma +m)^{2}e^{-z^{2}}{\sqrt {2}}\sigma dz=^{*}$

$z={\frac {x-m}{{\sqrt {2}}\sigma }},\quad x=z{\sqrt {2}}\sigma +m$

$dz={\frac {1}{{\sqrt {2}}\sigma }},\quad dx={\sqrt {2}}\sigma dz$

$=^{*}{\frac {1}{\sqrt {\pi }}}\int _{-\infty }^{+\infty }(2z^{2}\sigma ^{2}+2{\sqrt {2}}z\sigma m+m^{2})e^{-z^{2}}dz=$

$={\frac {1}{\sqrt {\pi }}}{\Bigg (}2\sigma ^{2}\left(\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz\right)+2{\sqrt {2}}\sigma m\left(\int _{-\infty }^{+\infty }ze^{-z^{2}}dz\right)_{=0}+m^{2}\left(\int _{-\infty }^{+\infty }e^{-z^{2}}dz\right)_{={\sqrt {\pi }}}{\Bigg )}=$

$={\frac {1}{\sqrt {\pi }}}\left(2\sigma ^{2}\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz+m^{2}{\sqrt {\pi }}\right)=^{*}$

$\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz=\int _{-\infty }^{+\infty }z\cdot ze^{-z^{2}}dz=-{\frac {1}{2}}ze^{-z^{2}}{\Bigg |}_{-\infty }^{+\infty }+{\frac {1}{2}}\int _{-\infty }^{+\infty }e^{-z^{2}}dz=$

$=-{\frac {1}{2}}ze^{-z^{2}}{\Bigg |}_{-\infty }^{+\infty }+{\frac {1}{2}}\int _{-\infty }^{+\infty }e^{-z^{2}}dz={\frac {\sqrt {\pi }}{2}}$

$=^{*}{\frac {1}{\sqrt {\pi }}}(2\sigma ^{2}{\frac {\sqrt {\pi }}{2}}+m^{2}{\sqrt {\pi }})=\sigma ^{2}+m^{2}$

$\operatorname {E} (X^{3})=\int _{-\infty }^{+\infty }x^{3}{\frac {1}{\sigma {\sqrt {2\pi }}}}e^{\frac {-(x-m)^{2}}{2\sigma ^{2}}}dx={\frac {1}{\sqrt {\pi }}}\int _{-\infty }^{+\infty }(z{\sqrt {2}}\sigma +m)^{3}e^{-z^{2}}dz=^{*}$

$z={\frac {x-m}{{\sqrt {2}}\sigma }},\quad x=z{\sqrt {2}}\sigma +m$

$dz={\frac {1}{{\sqrt {2}}\sigma }},\quad dx={\sqrt {2}}\sigma dz$

$=^{*}{\frac {1}{\sqrt {\pi }}}\int _{-\infty }^{+\infty }(z^{3}{\sqrt {2^{3}}}\sigma ^{3}+3z^{2}2\sigma ^{2}m+3z{\sqrt {2}}\sigma m^{2}+m^{3})e^{-z^{2}}dz=$

$={\frac {2{\sqrt {2}}\sigma ^{3}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{3}e^{-z^{2}}dz{\bigg )}_{=0}+{\frac {6\sigma ^{2}m}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz{\bigg )}_{={\frac {\sqrt {\pi }}{2}}}+$

$+{\frac {3{\sqrt {2}}\sigma m^{2}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }ze^{-z^{2}}dz{\bigg )}_{=0}+{\frac {m^{3}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }e^{-z^{2}}dz{\bigg )}_{={\sqrt {\pi }}}=$

$={\frac {6\sigma ^{2}m}{\sqrt {\pi }}}\cdot {\frac {\sqrt {\pi }}{2}}+m^{3}=3\sigma ^{2}m+m^{3}$

$\operatorname {E} (X^{4})=\int _{-\infty }^{+\infty }x^{4}{\frac {1}{\sigma {\sqrt {2\pi }}}}e^{\frac {-(x-m)^{2}}{2\sigma ^{2}}}dx={\frac {1}{\sqrt {\pi }}}\int _{-\infty }^{+\infty }(z{\sqrt {2}}\sigma +m)^{4}e^{-z^{2}}dz=^{*}$

$z={\frac {x-m}{{\sqrt {2}}\sigma }},\quad x=z{\sqrt {2}}\sigma +m$

$dz={\frac {1}{{\sqrt {2}}\sigma }},\quad dx={\sqrt {2}}\sigma dz$

$=^{*}{\frac {1}{\sqrt {\pi }}}\int _{-\infty }^{+\infty }(4\sigma ^{4}z^{4}+32{\sqrt {2}}\sigma ^{3}mz^{3}+12\sigma ^{2}m^{2}z^{2}+4{\sqrt {2}}\sigma m^{3}z+m^{4})e^{-z^{2}}dz=$

$={\frac {4\sigma ^{4}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{4}e^{-z^{2}}dz{\bigg )}+{\frac {32{\sqrt {2}}\sigma ^{3}m}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{3}e^{-z^{2}}dz{\bigg )}_{=0}+$

$+{\frac {12\sigma ^{2}m^{2}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz{\bigg )}_{={\frac {\sqrt {\pi }}{2}}}+{\frac {4{\sqrt {2}}\sigma m^{3}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }ze^{-z^{2}}dz{\bigg )}_{=0}+{\frac {m^{4}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }e^{-z^{2}}dz{\bigg )}_{={\sqrt {\pi }}}=$

$={\frac {4\sigma ^{4}}{\sqrt {\pi }}}{\bigg (}\int _{-\infty }^{+\infty }z^{4}e^{-z^{2}}dz{\bigg )}+0+{\frac {12\sigma ^{2}m^{2}}{\sqrt {\pi }}}(\sigma ^{2}+m^{2})+0+{\frac {m^{4}}{\sqrt {\pi }}}{\sqrt {\pi }}=^{**}$

$\int _{-\infty }^{+\infty }z^{4}e^{-z^{2}}dz=\int _{-\infty }^{+\infty }z\cdot z^{3}e^{-z^{2}}dz=^{*}$

$u=z,\quad u'=1,\quad v'=z^{3}e^{-z^{2}},\quad v=?$

$v=\int z^{3}e^{-z^{2}}dz={\frac {1}{2}}\int te^{-t}dt={\frac {1}{2}}(-te^{-t})+\int e^{-t}dt={\frac {1}{2}}(-z^{2}e^{-z^{2}}-e^{-z^{2}})$

$t=z^{2}\quad u=t,\quad u'=1,\quad dz={\frac {1}{2z}}dt,\quad v'=e^{-t},\quad v=-e^{-t}$

$=^{*}z\cdot {\frac {1}{2}}(-z^{2}e^{-z^{2}}-e^{-z^{2}}){\Big |}_{-\infty }^{+\infty }-\int _{-\infty }^{+\infty }{\frac {1}{2}}(-z^{2}e^{-z^{2}}-e^{-z^{2}})dz={\frac {1}{2}}\cdot 0+{\frac {1}{2}}\int _{-\infty }^{+\infty }z^{2}e^{-z^{2}}dz+{\frac {1}{2}}\int _{-\infty }^{+\infty }e^{-z^{2}}dz=$

$={\frac {1}{2}}\cdot {\frac {\sqrt {\pi }}{2}}+{\frac {1}{2}}\cdot {\sqrt {\pi }}={\frac {3{\sqrt {\pi }}}{4}}$

$=^{**}{\frac {4\sigma ^{4}}{\sqrt {\pi }}}\cdot {\frac {3{\sqrt {\pi }}}{4}}+{\frac {12\sigma ^{2}m^{2}}{\sqrt {\pi }}}\cdot {\frac {\sqrt {\pi }}{2}}+{\frac {m^{4}}{\sqrt {\pi }}}\cdot {\sqrt {\pi }}=3\sigma ^{4}+6\sigma ^{2}m^{2}+m^{4}\quad (=\operatorname {E} (X^{4}))$

Obliczone wartości:

$\operatorname {E} (X^{2})=\sigma ^{2}+m^{2}$

$\operatorname {E} (X^{3})=3\sigma ^{2}m+m^{3}$

$\operatorname {E} (X^{4})=3\sigma ^{4}+6\sigma ^{2}m^{2}+m^{4}$

podstawiamy do wzoru na czwarty moment centralny z punktu b):

${\begin{aligned}\mu _{4}(X)&=\operatorname {E} (X^{4})-4m\operatorname {E} (X^{3})+6m^{2}\operatorname {E} (X^{2})-3m^{4}\\&=(3\sigma ^{4}+6\sigma ^{2}m^{2}+m^{4})-4m(3\sigma ^{2}m+m^{3})+6m^{2}(\sigma ^{2}+m^{2})-3m^{4}\\&=3\sigma ^{4}+6\sigma ^{2}m^{2}+m^{4}-12\sigma ^{2}m^{2}-4m^{4}+6\sigma ^{2}m^{2}+6m^{4}-3m^{4}=3\sigma ^{4}\end{aligned}}$

Stąd kurtoza jest równa:

$\operatorname {Kurt} _{E}={\frac {(\mu _{4}(X))}{(\sigma ^{2})^{2}}}-3={\frac {3\sigma ^{4}}{\sigma ^{4}}}-3=3-3=0.$

Zobacz też

• miara rozkładu
• statystyka opisowa