Distribució khi quadrat no central

Plantilla:Distribució de probabilitat En Teoria de la Probabilitat i Estadística, la distribució khi quadrat no central (o distribució ${\chi }^2$ no central) és una generalització de la distribució khi quadrat incorporant un paràmetre que s'anomena de no centrament. Sovint sorgeix en l'anàlisi de potència de contrast d'hipòtesis estadístiques en què la distribució nul·la és (potser asimtòticament) una distribució khi quadrat; exemples importants d'aquestes proves són les prova de raó de versemblança.

Definicions

Com el el cas de la distribució ${\chi }^2$ordinària començarem pel cas que el nombre de graus de llibertat sigui un nombre enter positiu i després, mitjançant la funció de densitat ho estendrem a qualsevol nombre de graus de llibertat $k>0$ .Siguin $X_1,\dots ,X_k$ variables aleatòries independents, distribuïdes normalment amb mitjanes ${\mu }_1,\dots ,{\mu }_k$ respectivament i totes amb variància 1: $X_j\sim N({\mu }_j,1),j=1,\dots ,k$. Aleshores es diu que la variable aleatòria

$$J=\sum _{i=1}^k{X}_i^2$$

té una distribució khi-quadrat no central amb $k$ graus de llibertat i paràmetre de no centralitat $\lambda =\sum _{i=1}^k{\mu }_i^2.$ S'escriu $J\sim {\chi }_k^2(\lambda )$. Si $\lambda =0$ , aleshores $J$ té una distribució ${\chi }^2$ ordinària amb $k$ graus de llibertat: $X\sim {\chi }_k^2$ .

Equivalentment, es pot definir la distribució ${\chi }_k^2(\lambda )$ com la distribució de la suma $$(Z_1+{\mu }_1{)}^2+\cdots +(Z_k+{\mu }_k{)}^2,$$on $Z_1,\dots ,Z_k$ són variables aleatòries independents, totes amb distribució normal estàndard $𝒩(0,1)$ .

Nota: Algunes referències defineixen el paràmetre de no centralitat d'altres maneres, com la meitat de la suma $\sum _{i=1}^k{\mu }_i^2$ o la seva arrel quadrada.

Funció de densitat

La funció de densitat de probabilitat (pdf) ve donada per

$$f_J(x;k,\lambda )=\sum _{j=0}^{\mathrm{\infty }}\frac{e^{-\lambda /2}(\lambda /2{)}^j}{j!}{f}_{Y_{k+2j}}(x),x>0,(*)$$

on $Y_q$ es distribueix com una ${\chi }^2$amb $q$ graus de llibertat, $Y_q\sim {\chi }_q^2$ i $f_{Y_q}(x)$ és l seva funció de densitat: $$f_{Y_q}(x)=\frac{x^{q/2-1}{e}^{-x/2}}{2^{q/2}\mathrm{\Gamma }(q/2)},x>0,$$on $\mathrm{\Gamma }(u)$ és la funció gamma d'Euler.

És a dir, la distribució ${\chi }_k^2(\lambda )$ és una mixtura de distribucions ${\chi }_{k+2j}^2,j=0,1,2,\dots$ , amb pesos donats per una distribució de Poisson de paràmetre $\lambda /2$ .

Expressió alternativa de la funció de densitat

La funció de densitat també es pot escriure $$f_J(x)=\frac{1}{2}{e}^{-(x+\lambda )/2}(\frac{x}{\lambda }{)}^{k/4-1/2}{I}_{k/2-1}(\sqrt{\lambda x}),x>0,$$on $I_{\alpha }(y)$ és la funció de Bessel modificada de primer tipus, $$I_{\alpha }(y)=(\frac{y}{2}{)}^{\alpha }\sum _{j=0}^{\mathrm{\infty }}\frac{(y^2/4{)}^j}{j!\mathrm{\Gamma }(\alpha +j+1)}.$$

Extensió a un nombre de graus de llibertat no enter

La funció (*) està ben definida i és una funció de densitat per a qualsevol $k>0$ . Per tant, podem definir una variable ${\chi }_k^2(\delta )$ amb $k>0$ com aquella que té per funció de densitat (*). Naturalment, es perd la interpretació com al nombre de sumands independents.

Moments, funció generatriu de moments i funció característica

