我對回圈三個邏輯回歸分析的結果有疑問。我正在嘗試將 exp(b)/odds 比率及其 95% 置信區間并排放入資料框中。我想在列中包含模型名稱（例如模型 1、模型 2 和模型 3），在行中包含變數。我一直在嘗試修改從這個站點找到的代碼——（我想在發布我的問題之前我會在這里檢查一下，但沒有發現任何問題，但如果我錯過了，請指出我的方向）。不幸的是，這似乎有點超出我的知識基礎。這是我正在嘗試使用的代碼。它進行基本計算，并僅使用系數創建一個小標題，但我需要每個模型的 exp(b)/odds 比和 CI，正如我所說的，列中的模型名稱和行中的變數。
library(tidyverse)
mydata <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv")
forms <- list(model1 = admit ~ gre, model2 = admit ~ gre gpa, model3 = admit ~ gre gpa rank)
map_df(forms, ~coef(glm(.x, mydata, family = "binomial"))) %>%
select(-1)
我將如何修改腳本以收集我需要的資訊，以及我的研究需要的資料框樣式？
uj5u.com熱心網友回復：
broom::tidy()對于這些問題非常方便。嘗試
m1 <- (forms
|> map(glm, data = mydata, family = "binomial")
|> map_dfr(broom::tidy, exponentiate = TRUE, conf.int = TRUE, .id = "model")
)
print(m1)
select()并且pivot_wider()對于轉換很有用，例如
(m1
|> dplyr::select(model, term, estimate, conf.low, conf.high)
|> pivot_wider(names_from = "term",
values_from = c("estimate", "conf.low", "conf.high"))
)
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/518289.html
標籤：r循环字典逻辑回归