From b5d38a37b4a024f5bc9354d1b8a8b4d9d8c996f3 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 14:38:27 -0300 Subject: [PATCH 1/7] =?UTF-8?q?prot=C3=B3tipo=20de=20c=C3=B3digo=20para=20?= =?UTF-8?q?simplificar=20output=20simp=C3=A1tico=20dos=20modelos?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/Cap31-32/glm.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/Aulas/Cap31-32/glm.R b/Aulas/Cap31-32/glm.R index 8dc02c8..e949e4d 100644 --- a/Aulas/Cap31-32/glm.R +++ b/Aulas/Cap31-32/glm.R @@ -18,10 +18,10 @@ fisher.test(tc.idoso.osteo) glm.modelo5 <- glm(osteo ~ BMI + idade + vitD, binomial, dados.rlm) summary(glm.modelo5) -c(format.float(exp(coef(glm.modelo5)[2])), format.interval(exp(confint.default(glm.modelo5)[2, ]))) -c(format.float(exp(coef(glm.modelo5)[3])), format.interval(exp(confint.default(glm.modelo5)[3, ]))) -c(format.float(exp(coef(glm.modelo5)[4])), format.interval(exp(confint.default(glm.modelo5)[4, ]))) -c(format.float(exp(coef(glm.modelo5)[5])), format.interval(exp(confint.default(glm.modelo5)[5, ]))) +paste0("OR: ", format.float(exp(coef(glm.modelo5)[2])), ", IC: ", format.interval(exp(confint.default(glm.modelo5)[2, ]))) +paste0("OR: ", format.float(exp(coef(glm.modelo5)[3])), ", IC: ", format.interval(exp(confint.default(glm.modelo5)[3, ]))) +paste0("OR: ", format.float(exp(coef(glm.modelo5)[4]), 7), ", IC: ", format.interval(exp(confint.default(glm.modelo5)[4, ]), 7)) +paste0("OR: ", format.float(exp(coef(glm.modelo5)[5]), 7), ", IC: ", format.interval(exp(confint.default(glm.modelo5)[5, ]), 7)) # graficos ---------------------------------------------------------------- From da0996a743f7eeaa2d136e75bbce4b39afcaed00 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 14:40:39 -0300 Subject: [PATCH 2/7] =?UTF-8?q?c=C3=B3digo=20aplicado=20ao=20=C3=BAltimo?= =?UTF-8?q?=20modelo=20glm?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/Cap31-32/glm.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Aulas/Cap31-32/glm.R b/Aulas/Cap31-32/glm.R index e949e4d..f1a7f46 100644 --- a/Aulas/Cap31-32/glm.R +++ b/Aulas/Cap31-32/glm.R @@ -10,7 +10,7 @@ summary(dados.rlm) glm.modelo4 <- glm(osteo ~ idoso, binomial, dados.rlm) summary(glm.modelo4) -c(format.float(exp(coef(glm.modelo4)[2])), format.interval(exp(confint.default(glm.modelo4)[2, ]))) +paste0("OR: ", format.float(exp(coef(glm.modelo4)[2])), ", IC: ", format.interval(exp(confint.default(glm.modelo4)[2, ]))) tc.idoso.osteo <- with(dados.rlm, table(idoso, osteo)) fisher.test(tc.idoso.osteo) From aab6ff73160a733266dc71665b7b87d2cfe2a409 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 14:45:23 -0300 Subject: [PATCH 3/7] =?UTF-8?q?OR/IC=20do=20teste=20de=20Fisher=20(obs:=20?= =?UTF-8?q?diferen=C3=A7a=20num=C3=A9rica=202o=20d=C3=ADgito=20no=20IC)?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/Cap31-32/glm.R | 1 + 1 file changed, 1 insertion(+) diff --git a/Aulas/Cap31-32/glm.R b/Aulas/Cap31-32/glm.R index f1a7f46..924649b 100644 --- a/Aulas/Cap31-32/glm.R +++ b/Aulas/Cap31-32/glm.R @@ -14,6 +14,7 @@ paste0("OR: ", format.float(exp(coef(glm.modelo4)[2])), ", IC: ", format.interva tc.idoso.osteo <- with(dados.rlm, table(idoso, osteo)) fisher.test(tc.idoso.osteo) +paste0("OR: ", format.float(fisher.test(tc.idoso.osteo)$estimate), ", IC: ", format.interval(fisher.test(tc.idoso.osteo)$conf.int)) glm.modelo5 <- glm(osteo ~ BMI + idade + vitD, binomial, dados.rlm) summary(glm.modelo5) From c64baae6523fd0de3d2477e960876a2799ffc16e Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 15:03:55 -0300 Subject: [PATCH 4/7] =?UTF-8?q?N=C3=A3o=20carregar=20o=20script=20RLM=20no?= =?UTF-8?q?=20GLM=20(ggplots=20podem=20ser=20lentos)?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/Cap31-32/glm.R | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/Aulas/Cap31-32/glm.R b/Aulas/Cap31-32/glm.R index 924649b..b5f0efe 100644 --- a/Aulas/Cap31-32/glm.R +++ b/Aulas/Cap31-32/glm.R @@ -2,7 +2,9 @@ library(philsfmisc) # dados simulados --------------------------------------------------------- -source('~/Documents/Docencia/Bioestatistica/Aulas/Cap31-32/rlm.R') +dados.rlm <- fread("Aulas/Cap31-32/dados-rlm.csv", stringsAsFactors = TRUE) +dados.rlm$vitD <- factor(dados.rlm$vitD, levels = c("baixa", "media", "alta")) +dados.rlm$osteo <- relevel(dados.rlm$osteo, "Sadio") summary(dados.rlm) # modelos ----------------------------------------------------------------- From d93db291879ef486d9c70fa5f0b9fae454fa0336 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 15:54:05 -0300 Subject: [PATCH 5/7] =?UTF-8?q?c=C3=B3digo=20de=20output=20para=20modelos?= =?UTF-8?q?=20RLM?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/Cap31-32/rlm.R | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/Aulas/Cap31-32/rlm.R b/Aulas/Cap31-32/rlm.R index bd8c12c..7281b53 100644 --- a/Aulas/Cap31-32/rlm.R +++ b/Aulas/Cap31-32/rlm.R @@ -17,10 +17,14 @@ rlm.modelo2.2 <- lm(BMD ~ BMI + vitD, data = dados.rlm) rlm.modelo3 <- lm(BMD ~ BMI + idade + vitD, data = dados.rlm) print(summary(rlm.modelo1)) +paste0(format.float(coef(rlm.modelo1)[2]), ", IC: ", format.interval(confint.default(rlm.modelo1)[2, ])) print(summary(rlm.modelo2)) print(summary(rlm.modelo2.1)) +paste0(format.float(coef(rlm.modelo2.1)[2]), ", IC: ", format.interval(confint.default(rlm.modelo2.1)[2, ])) print(summary(rlm.modelo2.2)) +paste0(format.float(coef(rlm.modelo2.2)[2]), ", IC: ", format.interval(confint.default(rlm.modelo2.2)[2, ])) print(summary(rlm.modelo3)) +paste0(format.float(coef(rlm.modelo3)[2]), ", IC: ", format.interval(confint.default(rlm.modelo3)[2, ])) # graficos ---------------------------------------------------------------- From 3eb48a4445f6504844d5f27bf345b96da6144055 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 4 Nov 2019 15:54:15 -0300 Subject: [PATCH 6/7] =?UTF-8?q?frame=20de=20conclus=C3=A3o=20da=20se=C3=A7?= =?UTF-8?q?=C3=A3o=20RLM=20-=20todos=20os=20modelos=20juntos?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Aulas/BE-Regressao-Logistica-cap31-32.tex | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) diff --git a/Aulas/BE-Regressao-Logistica-cap31-32.tex b/Aulas/BE-Regressao-Logistica-cap31-32.tex index ebe4f57..21a8f1e 100644 --- a/Aulas/BE-Regressao-Logistica-cap31-32.tex +++ b/Aulas/BE-Regressao-Logistica-cap31-32.tex @@ -695,6 +695,29 @@ \subsection{Regressão Linear Múltipla} \end{exampleblock} \end{frame} +\begin{frame}{\scriptsize Conclusão} + \begin{center} + Qual das quatro estimativas de BMI é a melhor? + + \bigskip + Por que? + \end{center} +\end{frame} + +\begin{frame}{\scriptsize Qual das estimativas é melhor? Por que?} + \begin{exampleblock}{\scriptsize Modelo 1 -- BMI (estimativa ``crua'')} + -1.99, IC: [-2.40, -1.59] + \end{exampleblock} + \begin{exampleblock}{\scriptsize Modelo 2.1 -- BMI ajustado por idade} + 2.04, IC: [-2.37, -1.70] + \end{exampleblock} + \begin{exampleblock}{\scriptsize Modelo 2.2 -- BMI ajustado por vitamina D sérica} + -1.91, IC: [-2.14, -1.68] + \end{exampleblock} + \begin{exampleblock}{\scriptsize Modelo 3 -- BMI ajustado por idade + vitamina D sérica} + -1.95, IC: [-2.05, -1.85] + \end{exampleblock} +\end{frame} \section{Regressão Logística} \subsection{Regressão Logística} From 7e24aa0fc36e4abdf7f9d2dfcd5cea07cae797e4 Mon Sep 17 00:00:00 2001 From: Felipe Figueiredo Date: Mon, 18 Nov 2019 17:47:23 -0300 Subject: [PATCH 7/7] typo (sinal negativo) --- Aulas/BE-Regressao-Logistica-cap31-32.tex | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/Aulas/BE-Regressao-Logistica-cap31-32.tex b/Aulas/BE-Regressao-Logistica-cap31-32.tex index 21a8f1e..d7f8d86 100644 --- a/Aulas/BE-Regressao-Logistica-cap31-32.tex +++ b/Aulas/BE-Regressao-Logistica-cap31-32.tex @@ -709,7 +709,7 @@ \subsection{Regressão Linear Múltipla} -1.99, IC: [-2.40, -1.59] \end{exampleblock} \begin{exampleblock}{\scriptsize Modelo 2.1 -- BMI ajustado por idade} - 2.04, IC: [-2.37, -1.70] + -2.04, IC: [-2.37, -1.70] \end{exampleblock} \begin{exampleblock}{\scriptsize Modelo 2.2 -- BMI ajustado por vitamina D sérica} -1.91, IC: [-2.14, -1.68] @@ -718,6 +718,7 @@ \subsection{Regressão Linear Múltipla} -1.95, IC: [-2.05, -1.85] \end{exampleblock} \end{frame} + \section{Regressão Logística} \subsection{Regressão Logística}