{"id":484,"date":"2023-07-29T18:17:31","date_gmt":"2023-07-29T18:17:31","guid":{"rendered":"https:\/\/statorials.org\/nl\/ancova-in-r\/"},"modified":"2023-07-29T18:17:31","modified_gmt":"2023-07-29T18:17:31","slug":"ancova-in-r","status":"publish","type":"post","link":"https:\/\/statorials.org\/nl\/ancova-in-r\/","title":{"rendered":"Hoe voer je een ancova uit in r"},"content":{"rendered":"

<\/p>\n

\n

Deze tutorial geeft een voorbeeld van hoe u een ANCOVA<\/a> uitvoert in R.<\/span><\/p>\n

Voorbeeld: ANCOVA in R<\/strong><\/span><\/h3>\n

We voeren een ANCOVA uit om te testen of de studietechniek al dan niet invloed heeft op de examenresultaten, met behulp van de volgende variabelen:<\/span><\/p>\n

    \n
  • Technisch onderzoek<\/strong> : de onafhankelijke variabele die we willen analyseren<\/span><\/li>\n
  • Huidig leerlingcijfer<\/strong> : de covariabele waarmee we rekening willen houden<\/span><\/li>\n
  • Beoordelingsscore<\/strong> : de responsvariabelen die we willen analyseren<\/span><\/li>\n<\/ul>\n

    De volgende dataset bevat informatie over 90 studenten, willekeurig verdeeld in drie groepen van 30.<\/span><\/p>\n

    De dataset toont de studietechniek die elke leerling gebruikte (A, B<\/em> of C)<\/em> ,<\/em> hun huidige cijfer in de klas toen ze begonnen met het gebruik van de studietechniek, en het cijfer dat ze kregen op het examen nadat ze de studietechniek een maand lang ter voorbereiding hadden gebruikt. voor het examen. examen:<\/span><\/p>\n

     #make this example reproducible<\/span>\nset.seed(10)\n\n#create dataset<\/span>\ndata <- data.frame(technique = rep(c(\"A\", \"B\", \"C\"), each = 30),\n                   current_grade = runif(90, 65, 95),\n                   exam = c(runif(30, 80, 95), runif(30, 70, 95), runif(30, 70, 90)))\n\n#view first six lines of dataset<\/span>\nhead(data)\n\n# technical current_grade exam\n#1 A 80.22435 87.32759\n#2 A 74.20306 90.67114\n#3 A 77.80723 88.87902\n#4 A 85.79306 87.75735\n#5 A 67.55408 85.72442\n#6 A 71.76310 92.52167\n<\/strong><\/pre>\n
     Stap 1: Verken de gegevens<\/strong><\/span><\/h3>\n
     Voordat we het ANCOVA-model passen, moeten we eerst de gegevens verkennen om deze beter te begrijpen en te verifi\u00ebren dat er geen extreme uitschieters zijn die de resultaten zouden kunnen vertekenen.<\/span><\/p>\n
     Ten eerste kunnen we een samenvatting van elke variabele in de dataset weergeven:<\/span><\/p>\n
     summary(data)\n\n# technical current_grade exam      \n#A:30 Min. :65.43 Min. :71.17  \n# B:30 1st Qu.:71.79 1st Qu.:77.27  \n# C:30 Median:77.84 Median:84.69  \n# Mean:78.15 Mean:83.38  \n# 3rd Qu.:83.65 3rd Qu.:89.22  \n# Max. :93.84 Max. :94.76  \n<\/strong><\/pre>\n
     We kunnen zien dat elke studietechniekwaarde ( A, B<\/em> en C)<\/em> 30 keer in de gegevens voorkomt.<\/span><\/p>\n
     Ook kunnen we zien hoe de huidige scores van studenten bij aanvang van de studie verdeeld waren. De minimale score in de klas was 65,43, de maximale score 93,84 en het gemiddelde was 78,15.<\/span><\/p>\n
     Op dezelfde manier kunnen we zien dat de minimale score behaald op het examen 71,17 was, de maximale score 94,76 en het gemiddelde 83,38.<\/span><\/p>\n
     Dan kunnen we het dplyr-<\/strong> pakket gebruiken om eenvoudig de gemiddelde en standaardafwijking van huidige cijfers en examenresultaten voor elke studietechniek te vinden:<\/span><\/p>\n
     #load dplyr<\/em><\/span>\nlibrary(dplyr)\n\ndata %>%\n  group_by<\/span> (technical) %>%\n  summarize<\/span> (mean_grade = mean(current_grade),\n            sd_grade = sd(current_grade),\n            mean_exam = mean(exam),\n            sd_exam = sd(exam))\n\n# A tibble: 3 x 5\n# technique mean_grade sd_grade mean_exam sd_exam                      \n#1 A 79.0 7.00 88.5 3.88\n#2 B 78.5 8.33 81.8 7.62\n#3 C 76.9 8.24 79.9 5.71<\/strong><\/pre>\n
     We kunnen zien dat de gemiddelde en standaardafwijkingen van het huidige leerjaar van studenten die elke studietechniek gebruiken ongeveer vergelijkbaar zijn.<\/span><\/p>\n
     Ook zien we dat de gemiddelde examenscore significant hoger is voor studenten die gebruik maken van studietechniek A<\/em> vergeleken met technieken B<\/em> en C.<\/em><\/span><\/p>\n
     We kunnen de verdeling van examenresultaten ook visualiseren op basis van studietechniek met behulp van boxplots :<\/span><\/p>\n
     boxplot(exam ~ technique,\ndata = data,\nmain = \"Exam Score by Studying Technique\",\nxlab = \"Studying Technique\",\nylab = \"Exam Score\",\ncol = \"steelblue\",\nborder = \"black\"\n)<\/strong><\/pre>\n
    <\/h3>\n
     Op dezelfde manier kunnen we ook boxplots gebruiken om de verdeling van de huidige cijfers<\/em> te visualiseren op basis van de studietechniek:<\/span><\/p>\n
     boxplot(current_grade ~ technical,\ndata = data,\nmain = \"Current Grade by Studying Technique\",\nxlab = \"Studying Technique\",\nylab = \"Current Grade\",\ncol = \"steelblue\",\nborder = \"black\"\n)<\/strong><\/pre>\n
     Stap 2: Controleer de modelaannames<\/strong><\/span><\/h3>\n
     Zodra we wat basisgegevensverkenning hebben uitgevoerd en bekend zijn met de gegevens, moeten we verifi\u00ebren dat aan de volgende aannames voor ANCOVA wordt voldaan:<\/span><\/p>\n
      \n
    •  De covariabele en de behandeling zijn onafhankelijk<\/strong> \u2013 het is noodzakelijk om te verifi\u00ebren dat de covariabele ( huidige graad)<\/em> en de behandeling (onderzoekstechniek)<\/em> onafhankelijk van elkaar zijn, omdat het toevoegen van een covariabele term in het model alleen zinvol is als de covariabele en de behandeling de behandeling werkt onafhankelijk op de responsvariabele ( onderzoek<\/em> ).<\/span><\/li>\n
    •  Homogeniteit van variantie<\/strong> \u2013 we moeten controleren of de varianties tussen groepen gelijk zijn<\/span><\/li>\n<\/ul>\n
       Om te verifi\u00ebren dat de covariabele en de behandeling onafhankelijk zijn, kunnen we een ANOVA uitvoeren met de huidige graad<\/em> als responsvariabele en de onderzoekstechniek<\/em> als voorspellende variabele:<\/span><\/p>\n
       #fit anova model\n<\/span>anova_model <- aov(current_grade ~ technique, data = data)\n#view summary of anova model\n<\/span>summary(anova_model)\n\n# Df Sum Sq Mean Sq F value Pr(>F)\n#technical 2 74 37.21 0.599 0.552\n#Residuals 87 5406 62.14    \n<\/strong><\/pre>\n
       De p-waarde is groter dan 0,05, dus de covariabele ( huidige graad)<\/em> en behandeling ( onderzoekstechniek<\/em> ) lijken onafhankelijk.<\/span><\/p>\n
       Om vervolgens te controleren of er homogeniteit van variantie tussen groepen bestaat, kunnen we de Levene-test uitvoeren:<\/span><\/p>\n
       #load car<\/em> library to conduct Levene's Test<\/span>\nlibary(car)\n\n#conduct Levene's Test<\/span>\nleveneTest(technical exam, data = data)\n\n#Levene's Test for Homogeneity of Variance (center = median)\n# Df F value Pr(>F)    \n#group 2 9.4324 0.0001961 ***\n#87   \n<\/strong><\/pre>\n
       De p-waarde van de test is gelijk aan 0,0001961, wat aangeeft dat de varianties tussen de groepen niet gelijk zijn. Hoewel we een transformatie van de gegevens zouden kunnen proberen om dit probleem te corrigeren, zullen we ons op dit moment niet al te veel zorgen maken over variantieverschillen.<\/span><\/p>\n
       Stap 3: Pas het ANCOVA-model aan<\/strong><\/span><\/h3>\n
       Vervolgens passen we het ANCOVA-model aan met behulp van de examenscore<\/em> als de responsvariabele, de onderzoekstechniek<\/em> als de voorspellende (of \u2018behandelings\u2019) variabele en het huidige cijfer<\/em> als de covariabele.<\/span><\/p>\n
       We zullen de functie Anova() in het pakket gebruiken omdat we dit doen, gewoon om te kunnen specificeren dat we de kwadratensom van type III voor het model willen gebruiken, omdat de kwadratensom van type I afhangt van de volgorde waarin de Predictors worden in het model ingevoerd:<\/span><\/p>\n
       #load car<\/em> library<\/span>\nlibrary(car)\n\n#fit ANCOVA model\n<\/span>ancova_model <- aov(exam ~ technique + current_grade, data = data)\n\n#view summary of model\n<\/span>Anova(ancova_model, type=\"III\") \n\n#Answer: exam\n# Sum Sq Df F value Pr(>F)    \n#(Intercept) 7161.2 1 201.4621 < 2.2e-16 ***\n#technical 1242.9 2 17.4830 4.255e-07 ***\n#current_grade 12.3 1 0.3467 0.5576    \n#Residuals 3057.0 86         \n<\/strong><\/pre>\n
       We kunnen zien dat de p-waarde voor techniek<\/em> extreem laag is, wat aangeeft dat studietechniek een statistisch significant effect heeft op examenscores, zelfs na controle voor het huidige cijfer.<\/span><\/p>\n
       Stap 4: Post-hoc testen<\/strong><\/span><\/h3>\n
       Hoewel de ANCOVA-resultaten ons vertelden dat de studietechniek<\/em> een statistisch significant effect had op de examenscores, moeten we post-hoctesten uitvoeren om te bepalen welke studietechnieken van elkaar verschillen.<\/span><\/p>\n
       Om dit te doen, kunnen we de functie glht() in het multcomp-<\/strong> pakket in R gebruiken om de Tukey-test uit te voeren voor meerdere vergelijkingen:<\/span><\/p>\n
       #load the multcomp<\/em> library<\/span>\nlibrary(multicomp)\n\n#fit the ANCOVA model\n<\/span>ancova_model <- aov(exam ~ technique + current_grade, data = data)\n\n#define the post hoc comparisons to make\n<\/span>postHocs <- glht(ancova_model, linfct = mcp(technique = \"Tukey\"))\n\n#view a summary of the post hoc comparisons\n<\/span>summary(postHocs)\n\n#Multiple Comparisons of Means: Tukey Contrasts\n#\n#Fit: aov(formula = exam ~ technique + current_grade, data = data)\n#\n#Linear Assumptions:\n#Estimate Std. Error t value Pr(>|t|)    \n#B - A == 0 -6.711 1.540 -4.358 0.000109 ***\n#C - A == 0 -8.736 1.549 -5.640 < 1e-04 ***\n#C - B == 0 -2.025 1.545 -1.311 0.393089    \n\n#view the confidence intervals associated with the multiple comparisons\n<\/span>confint(postHocs)\n\n# Simultaneous Confidence Intervals\n#\n#Multiple Comparisons of Means: Tukey Contrasts\n#\n#Fit: aov(formula = exam ~ technique + current_grade, data = data)\n#\n#Quantile = 2.3845\n#95% family-wise confidence level\n#\n#Linear Assumptions:\n# Estimate lwr upr     \n#B - A == 0 -6.7112 -10.3832 -3.0392\n#C - A == 0 -8.7364 -12.4302 -5.0426\n#C - B == 0 -2.0252 -5.7091 1.6588\n<\/strong><\/pre>\n
       Uit het resultaat kunnen we zien dat er een statistisch significant verschil is (bij \u03b1 = 0,05) in de examenresultaten tussen het bestuderen van techniek A<\/em> en het bestuderen van techniek B<\/em> (p-waarde: 0,000109) en tussen techniek A<\/em> en techniek C<\/em> (p-waarde: <1e-04).<\/span><\/p>\n
       We kunnen ook zien dat er geen<\/em> statistisch significant verschil is (bij \u03b1 = 0,05) tussen technieken B<\/em> en C.<\/em><\/span> Betrouwbaarheidsintervallen tussen technieken ondersteunen deze bevindingen ook.<\/span><\/p>\n
       We kunnen dus concluderen dat het gebruik van studietechniek A<\/em> leidt tot een statistisch significant hoger examencijfer voor studenten vergeleken met technieken B<\/em> en C<\/em> , zelfs na controle voor het huidige cijfer van de student in de klas.<\/span><\/p>\n","protected":false},"excerpt":{"rendered":"
      Deze tutorial geeft een voorbeeld van hoe u een ANCOVA uitvoert in R. Voorbeeld: ANCOVA in R We voeren een ANCOVA uit om te testen of de studietechniek al dan niet invloed heeft op de examenresultaten, met behulp van de volgende variabelen: Technisch onderzoek : de onafhankelijke variabele die we willen analyseren Huidig leerlingcijfer : […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[11],"tags":[],"class_list":["post-484","post","type-post","status-publish","format-standard","hentry","category-gids"],"yoast_head":"\nHoe een ANCOVA uit te voeren in R - Statorials<\/title>\n<meta name=\"description\" content=\"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.\" \/>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/statorials.org\/nl\/ancova-in-r\/\" \/>\n<meta property=\"og:locale\" content=\"de_DE\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Hoe een ANCOVA uit te voeren in R - Statorials\" \/>\n<meta property=\"og:description\" content=\"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.\" \/>\n<meta property=\"og:url\" content=\"https:\/\/statorials.org\/nl\/ancova-in-r\/\" \/>\n<meta property=\"og:site_name\" content=\"Statorials\" \/>\n<meta property=\"article:published_time\" content=\"2023-07-29T18:17:31+00:00\" \/>\n<meta name=\"author\" content=\"Dr.benjamin anderson\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Verfasst von\" \/>\n\t<meta name=\"twitter:data1\" content=\"Dr.benjamin anderson\" \/>\n\t<meta name=\"twitter:label2\" content=\"Gesch\u00e4tzte Lesezeit\" \/>\n\t<meta name=\"twitter:data2\" content=\"6\u00a0Minuten\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\/\/schema.org\",\"@graph\":[{\"@type\":\"WebPage\",\"@id\":\"https:\/\/statorials.org\/nl\/ancova-in-r\/\",\"url\":\"https:\/\/statorials.org\/nl\/ancova-in-r\/\",\"name\":\"Hoe een ANCOVA uit te voeren in R - Statorials\",\"isPartOf\":{\"@id\":\"https:\/\/statorials.org\/nl\/#website\"},\"datePublished\":\"2023-07-29T18:17:31+00:00\",\"dateModified\":\"2023-07-29T18:17:31+00:00\",\"author\":{\"@id\":\"https:\/\/statorials.org\/nl\/#\/schema\/person\/d4b8842173cca1bb62cdec41860e4219\"},\"description\":\"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.\",\"breadcrumb\":{\"@id\":\"https:\/\/statorials.org\/nl\/ancova-in-r\/#breadcrumb\"},\"inLanguage\":\"de\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\/\/statorials.org\/nl\/ancova-in-r\/\"]}]},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\/\/statorials.org\/nl\/ancova-in-r\/#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Thuis\",\"item\":\"https:\/\/statorials.org\/nl\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"Hoe voer je een ancova uit in r\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\/\/statorials.org\/nl\/#website\",\"url\":\"https:\/\/statorials.org\/nl\/\",\"name\":\"Statorials\",\"description\":\"Uw gids voor statistische competentie\",\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\/\/statorials.org\/nl\/?s={search_term_string}\"},\"query-input\":\"required name=search_term_string\"}],\"inLanguage\":\"de\"},{\"@type\":\"Person\",\"@id\":\"https:\/\/statorials.org\/nl\/#\/schema\/person\/d4b8842173cca1bb62cdec41860e4219\",\"name\":\"Dr.benjamin anderson\",\"image\":{\"@type\":\"ImageObject\",\"inLanguage\":\"de\",\"@id\":\"https:\/\/statorials.org\/nl\/#\/schema\/person\/image\/\",\"url\":\"http:\/\/statorials.org\/nl\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg\",\"contentUrl\":\"http:\/\/statorials.org\/nl\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg\",\"caption\":\"Dr.benjamin anderson\"},\"description\":\"Ik ben Benjamin, een gepensioneerde hoogleraar statistiek die nu een toegewijde Statorials-lesgever is. Ik heb uitgebreide ervaring en expertise op het gebied van statistiek en ik ben vastbesloten om mijn kennis te delen met studenten via Statorials. Lees verder\",\"sameAs\":[\"http:\/\/statorials.org\/nl\"]}]}<\/script>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Hoe een ANCOVA uit te voeren in R - Statorials","description":"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/statorials.org\/nl\/ancova-in-r\/","og_locale":"de_DE","og_type":"article","og_title":"Hoe een ANCOVA uit te voeren in R - Statorials","og_description":"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.","og_url":"https:\/\/statorials.org\/nl\/ancova-in-r\/","og_site_name":"Statorials","article_published_time":"2023-07-29T18:17:31+00:00","author":"Dr.benjamin anderson","twitter_card":"summary_large_image","twitter_misc":{"Verfasst von":"Dr.benjamin anderson","Gesch\u00e4tzte Lesezeit":"6\u00a0Minuten"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"WebPage","@id":"https:\/\/statorials.org\/nl\/ancova-in-r\/","url":"https:\/\/statorials.org\/nl\/ancova-in-r\/","name":"Hoe een ANCOVA uit te voeren in R - Statorials","isPartOf":{"@id":"https:\/\/statorials.org\/nl\/#website"},"datePublished":"2023-07-29T18:17:31+00:00","dateModified":"2023-07-29T18:17:31+00:00","author":{"@id":"https:\/\/statorials.org\/nl\/#\/schema\/person\/d4b8842173cca1bb62cdec41860e4219"},"description":"In deze tutorial wordt uitgelegd hoe u een ANCOVA (covariantieanalyse) uitvoert in R.","breadcrumb":{"@id":"https:\/\/statorials.org\/nl\/ancova-in-r\/#breadcrumb"},"inLanguage":"de","potentialAction":[{"@type":"ReadAction","target":["https:\/\/statorials.org\/nl\/ancova-in-r\/"]}]},{"@type":"BreadcrumbList","@id":"https:\/\/statorials.org\/nl\/ancova-in-r\/#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Thuis","item":"https:\/\/statorials.org\/nl\/"},{"@type":"ListItem","position":2,"name":"Hoe voer je een ancova uit in r"}]},{"@type":"WebSite","@id":"https:\/\/statorials.org\/nl\/#website","url":"https:\/\/statorials.org\/nl\/","name":"Statorials","description":"Uw gids voor statistische competentie","potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/statorials.org\/nl\/?s={search_term_string}"},"query-input":"required name=search_term_string"}],"inLanguage":"de"},{"@type":"Person","@id":"https:\/\/statorials.org\/nl\/#\/schema\/person\/d4b8842173cca1bb62cdec41860e4219","name":"Dr.benjamin anderson","image":{"@type":"ImageObject","inLanguage":"de","@id":"https:\/\/statorials.org\/nl\/#\/schema\/person\/image\/","url":"http:\/\/statorials.org\/nl\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg","contentUrl":"http:\/\/statorials.org\/nl\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg","caption":"Dr.benjamin anderson"},"description":"Ik ben Benjamin, een gepensioneerde hoogleraar statistiek die nu een toegewijde Statorials-lesgever is. Ik heb uitgebreide ervaring en expertise op het gebied van statistiek en ik ben vastbesloten om mijn kennis te delen met studenten via Statorials. Lees verder","sameAs":["http:\/\/statorials.org\/nl"]}]}},"yoast_meta":{"yoast_wpseo_title":"","yoast_wpseo_metadesc":"","yoast_wpseo_canonical":""},"_links":{"self":[{"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/posts\/484","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/comments?post=484"}],"version-history":[{"count":0,"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/posts\/484\/revisions"}],"wp:attachment":[{"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/media?parent=484"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/categories?post=484"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/statorials.org\/nl\/wp-json\/wp\/v2\/tags?post=484"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}