Dispensa di modelli lineari in R - Dipartimento di Statistica
Dispensa di modelli lineari in R - Dipartimento di Statistica
Dispensa di modelli lineari in R - Dipartimento di Statistica
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CAPITOLO 3. DISTRIBUZIONI E STUDI DI SIMULAZIONE 33<br />
Density<br />
Histogram of beta2.hat.sim<br />
0 1 2 3 4<br />
2.7 2.9 3.1 3.3<br />
beta2.hat.sim<br />
Figura 3.10: Istogramma e <strong>di</strong>agramma a scatola delle stime <strong>di</strong> β2 ottenute da 1000<br />
campioni simulati.<br />
Il valore è prossimo al livello <strong>di</strong> copertura nom<strong>in</strong>ale 0.95.<br />
Si può fare un grafico <strong>di</strong> alcuni <strong>di</strong> questi <strong>in</strong>tervalli <strong>di</strong> confidenza. Ad esempio, per i<br />
primi 100 <strong>in</strong>tervalli simulati si ha:<br />
> plot(1:100, beta2.ic[1:100, 1], ylim = c(2.5,<br />
+ 3.5), ylab = "beta2")<br />
> po<strong>in</strong>ts(1:100, beta2.ic[1:100, 2])<br />
> segments(1:100, beta2.ic[1:100, 2], 1:100, beta2.ic[1:100,<br />
+ 1])<br />
> abl<strong>in</strong>e(h = 3, col = 2)<br />
Esercizio. Si valuti il livello <strong>di</strong> copertura effettivo dell’<strong>in</strong>tervallo <strong>di</strong> confidenza con<br />
livello nom<strong>in</strong>ale 0.9 per il coefficiente β1. ✸<br />
2.8 2.9 3.0 3.1 3.2 3.3<br />
●<br />
●<br />
●