Questo tutorial spiega come tracciare un intervallo di confidenza per un set di dati in R.<\/span><\/p>\n Supponiamo di avere il seguente set di dati in R con 100 righe e 2 colonne:<\/span><\/p>\n Per creare un grafico della relazione tra x e y, possiamo prima adattare un modello di regressione lineare:<\/span><\/p>\n Successivamente, possiamo creare un grafico della retta di regressione lineare stimata utilizzando la funzione abline() e la funzione lines() per creare le bande di confidenza effettive:<\/span> <\/p>\n La linea nera mostra la retta di regressione lineare adattata mentre le due linee blu tratteggiate mostrano gli intervalli di confidenza.<\/span><\/p>\n Facoltativamente, \u00e8 anche possibile riempire l’area tra le linee dell’intervallo di confidenza e la linea di regressione lineare stimata utilizzando il seguente codice:<\/span> <\/p>\n Esempio: tracciare un intervallo di confidenza in R<\/strong><\/h3>\n
#make this example reproducible<\/span>\nset.seed(0)\n\n#create dataset<\/span>\nx <- rnorm(100)\ny <- x*2 + rnorm(100)\ndf <- data.frame(x = x, y = y)\n\n#view first six rows of dataset<\/span>\nhead(df)\n\n xy\n1 1.2629543 3.3077678\n2 -0.3262334 -1.4292433\n3 1.3297993 2.0436086\n4 1.2724293 2.5914389\n5 0.4146414 -0.3011029\n6 -1.5399500 -2.5031813\n<\/strong><\/pre>\n model <- lm(y ~ x, data = df)\n<\/strong><\/pre>\n
#get predicted y values using regression equation<\/span>\nnewx <- seq(min(df$x), max(df$x), length.out=100)\npreds <- predict(model, newdata = data.frame(x=newx), interval = 'confidence')\n\n#create plot of x vs. y, but don't display individual points (type='n')<\/span> \nplot(y ~ x, data = df, type = 'n')\n\n#add fitted regression line<\/span>\nabline(model)\n\n#add dashed lines for confidence bands\n<\/span>lines(newx, preds[,3], lty = 'dashed', col = 'blue')\nlines(newx, preds[,2], lty = 'dashed', col = 'blue')<\/strong> <\/pre>\n![\"Tracciare](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/plotconfintr1.png\")
<\/p>\n #create plot of x vs. y<\/span>\nplot(y ~ x, data = df, type = 'n')\n\n#fill in area between regression line and confidence interval\npolygon(c(rev(newx), newx), c(rev(preds[,3]), preds[,2]), col = 'grey', border = NA)<\/span>\n\n#add fitted regression line<\/span>\nabline(model)\n\n#add dashed lines for confidence bands\n<\/span>lines(newx, preds[,3], lty = 'dashed', col = 'blue')\nlines(newx, preds[,2], lty = 'dashed', col = 'blue')<\/strong> <\/pre>\n![\"Grafico](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/plotconfintr2.png\")
<\/p>\n