{"id":1176,"date":"2023-07-27T09:44:39","date_gmt":"2023-07-27T09:44:39","guid":{"rendered":"https:\/\/statorials.org\/pt\/validacao-cruzada-k-fold-em-r\/"},"modified":"2023-07-27T09:44:39","modified_gmt":"2023-07-27T09:44:39","slug":"validacao-cruzada-k-fold-em-r","status":"publish","type":"post","link":"https:\/\/statorials.org\/pt\/validacao-cruzada-k-fold-em-r\/","title":{"rendered":"Valida\u00e7\u00e3o cruzada k-fold em r (passo a passo)"},"content":{"rendered":"

<\/p>\n

\n

Para avaliar o desempenho de um modelo em um conjunto de dados, precisamos medir at\u00e9 que ponto as previs\u00f5es feitas pelo modelo correspondem aos dados observados.<\/span><\/p>\n

Um m\u00e9todo comumente usado para fazer isso \u00e9 conhecido como valida\u00e7\u00e3o cruzada k-fold<\/a> , que usa a seguinte abordagem:<\/span><\/p>\n

1.<\/strong> Divida aleatoriamente um conjunto de dados em k<\/em> grupos, ou \u201cdobras\u201d, de tamanho aproximadamente igual.<\/span><\/p>\n

2.<\/strong> Escolha uma das dobras como conjunto de conten\u00e7\u00e3o. Ajuste o modelo \u00e0s dobras k-1 restantes. Calcule o teste MSE nas observa\u00e7\u00f5es da camada que foi tensionada.<\/span><\/p>\n

3.<\/strong> Repita esse processo k<\/em> vezes, cada vez usando um conjunto diferente como conjunto de exclus\u00e3o.<\/span><\/p>\n

4.<\/strong> Calcule o MSE geral do teste como a m\u00e9dia dos k<\/em> MSEs do teste.<\/span><\/p>\n

A maneira mais f\u00e1cil de realizar a valida\u00e7\u00e3o cruzada k-fold em R \u00e9 usar a fun\u00e7\u00e3o trainControl()<\/a> da biblioteca caret<\/strong> em R.<\/span><\/p>\n

Este tutorial fornece um exemplo r\u00e1pido de como usar esta fun\u00e7\u00e3o para realizar valida\u00e7\u00e3o cruzada k-fold para um determinado modelo em R.<\/span><\/p>\n

Exemplo: valida\u00e7\u00e3o cruzada K-Fold em R<\/strong><\/span><\/h3>\n

Suponha que temos o seguinte conjunto de dados em R:<\/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
 O c\u00f3digo a seguir mostra como ajustar um modelo de regress\u00e3o linear m\u00faltipla<\/a> a este conjunto de dados em R e realizar valida\u00e7\u00e3o cruzada k-fold com k = 5 vezes para avaliar o desempenho do modelo:<\/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
 Veja como interpretar o resultado:<\/span><\/p>\n