dove B {\displaystyle \mathrm {B} } è la funzione beta
dove ϝ {\displaystyle \digamma } è la funzione digamma e B {\displaystyle \mathrm {B} } è la funzione beta
dove K n ( x ) {\displaystyle K_{n}(x)} è una funzione di Bessel
Nella teoria delle probabilità la distribuzione di Student, o t di Student, è una distribuzione di probabilità continua che governa il rapporto tra due variabili aleatorie, la prima con distribuzione normale standard e la seconda, al quadrato, segue una distribuzione chi quadrato.
Questa distribuzione interviene nella stima della media di una popolazione che segue la distribuzione normale, e viene utilizzata negli omonimi test t di Student per la significatività e per ogni intervallo di confidenza della differenza tra due medie.
La distribuzione venne descritta nel 1908 da William Sealy Gosset, che pubblicò il suo risultato sotto lo pseudonimo "Student" perché la fabbrica di birra Guinness presso la quale era impiegato vietava ai propri dipendenti di pubblicare articoli affinché questi non divulgassero segreti di produzione. Il nome distribuzione di Student venne successivamente introdotto da Ronald Fisher.[2][3]
La distribuzione di Student con parametro n {\displaystyle n} (gradi di libertà) governa la variabile aleatoria
dove Z {\displaystyle Z} e K {\displaystyle K} sono due variabili aleatorie indipendenti che seguono rispettivamente la distribuzione normale standard N ( 0 , 1 ) {\displaystyle {\mathcal {N}}(0,1)} e la distribuzione chi quadro χ 2 ( n ) {\displaystyle \chi ^{2}(n)} con n {\displaystyle n} gradi di libertà.
La media μ {\displaystyle \mu } e la varianza σ 2 {\displaystyle \sigma ^{2}} di una popolazione X {\displaystyle X} possono essere stimate tramite un suo campione di N {\displaystyle N} elementi, X 1 , … , X N {\displaystyle X_{1},\ldots ,X_{N}} con gli stimatori
Supponiamo che le variabili aleatorie X 1 , … , X N {\displaystyle X_{1},\ldots ,X_{N}} che compongono il campione siano indipendenti e distribuite normalmente, allora X ¯ {\displaystyle {\bar {X}}} è una variabile normale N ( μ , σ 2 N ) {\displaystyle {\mathcal {N}}\left(\mu ,{\frac {\sigma ^{2}}{N}}\right)} con valore atteso μ {\displaystyle \mu } e varianza σ 2 N {\displaystyle {\frac {\sigma ^{2}}{N}}} . Pertanto la variabile Z {\displaystyle Z} così definita
seguirà una distribuzione normale standard, N ( 0 , 1 ) {\displaystyle {\mathcal {N}}(0,1)} . Il problema è che spesso non si conosce σ 2 {\displaystyle \sigma ^{2}} , pertanto dovremo avere a che fare con uno stimatore della varianza come S 2 {\displaystyle S^{2}} .
Dimostreremo che la seguente variabile aleatoria
segue una distribuzione chi-quadro con N − 1 {\displaystyle N-1} gradi di libertà, χ ( N − 1 ) 2 {\displaystyle \chi _{(N-1)}^{2}} .
Le due variabili aleatorie Z {\displaystyle Z} e k {\displaystyle k} sono indipendenti, per il teorema di Cochran.
Pertanto si definisce la variabile aleatoria
Tale variabile aleatoria segue una distribuzione di probabilità detta "t di Student".
Cominciamo con il dimostrare che k {\displaystyle k} è una variabile aleatoria di tipo chi-quadro. Ricordiamo che una distribuzione χ 2 ( n ) {\displaystyle \chi ^{2}(n)} è una particolare variabile di tipo gamma definita come segue
dove Γ ( x ) {\displaystyle \Gamma (x)} è la funzione Gamma di Eulero definita come
Una variabile chi-quadro con n {\displaystyle n} gradi di libertà si ottiene sommando n {\displaystyle n} variabili normali standard N ( 0 , 1 ) {\displaystyle {\mathcal {N}}(0,1)} elevate al quadrato. Detto ciò partiamo dalla definizione della varianza campionaria e aggiungiamo e sottraiamo nell'argomento della sommatoria μ {\displaystyle \mu } , il valore aspettato della variabile aleatoria X i {\displaystyle X_{i}} che coincide con quello della variabile aleatoria X ¯ {\displaystyle {\bar {X}}} .
Definiamo i parametri a i {\displaystyle a_{i}} e b {\displaystyle b} come a i = X i − μ , b = X ¯ − μ {\displaystyle a_{i}=X_{i}-\mu ,b={\bar {X}}-\mu } e riscriviamo la formula precedente
Ora possiamo esplicitare fuori dalle sommatorie tutti i termini che non dipendono da i {\displaystyle i} , ossia X ¯ {\displaystyle {\bar {X}}} e μ {\displaystyle \mu }
sapendo che la somma su tutti gli X i {\displaystyle X_{i}} è uguale a N X ¯ {\displaystyle N{\bar {X}}} . Dividendo ora a destra e a sinistra per σ 2 {\displaystyle \sigma ^{2}} otteniamo a destra delle variabili normali
Abbiamo quindi ottenuto a sinistra una variabile che precedentemente avevamo indicato con k {\displaystyle k} , mentre a destra abbiamo somme di variabili normali standard al quadrato, coincidenti con una variabile chi quadro con N {\displaystyle N} gradi di libertà e un'altra variabile normale anch'essa standard elevata al quadrato, ossia una variabile chi-quadro ad un solo grado di libertà. Sapendo che somme di variabili di tipo chi-quadro con n {\displaystyle n} e m {\displaystyle m} gradi di libertà corrispondono ancora ad una variabile chi-quadro con n + m {\displaystyle n+m} gradi di libertà otteniamo che la funzione di densità di probabilità di k {\displaystyle k} è di tipo chi-quadro con N − 1 {\displaystyle N-1} gradi di libertà.
Pertanto ora iniziamo a dire che
dove n = N − 1 {\displaystyle n=N-1} è il numero di gradi di libertà, e che
Conosciuta la variabile aleatoria k {\displaystyle k} , essa si riduce difatti ad un parametro moltiplicativo per la normale. Dalla definizione di probabilità condizionata si ha
dove
è una distribuzione chi-quadro con n = N − 1 {\displaystyle n=N-1} gradi di libertà. Quindi
Notiamo che la funzione di distribuzione cercata non è altro che una funzione marginale di f ( t n , k ) {\displaystyle f(t_{n},k)} , pertanto si ha
Ponendo una sostituzione con l'argomento dell'esponenziale, mantenendolo però negativo
otteniamo
l'integrale definito ha come risultato la funzione Gamma di Eulero stessa
Pertanto otteniamo al fine il nostro risultato
Notiamo che il limite di questa successione di funzioni per n → ∞ {\displaystyle n\rightarrow \infty } è
Sapendo che il primo limite ha come risultato 1 2 {\displaystyle {\frac {1}{\sqrt {2}}}} e il secondo tende a e − t 2 2 {\displaystyle e^{-{\frac {t^{2}}{2}}}} .
In pratica, prendendo una popolazione di numerosità N {\displaystyle N} molto grande, la variabile aleatoria t tende ad essere una normale standard.
La distribuzione di Student con n {\displaystyle n} gradi di libertà è simmetrica, perché lo è la distribuzione normale standard mentre la distribuzione chi quadrato che funge da "parametro casuale di scala" non produce effetti di distorsione di tale simmetria.
La sua funzione di densità di probabilità è
dove B {\displaystyle \mathrm {B} } la funzione beta.
La sua funzione di ripartizione è
dove I x ( a , b ) = B ( x , a , b ) B ( a , b ) {\displaystyle I_{x}(a,b)={\frac {\mathrm {B} (x,a,b)}{\mathrm {B} (a,b)}}} è la funzione beta incompleta regolarizzata con
Per k < n {\displaystyle k<n} i momenti (semplici o centrali, in quanto coincidono per una PDF simmetrica) di ordine k {\displaystyle k} della distribuzione sono
In particolare, oltre alla speranza matematica E ( t ) = 0 {\displaystyle E(t)=0} e all'indice di asimmetria γ 1 = 0 {\displaystyle \gamma _{1}=0} (per n > 3 {\displaystyle n>3} ) predetti dalla simmetria della distribuzione, si trovano:
Consideriamo infine un ultimo parametro: la larghezza a metà altezza. Per una variabile t {\displaystyle t} di Student abbiamo che il picco della funzione è nel suo valore atteso, ossia in 0 {\displaystyle 0} , dove la distribuzione ha valore massimo
Per cui troviamo i valori di t {\displaystyle t} per i quali f ( t n ) {\displaystyle f(t_{n})} assume altezza uguale a metà della massima assoluta.
Per cui
dove t {\displaystyle t} ha due soluzioni, come ci aspettavamo dalla simmetria della funzione, coincidenti a
Per cui la larghezza a mezza altezza della funzione è data da
Eseguendo il limite per n → ∞ {\displaystyle n\rightarrow \infty } troviamo un'espressione convergente a
che è l'equivalente della larghezza a metà altezza (FWHM) della normale standard. Viceversa per n = 1 {\displaystyle n=1} otteniamo un FWHM = 2. Difatti per n = 1 {\displaystyle n=1} la distribuzione t di Student coincide con una distribuzione di Lorentz-Cauchy di parametri ( 0 , 1 ) {\displaystyle (0,1)} dove la FWHM è per l'appunto uguale a 2 {\displaystyle 2} .
La distribuzione di Student viene utilizzata per definire degli intervalli di confidenza per la media di una popolazione, sulla base degli stimatori puntuali X ¯ {\displaystyle {\bar {X}}} e S n 2 {\displaystyle S_{n}^{2}} della sua media e della sua varianza. Dall'equazione
si ha infatti
Scegliendo quindi dei quantili q α < q β {\displaystyle q_{\alpha }<q_{\beta }} per la distribuzione di Student con n {\displaystyle n} gradi di libertà, si ha
cioè un intervallo di confidenza per la media μ {\displaystyle \mu } con livello di confidenza β − α {\displaystyle \beta -\alpha } è:
Qualora si considerino intervalli simmetrici si può utilizzare l'indice z α {\displaystyle z_{\alpha }} definito da
ossia
e si ottiene l'intervallo di confidenza per μ {\displaystyle \mu } con livello di confidenza α {\displaystyle \alpha }
La distribuzione di Student con parametro n = 1 {\displaystyle n=1} corrisponde alla distribuzione di Cauchy di parametri ( 0 , 1 ) {\displaystyle (0,1)} : entrambe regolano il rapporto X / Y {\displaystyle X/Y} tra due variabili aleatorie indipendenti aventi distribuzione normale standard.
Al tendere di n {\displaystyle n} a infinito la distribuzione di Student con n {\displaystyle n} gradi di libertà converge alla distribuzione normale standard N ( 0 , 1 ) {\displaystyle {\mathcal {N}}(0,1)} .
Se T {\displaystyle T} è una variabile aleatoria con distribuzione t di Student di parametro n {\displaystyle n} , allora F = T 2 {\displaystyle F=T^{2}} segue la distribuzione di Fisher-Snedecor di parametri ( 1 , n ) {\displaystyle (1,n)} .
La seguente tabella[4] esprime, in funzione del parametro n {\displaystyle n} (riga) e di particolari valori di α {\displaystyle \alpha } (colonna), i quantili q α {\displaystyle q_{\alpha }} per la distribuzione di Student di parametro n {\displaystyle n} :
L'ultima riga, indicata con "∞", si riferisce a una distribuzione normale standard.
Altri progetti