20.06.2013 Views

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

SHOW MORE
SHOW LESS

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 />

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!