{"id":2948,"date":"2023-07-19T23:07:53","date_gmt":"2023-07-19T23:07:53","guid":{"rendered":"https:\/\/statorials.org\/pt\/glm-r-carre\/"},"modified":"2023-07-19T23:07:53","modified_gmt":"2023-07-19T23:07:53","slug":"glm-r-carre","status":"publish","type":"post","link":"https:\/\/statorials.org\/pt\/glm-r-carre\/","title":{"rendered":"Como calcular r-quadrado para glm em r"},"content":{"rendered":"<p><\/p>\n<hr>\n<p><span style=\"color: #000000;\">Freq\u00fcentemente, quando ajustamos um modelo de regress\u00e3o linear, usamos <strong>R-quadrado<\/strong> para avaliar qu\u00e3o bem um modelo se ajusta aos dados.<\/span><\/p>\n<p> <span style=\"color: #000000;\">R ao quadrado representa a propor\u00e7\u00e3o da vari\u00e2ncia na <a href=\"https:\/\/statorials.org\/pt\/respostas-explicativas-das-variaveis\/\" target=\"_blank\" rel=\"noopener\">vari\u00e1vel resposta<\/a> que pode ser explicada pelas vari\u00e1veis preditoras em um modelo de regress\u00e3o.<\/span><\/p>\n<p> <span style=\"color: #000000;\">Esse n\u00famero varia de 0 a 1, com valores mais altos indicando melhor ajuste do modelo.<\/span><\/p>\n<p> <span style=\"color: #000000;\">No entanto, n\u00e3o h\u00e1 valor de R ao quadrado para modelos lineares gerais, como modelos <a href=\"https:\/\/statorials.org\/pt\/regressao-logistica-1\/\" target=\"_blank\" rel=\"noopener\">de regress\u00e3o log\u00edstica<\/a> e modelos <a href=\"https:\/\/statorials.org\/pt\/regressao-de-peixe\/\" target=\"_blank\" rel=\"noopener\">de regress\u00e3o de Poisson<\/a> .<\/span><\/p>\n<p> <span style=\"color: #000000;\">Em vez disso, podemos calcular uma m\u00e9trica conhecida como <strong>R-Squared de McFadden<\/strong> , que varia de 0 a pouco menos de 1, com valores mais altos indicando melhor ajuste do modelo.<\/span><\/p>\n<p> <span style=\"color: #000000;\">Usamos a seguinte f\u00f3rmula para calcular o R ao quadrado de McFadden:<\/span><\/p>\n<p> <span style=\"color: #000000;\">R-quadrado de McFadden = 1 &#8211; ( <sub>modelo<\/sub> de log de verossimilhan\u00e7a \/ log de probabilidade <sub>zero<\/sub> )<\/span><\/p>\n<p> <span style=\"color: #000000;\">Ouro:<\/span><\/p>\n<ul>\n<li> <span style=\"color: #000000;\"><strong><sub>modelo<\/sub> de log de verossimilhan\u00e7a<\/strong> : valor de log de verossimilhan\u00e7a do modelo ajustado atual<\/span><\/li>\n<li> <span style=\"color: #000000;\"><strong>log de probabilidade <sub>zero<\/sub><\/strong> : valor de log de probabilidade do modelo nulo (modelo com intercepta\u00e7\u00e3o apenas)<\/span><\/li>\n<\/ul>\n<p> <span style=\"color: #000000;\">Na pr\u00e1tica, valores acima de 0,40 indicam que um modelo se ajusta muito bem aos dados.<\/span><\/p>\n<p> <span style=\"color: #000000;\">O exemplo a seguir mostra como calcular o R-quadrado de McFadden para um modelo de regress\u00e3o log\u00edstica em R.<\/span><\/p>\n<h3> <span style=\"color: #000000;\"><strong>Exemplo: Calculando R-quadrado de McFadden em R<\/strong><\/span><\/h3>\n<p> <span style=\"color: #000000;\">Para este exemplo, usaremos o conjunto de dados <strong>padr\u00e3o<\/strong> do pacote ISLR. Podemos usar o seguinte c\u00f3digo para carregar e exibir um resumo do conjunto de dados:<\/span><\/p>\n<pre style=\"background-color: #ececec; font-size: 15px;\"> <strong><span style=\"color: #008080;\">#install and load ISLR package\n<span style=\"color: #000000;\">install. <span style=\"color: #3366ff;\">packages<\/span> (' <span style=\"color: #ff0000;\">ISLR<\/span> ')<\/span>\n<span style=\"color: #000000;\"><span style=\"color: #008000;\">library<\/span> (ISLR)<\/span>\n\n#define dataset<\/span>\ndata &lt;- ISLR::Default\n\n<span style=\"color: #008080;\">#view summary of dataset\n<\/span>summary(data)\n\n default student balance income     \n No:9667 No:7056 Min. : 0.0 Min. : 772  \n Yes: 333 Yes:2944 1st Qu.: 481.7 1st Qu.:21340  \n                       Median: 823.6 Median: 34553  \n                       Mean: 835.4 Mean: 33517  \n                       3rd Qu.:1166.3 3rd Qu.:43808  \n                       Max. :2654.3 Max. :73554  \n\n<span style=\"color: #008080;\">#find total observations in dataset<\/span>\nnrow(data)\n\n[1] 10000\n<\/strong><\/pre>\n<p> <span style=\"color: #000000;\">Este conjunto de dados cont\u00e9m as seguintes informa\u00e7\u00f5es sobre 10.000 indiv\u00edduos:<\/span><\/p>\n<ul>\n<li> <span style=\"color: #000000;\"><strong>inadimpl\u00eancia:<\/strong> indica se um indiv\u00edduo est\u00e1 inadimplente ou n\u00e3o.<\/span><\/li>\n<li> <span style=\"color: #000000;\"><strong>estudante:<\/strong> indica se um indiv\u00edduo \u00e9 estudante ou n\u00e3o.<\/span><\/li>\n<li> <span style=\"color: #000000;\"><strong>saldo:<\/strong> Saldo m\u00e9dio mantido por um indiv\u00edduo.<\/span><\/li>\n<li> <span style=\"color: #000000;\"><strong>renda:<\/strong> Renda da pessoa f\u00edsica.<\/span><\/li>\n<\/ul>\n<p> <span style=\"color: #000000;\">Usaremos a situa\u00e7\u00e3o de estudante, o saldo banc\u00e1rio e a renda para construir um modelo de regress\u00e3o log\u00edstica que prev\u00ea a probabilidade de um determinado indiv\u00edduo entrar em default:<\/span><\/p>\n<pre style=\"background-color: #ececec; font-size: 15px;\"> <strong><span style=\"color: #008080;\">#fit logistic regression model<\/span>\nmodel &lt;- glm(default~student+balance+income, family=' <span style=\"color: #ff0000;\">binomial<\/span> ', data=data)\n\n<span style=\"color: #008080;\">#view model summary<\/span>\nsummary(model)\n\nCall:\nglm(formula = default ~ balance + student + income, family = \"binomial\", \n    data = data)\n\nDeviance Residuals: \n    Min 1Q Median 3Q Max  \n-2.4691 -0.1418 -0.0557 -0.0203 3.7383  \n\nCoefficients:\n              Estimate Std. Error z value Pr(&gt;|z|)    \n(Intercept) -1.087e+01 4.923e-01 -22.080 &lt; 2e-16 ***\nbalance 5.737e-03 2.319e-04 24.738 &lt; 2e-16 ***\nstudentYes -6.468e-01 2.363e-01 -2.738 0.00619 ** \nincome 3.033e-06 8.203e-06 0.370 0.71152    \n---\nSignificant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n\n(Dispersion parameter for binomial family taken to be 1)\n\n    Null deviance: 2920.6 on 9999 degrees of freedom\nResidual deviance: 1571.5 on 9996 degrees of freedom\nAIC: 1579.5\n\nNumber of Fisher Scoring iterations: 8\n<\/strong><\/pre>\n<p> <span style=\"color: #000000;\"><span style=\"color: #000000;\">A seguir, usaremos a seguinte f\u00f3rmula para calcular o valor R ao quadrado de McFadden para este modelo:<\/span><\/span><\/p>\n<pre style=\"background-color: #ececec; font-size: 15px;\"> <strong><span style=\"color: #008080;\">#calculate McFadden's R-squared for model<\/span>\nwith(summary(model), 1 - deviance\/null. <span style=\"color: #3366ff;\">deviance<\/span> )\n\n[1] 0.4619194\n<\/strong><\/pre>\n<p> <span style=\"color: #000000;\">O valor R ao quadrado de McFadden \u00e9 <strong>0,4619194<\/strong> . Este valor \u00e9 bastante elevado, indicando que nosso modelo se ajusta bem aos dados e possui alto poder preditivo.<\/span><\/p>\n<p> <span style=\"color: #000000;\">Observe tamb\u00e9m que tamb\u00e9m poder\u00edamos usar a fun\u00e7\u00e3o <strong>pR2()<\/strong> do pacote <strong>pscl<\/strong> para calcular o valor R-quadrado de McFadden para o modelo:<\/span><\/p>\n<pre style=\"background-color: #ececec; font-size: 15px;\"> <strong><span style=\"color: #008080;\">#install and load pscl package<\/span>\ninstall. <span style=\"color: #3366ff;\">packages<\/span> (' <span style=\"color: #ff0000;\">pscl<\/span> ')\n<span style=\"color: #008000;\">library<\/span> (pscl)\n\n<span style=\"color: #008080;\">#calculate McFadden's R-squared for model<\/span>\npR2(model)[' <span style=\"color: #ff0000;\">McFadden<\/span> ']\n\n McFadden \n0.4619194\n<\/strong><\/pre>\n<p> <span style=\"color: #000000;\"><span style=\"color: #000000;\">Observe que este valor corresponde ao calculado anteriormente.<\/span><\/span><\/p>\n<h3> <span style=\"color: #000000;\"><strong>Recursos adicionais<\/strong><\/span><\/h3>\n<p> <span style=\"color: #000000;\">Os tutoriais a seguir explicam como realizar outras tarefas comuns em R:<\/span><\/p>\n<p> <a href=\"https:\/\/statorials.org\/pt\/r-elevado-ao-quadrado-de-r\/\" target=\"_blank\" rel=\"noopener\">Como calcular R-quadrado em R<\/a><br \/> <a href=\"https:\/\/statorials.org\/pt\/r-quadrados-em-r-se-ajusta\/\" target=\"_blank\" rel=\"noopener\">Como calcular R-quadrado ajustado em R<\/a><br \/> <a href=\"https:\/\/statorials.org\/pt\/bom-valor-de-r-ao-quadrado\/\" target=\"_blank\" rel=\"noopener\">O que \u00e9 um bom valor de R ao quadrado?<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Freq\u00fcentemente, quando ajustamos um modelo de regress\u00e3o linear, usamos R-quadrado para avaliar qu\u00e3o bem um modelo se ajusta aos dados. R ao quadrado representa a propor\u00e7\u00e3o da vari\u00e2ncia na vari\u00e1vel resposta que pode ser explicada pelas vari\u00e1veis preditoras em um modelo de regress\u00e3o. Esse n\u00famero varia de 0 a 1, com valores mais altos indicando [&hellip;]<\/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-2948","post","type-post","status-publish","format-standard","hentry","category-guia"],"yoast_head":"<!-- This site is optimized with the Yoast SEO plugin v21.5 - https:\/\/yoast.com\/wordpress\/plugins\/seo\/ -->\n<title>Como calcular R-quadrado para glm em R - Estatologia<\/title>\n<meta name=\"description\" content=\"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.\" \/>\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\/pt\/glm-r-carre\/\" \/>\n<meta property=\"og:locale\" content=\"pt_PT\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Como calcular R-quadrado para glm em R - Estatologia\" \/>\n<meta property=\"og:description\" content=\"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.\" \/>\n<meta property=\"og:url\" content=\"https:\/\/statorials.org\/pt\/glm-r-carre\/\" \/>\n<meta property=\"og:site_name\" content=\"Statorials\" \/>\n<meta property=\"article:published_time\" content=\"2023-07-19T23:07:53+00:00\" \/>\n<meta name=\"author\" content=\"Dr. benjamim anderson\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Escrito por\" \/>\n\t<meta name=\"twitter:data1\" content=\"Dr. benjamim anderson\" \/>\n\t<meta name=\"twitter:label2\" content=\"Tempo estimado de leitura\" \/>\n\t<meta name=\"twitter:data2\" content=\"3 minutos\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\/\/schema.org\",\"@graph\":[{\"@type\":\"WebPage\",\"@id\":\"https:\/\/statorials.org\/pt\/glm-r-carre\/\",\"url\":\"https:\/\/statorials.org\/pt\/glm-r-carre\/\",\"name\":\"Como calcular R-quadrado para glm em R - Estatologia\",\"isPartOf\":{\"@id\":\"https:\/\/statorials.org\/pt\/#website\"},\"datePublished\":\"2023-07-19T23:07:53+00:00\",\"dateModified\":\"2023-07-19T23:07:53+00:00\",\"author\":{\"@id\":\"https:\/\/statorials.org\/pt\/#\/schema\/person\/e08f98e8db95e0aa9c310e1b27c9c666\"},\"description\":\"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.\",\"breadcrumb\":{\"@id\":\"https:\/\/statorials.org\/pt\/glm-r-carre\/#breadcrumb\"},\"inLanguage\":\"pt-PT\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\/\/statorials.org\/pt\/glm-r-carre\/\"]}]},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\/\/statorials.org\/pt\/glm-r-carre\/#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Lar\",\"item\":\"https:\/\/statorials.org\/pt\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"Como calcular r-quadrado para glm em r\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\/\/statorials.org\/pt\/#website\",\"url\":\"https:\/\/statorials.org\/pt\/\",\"name\":\"Statorials\",\"description\":\"O seu guia para a literacia estat\u00edstica!\",\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\/\/statorials.org\/pt\/?s={search_term_string}\"},\"query-input\":\"required name=search_term_string\"}],\"inLanguage\":\"pt-PT\"},{\"@type\":\"Person\",\"@id\":\"https:\/\/statorials.org\/pt\/#\/schema\/person\/e08f98e8db95e0aa9c310e1b27c9c666\",\"name\":\"Dr. benjamim anderson\",\"image\":{\"@type\":\"ImageObject\",\"inLanguage\":\"pt-PT\",\"@id\":\"https:\/\/statorials.org\/pt\/#\/schema\/person\/image\/\",\"url\":\"https:\/\/statorials.org\/pt\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg\",\"contentUrl\":\"https:\/\/statorials.org\/pt\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg\",\"caption\":\"Dr. benjamim anderson\"},\"description\":\"Ol\u00e1, sou Benjamin, um professor aposentado de estat\u00edstica que se tornou professor dedicado na Statorials. Com vasta experi\u00eancia e conhecimento na \u00e1rea de estat\u00edstica, estou empenhado em compartilhar meu conhecimento para capacitar os alunos por meio de Statorials. Saber mais\",\"sameAs\":[\"https:\/\/statorials.org\/pt\"]}]}<\/script>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Como calcular R-quadrado para glm em R - Estatologia","description":"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.","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\/pt\/glm-r-carre\/","og_locale":"pt_PT","og_type":"article","og_title":"Como calcular R-quadrado para glm em R - Estatologia","og_description":"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.","og_url":"https:\/\/statorials.org\/pt\/glm-r-carre\/","og_site_name":"Statorials","article_published_time":"2023-07-19T23:07:53+00:00","author":"Dr. benjamim anderson","twitter_card":"summary_large_image","twitter_misc":{"Escrito por":"Dr. benjamim anderson","Tempo estimado de leitura":"3 minutos"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"WebPage","@id":"https:\/\/statorials.org\/pt\/glm-r-carre\/","url":"https:\/\/statorials.org\/pt\/glm-r-carre\/","name":"Como calcular R-quadrado para glm em R - Estatologia","isPartOf":{"@id":"https:\/\/statorials.org\/pt\/#website"},"datePublished":"2023-07-19T23:07:53+00:00","dateModified":"2023-07-19T23:07:53+00:00","author":{"@id":"https:\/\/statorials.org\/pt\/#\/schema\/person\/e08f98e8db95e0aa9c310e1b27c9c666"},"description":"Este tutorial explica como calcular um valor pseudo R-quadrado para modelos glm em R, incluindo um exemplo completo.","breadcrumb":{"@id":"https:\/\/statorials.org\/pt\/glm-r-carre\/#breadcrumb"},"inLanguage":"pt-PT","potentialAction":[{"@type":"ReadAction","target":["https:\/\/statorials.org\/pt\/glm-r-carre\/"]}]},{"@type":"BreadcrumbList","@id":"https:\/\/statorials.org\/pt\/glm-r-carre\/#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Lar","item":"https:\/\/statorials.org\/pt\/"},{"@type":"ListItem","position":2,"name":"Como calcular r-quadrado para glm em r"}]},{"@type":"WebSite","@id":"https:\/\/statorials.org\/pt\/#website","url":"https:\/\/statorials.org\/pt\/","name":"Statorials","description":"O seu guia para a literacia estat\u00edstica!","potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/statorials.org\/pt\/?s={search_term_string}"},"query-input":"required name=search_term_string"}],"inLanguage":"pt-PT"},{"@type":"Person","@id":"https:\/\/statorials.org\/pt\/#\/schema\/person\/e08f98e8db95e0aa9c310e1b27c9c666","name":"Dr. benjamim anderson","image":{"@type":"ImageObject","inLanguage":"pt-PT","@id":"https:\/\/statorials.org\/pt\/#\/schema\/person\/image\/","url":"https:\/\/statorials.org\/pt\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg","contentUrl":"https:\/\/statorials.org\/pt\/wp-content\/uploads\/2023\/10\/Dr.-Benjamin-Anderson-96x96.jpg","caption":"Dr. benjamim anderson"},"description":"Ol\u00e1, sou Benjamin, um professor aposentado de estat\u00edstica que se tornou professor dedicado na Statorials. Com vasta experi\u00eancia e conhecimento na \u00e1rea de estat\u00edstica, estou empenhado em compartilhar meu conhecimento para capacitar os alunos por meio de Statorials. Saber mais","sameAs":["https:\/\/statorials.org\/pt"]}]}},"yoast_meta":{"yoast_wpseo_title":"","yoast_wpseo_metadesc":"","yoast_wpseo_canonical":""},"_links":{"self":[{"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/posts\/2948","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/comments?post=2948"}],"version-history":[{"count":0,"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/posts\/2948\/revisions"}],"wp:attachment":[{"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/media?parent=2948"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/categories?post=2948"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/statorials.org\/pt\/wp-json\/wp\/v2\/tags?post=2948"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}