In een notendop vindt deze techniek een lijn die het beste bij de gegevens past en neemt deze de volgende vorm aan:<\/span><\/p>\n \u0177 = b0<\/sub> + b1<\/sub> x<\/strong><\/span><\/p>\n Goud:<\/span><\/p>\n Deze vergelijking kan ons helpen de relatie tussen de verklarende variabele en de responsvariabele te begrijpen en kan (ervan uitgaande dat deze statistisch significant is) worden gebruikt om de waarde van een responsvariabele te voorspellen gegeven de waarde van de verklarende variabele.<\/span><\/p>\n Deze tutorial biedt stapsgewijze uitleg over het uitvoeren van eenvoudige lineaire regressie in R.<\/span><\/p>\n Voor dit voorbeeld maken we een nep-dataset met de volgende twee variabelen voor 15 studenten:<\/span><\/p>\n We zullen proberen een eenvoudig lineair regressiemodel in te passen met uren<\/em> als verklarende variabele en onderzoeksresultaten<\/em> als responsvariabele.<\/span><\/p>\n De volgende code laat zien hoe u deze nep-dataset in R kunt maken:<\/span><\/p>\n Voordat we een eenvoudig lineair regressiemodel kunnen toepassen, moeten we eerst de gegevens visualiseren om deze te begrijpen.<\/span><\/p>\n Ten eerste willen we er zeker van zijn dat de relatie tussen uren<\/em> en score<\/em> ongeveer lineair is, aangezien dit een enorme onderliggende aanname<\/a> is van eenvoudige lineaire regressie. We kunnen een eenvoudig spreidingsdiagram maken om de relatie tussen de twee variabelen te visualiseren:<\/span> <\/p>\n Uit de grafiek kunnen we zien dat de relatie lineair lijkt te zijn. Naarmate het aantal uren<\/em> toeneemt, stijgt de score<\/em> ook lineair.<\/span><\/p>\n Vervolgens kunnen we een boxplot maken om de verdeling van de examenresultaten te visualiseren en te controleren op uitschieters<\/a> . Standaard definieert R een waarneming als een uitbijter als deze 1,5 keer de interkwartielafstand boven het derde kwartiel (Q3) of 1,5 keer de interkwartielafstand onder het eerste kwartiel (Q1) ligt.<\/span><\/p>\n Als een waarneming een uitbijter is, verschijnt er een kleine cirkel in de boxplot:<\/span> <\/p>\n Er zijn geen kleine cirkels in de boxplot, wat betekent dat er geen uitschieters in onze dataset voorkomen.<\/span><\/p>\n Zodra we hebben bevestigd dat de relatie tussen onze variabelen lineair is en er geen uitschieters zijn, kunnen we overgaan tot het fitten van een eenvoudig lineair regressiemodel met uren<\/em> als verklarende variabele en de score<\/em> als responsvariabele:<\/span><\/p>\n Uit de samenvatting van het model kunnen we zien dat de aangepaste regressievergelijking is:<\/span><\/p>\n Score = 65.334 + 1.982*(uur)<\/strong><\/span><\/p>\n Dit betekent dat elk extra bestudeerd uur gepaard gaat met een gemiddelde stijging van de examenscore van 1.982<\/strong> punten. En de oorspronkelijke waarde van 65.334<\/strong> vertelt ons de gemiddelde verwachte examenscore voor een student die nul uur studeert.<\/span><\/p>\n We kunnen deze vergelijking ook gebruiken om de verwachte examenscore te vinden op basis van het aantal uren dat een student studeert. Een student die bijvoorbeeld 10 uur studeert, moet een examenscore van 85,15<\/strong> behalen:<\/span><\/p>\n Score = 65,334 + 1,982*(10) = 85,15<\/strong><\/span><\/p>\n Zo interpreteert u de rest van de modelsamenvatting:<\/span><\/p>\n Nadat het eenvoudige lineaire regressiemodel aan de gegevens is aangepast, is de laatste stap het maken van residuele plots.<\/span><\/p>\n Een van de belangrijkste aannames van lineaire regressie is dat de residuen van een regressiemodel bij benadering normaal verdeeld zijn en homoscedastisch<\/a> zijn op elk niveau van de verklarende variabele. Als niet aan deze aannames wordt voldaan, kunnen de resultaten van ons regressiemodel misleidend of onbetrouwbaar zijn.<\/span><\/p>\n Om te verifi\u00ebren dat aan deze aannames wordt voldaan, kunnen we de volgende restplots maken:<\/span><\/p>\n Plot van residuen versus aangepaste waarden:<\/strong> deze plot is nuttig voor het bevestigen van homoscedasticiteit. Op de x-as worden de gefitte waarden weergegeven en op de y-as de residuen. Zolang de residuen willekeurig en uniform verspreid lijken te zijn over de grafiek rond de nulwaarde, kunnen we aannemen dat de homoscedasticiteit niet wordt geschonden:<\/span> <\/p>\n\n
Stap 1: Gegevens laden<\/b><\/span><\/h3>\n
\n
#create dataset<\/span>\ndf <- data.frame(hours=c(1, 2, 4, 5, 5, 6, 6, 7, 8, 10, 11, 11, 12, 12, 14),\n score=c(64, 66, 76, 73, 74, 81, 83, 82, 80, 88, 84, 82, 91, 93, 89))\n\n#view first six rows of dataset\n<\/span>head(df)\n\n hours score\n1 1 64\n2 2 66\n3 4 76\n4 5 73\n5 5 74\n6 6 81\n\n#attach dataset to make it more convenient to work with\n<\/span>attach(df)\n<\/strong><\/pre>\n Stap 2: Visualiseer de gegevens<\/b><\/span><\/h3>\n
scatter.smooth(hours, score, main=' Hours studied vs. Exam Score<\/span> ')\n<\/strong><\/pre>\n![\"Spreidingsplot](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/simpleregr1.png\")
<\/p>\n boxplot(score)<\/strong> <\/pre>\n
![\"Boxplot](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/simpleregr2.png\")
<\/p>\n Stap 3: Voer een eenvoudige lineaire regressie uit<\/b><\/span><\/h3>\n
#fit simple linear regression model\n<\/span>model <- lm(score~hours)\n\n#view model summary<\/span>\nsummary(model)\n\nCall:\nlm(formula = score ~ hours)\n\nResiduals:\n Min 1Q Median 3Q Max \n-5,140 -3,219 -1,193 2,816 5,772 \n\nCoefficients:\n Estimate Std. Error t value Pr(>|t|) \n(Intercept) 65,334 2,106 31,023 1.41e-13 ***\nhours 1.982 0.248 7.995 2.25e-06 ***\n---\nSignificant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n\nResidual standard error: 3.641 on 13 degrees of freedom\nMultiple R-squared: 0.831, Adjusted R-squared: 0.818 \nF-statistic: 63.91 on 1 and 13 DF, p-value: 2.253e-06\n<\/strong><\/pre>\n\n
Stap 4: Maak resterende plots<\/strong><\/span><\/h3>\n
#define residuals\n<\/span>res <- resid(model)\n\n#produce residual vs. fitted plot\n<\/span>plot(fitted(model), res)\n\n#add a horizontal line at 0 \n<\/span>abline(0,0)\n<\/strong><\/pre>\n![\"Residuele](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/simpleregr3.png\")
<\/p>\n
![\"Normale](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/simpleregr4.png\")
<\/p>\n