![](/img/trans.png)
[英]How to include confidence intervals from multiple models in tidy output using broom?
[英]Exponentiating GAM coefficients and confidence intervals using broom::tidy
我正在使用mgcv::gam
運行廣義加法模型,並嘗試使用broom:tidy
組織我的結果,但tidy
顯然不會對 GAM 的系數或置信區間進行指數化,盡管它適用於常規glm
模型。 是否有broom::tidy
方法對 GAM 中的系數和 CI 求冪? 我特別詢問tidy
因為我想在gtsummary
創建的回歸表中使用結果。
library(tidyverse)
library(magrittr)
library(mgcv)
library(parameters)
library(gtsummary)
library(broom)
# sample data
id <- 1:2000
gender <- sample(0:1, 2000, replace = T)
age <- sample(17:64, 2000, replace = T)
race <- sample(0:1, 2000, replace = T)
health_score <- sample(0:25, 2000, replace = T)
dead <- sample(0:1, 2000, replace = T)
days_enrolled <- sample(30:3000, 2000, replace = T)
df <- data.frame(id, gender, age, race, health_score, dead, days_enrolled)
# model
model <- gam(dead ~ gender + s(age) + race + s(health_score) + offset(log(days_enrolled)),
data = df, method = "REML", family = nb())
# both give the same output:
tidy(model, parametric = T, conf.int = T)
tidy(model, parametric = T, conf.int = T, exponentiate = T)
您可以直接使用tbl_regression()
對結果求冪。 如果這不是你所追求的,請告訴我。
library(tidyverse)
library(mgcv)
library(parameters)
library(gtsummary)
library(broom)
# sample data
id <- 1:2000
gender <- sample(0:1, 2000, replace = T)
age <- sample(17:64, 2000, replace = T)
race <- sample(0:1, 2000, replace = T)
health_score <- sample(0:25, 2000, replace = T)
dead <- sample(0:1, 2000, replace = T)
days_enrolled <- sample(30:3000, 2000, replace = T)
df <- data.frame(id, gender, age, race, health_score, dead, days_enrolled)
# model
model <- gam(dead ~ gender + s(age) + race + s(health_score) + offset(log(days_enrolled)),
data = df, method = "REML", family = nb())
tbl_regression(model, exponentiate = TRUE)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.