In deze zelfstudie wordt uitgelegd hoe u een betrouwbaarheidsinterval voor een gegevensset in R kunt plotten.<\/span><\/p>\n Stel dat we de volgende dataset in R hebben met 100 rijen en 2 kolommen:<\/span><\/p>\n Om een grafiek te maken van de relatie tussen x en y, kunnen we eerst een lineair regressiemodel passen:<\/span><\/p>\n Vervolgens kunnen we een grafiek maken van de geschatte lineaire regressielijn met behulp van de functie abline() en de functie lines() om de feitelijke betrouwbaarheidsbanden te cre\u00ebren:<\/span> <\/p>\n De zwarte lijn toont de gepaste lineaire regressielijn, terwijl de twee gestreepte blauwe lijnen de betrouwbaarheidsintervallen weergeven.<\/span><\/p>\n Optioneel kunt u ook het gebied tussen de betrouwbaarheidsintervallijnen en de geschatte lineaire regressielijn opvullen met behulp van de volgende code:<\/span> <\/p>\n Voorbeeld: een betrouwbaarheidsinterval in R uitzetten<\/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![\"Een](\"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![\"Betrouwbaarheidsintervalgrafiek](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/plotconfintr2.png\")
<\/p>\n