{"id":1177,"date":"2023-07-27T09:44:39","date_gmt":"2023-07-27T09:44:39","guid":{"rendered":"https:\/\/statorials.org\/nl\/k-voudige-kruisvalidatie-in-r\/"},"modified":"2023-07-27T09:44:39","modified_gmt":"2023-07-27T09:44:39","slug":"k-voudige-kruisvalidatie-in-r","status":"publish","type":"post","link":"https:\/\/statorials.org\/nl\/k-voudige-kruisvalidatie-in-r\/","title":{"rendered":"K-voudige kruisvalidatie in r (stap voor stap)"},"content":{"rendered":"

<\/p>\n

\n

Om de prestaties van een model op een dataset te evalueren, moeten we meten hoe goed de voorspellingen van het model overeenkomen met de waargenomen gegevens.<\/span><\/p>\n

Een veelgebruikte methode om dit te doen staat bekend als k-fold cross-validatie<\/a> , waarbij de volgende aanpak wordt gebruikt:<\/span><\/p>\n

1.<\/strong> Verdeel een dataset willekeurig in k<\/em> groepen, of \u2018vouwen\u2019, van ongeveer gelijke grootte.<\/span><\/p>\n

2.<\/strong> Kies een van de vouwen als bevestigingsset. Pas de sjabloon aan de resterende k-1-vouwen aan. Bereken de MSE-proef op de waarnemingen in de gespannen lamel.<\/span><\/p>\n

3.<\/strong> Herhaal dit proces k<\/em> keer, telkens met een andere set als uitsluitingsset.<\/span><\/p>\n

4.<\/strong> Bereken de totale test-MSE als het gemiddelde van de k-<\/em> test-MSE’s.<\/span><\/p>\n

De eenvoudigste manier om k-voudige kruisvalidatie uit te voeren in R is door de functie trainControl()<\/a> uit de caret-<\/strong> bibliotheek in R te gebruiken.<\/span><\/p>\n

Deze tutorial geeft een snel voorbeeld van hoe u deze functie kunt gebruiken om k-voudige kruisvalidatie uit te voeren voor een bepaald model in R.<\/span><\/p>\n

Voorbeeld: K-voudige kruisvalidatie in R<\/strong><\/span><\/h3>\n

Stel dat we de volgende dataset in R hebben:<\/span><\/p>\n

 #create data frame<\/span>\ndf <- data.frame(y=c(6, 8, 12, 14, 14, 15, 17, 22, 24, 23),\n                 x1=c(2, 5, 4, 3, 4, 6, 7, 5, 8, 9),\n                 x2=c(14, 12, 12, 13, 7, 8, 7, 4, 6, 5))\n\n#view data frame\n<\/span>df\n\ny x1 x2\n6 2 14\n8 5 12\n12 4 12\n14 3 13\n14 4 7\n15 6 8\n17 7 7\n22 5 4\n24 8 6\n23 9 5\n<\/strong><\/pre>\n
 De volgende code laat zien hoe u een meervoudig lineair regressiemodel<\/a> aan deze gegevensset in R kunt aanpassen en k-voudige kruisvalidatie met k = 5 keer kunt uitvoeren om de prestaties van het model te evalueren:<\/span><\/p>\n
 library<\/span> (caret)<\/span>\n\n#specify the cross-validation method<\/span>\nctrl <- trainControl(method = \" cv<\/span> \", number = 5)\n\n#fit a regression model and use k-fold CV to evaluate performance\n<\/span>model <- train(y ~ x1 + x2, data = df, method = \" lm<\/span> \", trControl = ctrl)\n\n#view summary of k-fold CV               \n<\/span>print(model)\n\nLinear Regression \n\n10 samples\n 2 predictors\n\nNo pre-processing\nResampling: Cross-Validated (5 fold) \nSummary of sample sizes: 8, 8, 8, 8, 8 \nResampling results:\n\n  RMSE Rsquared MAE     \n  3.018979 1 2.882348\n\nTuning parameter 'intercept' was held constant at a value of TRUE\n<\/strong><\/pre>\n
 Zo interpreteert u het resultaat:<\/span><\/p>\n