Els moments es poden calcular utilitzant les propietats de les mixtures de distribucions. Sigui $J\sim {\chi }_k^2(\lambda )$ i $Y_q\sim {\chi }^2(q)$. Designem per $p_j$ els pesos donats per una distribució de Poisson de paràmetre $\lambda /2$: $$p_j=e^{-\lambda /2}\frac{(\lambda /2{)}^j}{j!},j=0,1,\dots$$Aleshores, atès que una distribució khi-quadrat té moments de tots els ordres, tindrem que la distribució khi-quadrat no central també, i $$E\left[J^n\right]=\sum _{j=0}^{\mathrm{\infty }}{p}_{j}E\left[Y_{k+2j}^n\right].$$Per exemple, per a $n=1$, $E[Y_q]=q$ i llavors $$E[J]=\sum _{j=0}^{\mathrm{\infty }}{p}_j(k+2j)=k\sum _{j=0}^{\mathrm{\infty }}{p}_j+2\sum _{j=0}^{\mathrm{\infty }}j{p}_j=k+2\frac{\lambda }{2}=k+\lambda .$$De manera anàloga, es calcula $$E[J^2]=(k+\lambda {)}^2+2(k+2\lambda ),$$d'on $$\text{Var}(J)=2(k+2\lambda ).$$ La funció generatriu de moments també es pot calcular de la mateixa forma: Designem per $M_Y$ la funció generatriu de moments d'una variable aleatòria $Y$, $$M_Y(t)=E[e^{tY}].$$ Si $Y_q\sim {\chi }^2(q)$, $$M_{Y_q}(t)=\frac{1}{(1-2t{)}^{q/2}},t\in (-\mathrm{\infty },1/2).$$Llavors, per a $t\in (-\mathrm{\infty },1/2)$,$$M_J(t)=\sum _{j=0}^{\mathrm{\infty }}{p}_j{M}_{k+2j}(t)=e^{-\lambda /2}\sum _{j=0}^{\mathrm{\infty }}\frac{(\lambda /2{)}^j}{j!}\frac{1}{(1-2t{)}^{k/2+j}}=\frac{e^{-\lambda /2}}{(1-2t{)}^{k/2}}\sum _{j=0}^{\mathrm{\infty }}\frac{1}{j!}(\frac{\lambda }{2(1-2t)}{)}^j=\frac{e^{\lambda t/(1-2t)}}{(1-2t{)}^{k/2}}.$$Anàlogament, la funció característica dona $${\phi }_J(t)=\frac{e^{i\lambda t/(1-2it)}}{(1-2it{)}^{k/2}}.$$

Una propietat de les formes quadràtiques en variables normals

Aquesta propietat, que té interès per ella mateixa, és el fonament de la utilització de la distribució ${\chi }_k^2(\lambda )$ en l'estudi de la potència d'un test sobre la mitjana d'una població normal multivariable, segon veurem en un exemple.

Propietat. Considerem un vector aleatori normal multivariable $𝑿\sim {𝒩}_k(\mathit{\mu },\mathrm{\Sigma })$ amb $\text{det}(\mathrm{\Sigma })>0$. Aleshores:

1. $\mathit{(}𝑿\mathit{-}\mathit{\mu }{\mathit{)}}^{\prime }{\mathrm{\Sigma }}^{\mathit{-}1}\mathit{(}𝑿\mathit{-}\mathit{\mu }\mathit{)}\sim {\chi }_k^2$.
2. ${𝑿}^{\prime }{\mathrm{\Sigma }}^{\mathit{-}1}𝑿\sim {\chi }_k^2(\lambda )$, amb $\lambda ={\mathit{\mu }}^{\prime }{\mathrm{\Sigma }}^{\mathit{-}1}\mathit{\mu }$.

Plantilla:Caixa desplegable Exemple. Muirhead. Considerem un contrast d'hipòtesis sobre la mitjana d'una població normal multivariable amb matriu de variàncies-covariàncies coneguda. Sigui ${𝑿}_1,\dots ,{𝑿}_n$ una mostra d'una distribució ${𝒩}_k(\mathit{\mu },\mathrm{\Sigma })$. Llavors $$\overline{𝑿}=\frac{1}{n}\sum _{i=1}^n{𝑿}_i\sim {𝒩}_k(\mathit{\mu },\frac{1}{n}\mathrm{\Sigma }).$$Fixem ${\mathit{\mu }}_0\in {ℝ}^k$. Anem a fer el contrast $${\text{H}}_0:\mathit{\mu }={\mathit{\mu }}_0\text{contra}{\text{H}}_1:\mathit{\mu }\ne {\mathit{\mu }}_0.$$Com a estadístic de contrast utilitzarem $$W=\mathit{(}𝑿\mathit{-}{\mathit{\mu }}_0{\mathit{)}}^{\prime }{\mathrm{\Sigma }}^{\mathit{-}1}\mathit{(}𝑿\mathit{-}{\mathit{\mu }}_0\mathit{)}.$$Fixem un nivell de significació del test $\alpha .$ Atès que, per la primera part de la propietat anterior, sota ${\text{H}}_0$ , $W\sim {\chi }_k^2$ , rebutjarem ${\text{H}}_0$ si $$W>C_{\alpha },$$on $C_{\alpha }$ és el nombre tal que $$P({\chi }_k^2>C_{\alpha })=\alpha .$$Si ${\text{H}}_0$ no és veritat, $$\overline{𝑿}-{\mathit{\mu }}_0\sim {𝒩}_k(\mathit{\mu }-{\mathit{\mu }}_0,\frac{1}{n}\mathrm{\Sigma }).$$Per tant, per la segona part de la propietat anterior, $$W\sim {\chi }_k^2(\lambda ),\text{amb}\lambda =n(\mathit{\mu }-{\mathit{\mu }}_0{)}^{\prime }{\mathrm{\Sigma }}^{-1}(\mathit{\mu }-{\mathit{\mu }}_0).$$Per tant, la potència del test és funció de $\lambda$: $$\beta (\lambda )=P_{\lambda }(W>C)=P({\chi }_k^2(\lambda )>C_{\alpha }).$$

Ocurrència i aplicacions

Es poden obtenir intervals de tolerància de regressió normal a dues cares basant-se en la distribució khi quadrat no central. Això permet calcular un interval estadístic dins del qual, amb un cert nivell de confiança, es troba una proporció especificada d'una població mostrada.