da Davimal » 16/02/2009, 01:49
Quello che Ada dice, nonostante la buona volontà, non è esatto.
Faccio essenzialmente due correzioni:
1) Dato che la distrib. esponenziale non ammette valori negativi per la var. aleatoria è sbagliato nell'integrale usare come estremo di integrazione un valore $1/lambda-x$ che potrebbe portare ad un valore di x negativo, e in quanto tale privo di senso. Semmai, l'integrale va fatto usando come estremi di integrazione 0 e x(99) il che porta (vista la forma della funzione di ripartizione) all'equazione esponenziale in x :
$1-exp(-\lambda x)=0.99$ ovvero
$x_99= 2log(10)/\lambda$
dove log(10)=2.3 circa e $\lambda$, come giustamente indicato da Ada, va stimato come "Media Aritmetica" fatta sul campione.
2) Ma c'è un senso più profondo in cui è sbagliato il suggerimento di Ada. Il suddetto intervallo fiduciale al 99% $(0, x_99)$ NON è un intervallo fiduciale per la media della distribuzione in oggetto. Esso è un intervallo fiduciale al 99% per una successiva (singola) estrazione della variabile aleatoria in esame.
Per trovare un intervallo fiduciale per la media, intanto si calcola la media aritmetica sul campione, che fa m=40. Poi si sfrutta il fatto che per una esponenziale la dev. standard coincide con la media e quindi s=40. (Se si stima anche questa dal campione si ottiene quasi la stessa quantità, cioè 33,5 circa).
Ora se una X ha varianza $\sigma^2$, la media di N osservazioni campionarie di X ha varianza $\sigma_m^2=\sigma^2/N$. Quindi al valore m=40 testè ottenuto dobbiamo attribuire una dev. standard
$\sigma_m=(frac{40^2}{10})^{1/2}=12.65$
A questo punto, nonostante il campione sia di piccola taglia (N=10), possiamo usare comunque in prima approssimazione il teorema centrale del limite per concludere che:
Mentre la X ha distribuzione esponenziale, la media campionaria di X fatta su 10 osservazioni, è approssimativamente distribuita come una gaussiana di valor medio pari a m=40 e di dev.stand. pari a $\sigma_m=12.65$.
Pertanto, l'intervallo fiduciale bilatero al 99% per la media teorica $\mu=1/\lambda$ della variabile originaria X è dato da:
$m-z_{0.495}\sigma_m<\mu<m+z_{0.495}\sigma_m$
cioè, tenendo conto dei valori trovati sopra e del fatto che $z_{0.495}=2.575$:
$40-32.57 < \mu < 40+32.57$ ovvero l'intervallo fiduciale al 99% per $\mu$ è [7.43, 72.57].
Se vi vuole essere più raffinati, si dovrà ricorrere alla T di Student, dato che la varianza viene stimata dal campione e non è nota "in advance".
Ultima modifica di Davimal il 16/02/2009, 10:35, modificato 1 volta in totale.