使用 modelsummary 報告 mgcv::gam 摘要

Question

我正在嘗試使用模型摘要 package 從mgcv::gam()報告 model 摘要。 靈活表 package 提供的摘要與 R 中的摘要 output 以及出版物中經常出現的內容一致。 它將固定/參數效應和平滑項的報告分開。 雖然 flextable 效果很好，但我還是想用 modelsummary（主要是因為它能夠 output 到 gt、kable 等）。 我的計划是生成兩個單獨的表格，並分別報告參數和平滑項的適當數據（可能有更好的方法？）。 但是，我在嘗試省略modelsummary()中的系數時掛斷了電話。

彈性表示例：

library(mgcv)
library(flextable)
library(modelsummary)

dat <- gamSim(1, n = 4000, dist = "normal", scale = 2)
mod <- gam(y ~ s(x0) + s(x1) + s(x2), data = dat)
flextable::as_flextable(mod)

我使用modelsummary()獲取參數項摘要的第一步：

modelsummary(mod,
             estimate = "estimate",
             statistic = c("Std.Error" = "std.error", 
                           "t-value" = "statistic", 
                           "p-value" = "p.value"),
             shape = term ~ model + statistic,
             gof_map = NA)

我想刪除平滑術語並將它們包含在不同的表或組中，所以我嘗試了 coef_omit 參數：

modelsummary(mod,
             estimate = "estimate",
             statistic = c("Std.Error" = "std.error", 
                           "t-value" = "statistic", 
                           "p-value" = "p.value"),
             coef_omit = "^(?!.*Intercept)", #this should retain the intercept term
             omit = ".*",
             shape = term ~ model + statistic,
             gof_map = NA)

Error in if (dat$part[i] == "estimates" && dat[[column]][i - 1] == dat[[column]][i]) { : 
  missing value where TRUE/FALSE needed

有趣的是，如果我刪除 shape 參數以“長格式”報告統計信息，錯誤就會消失。 我可能正在接近完全錯誤地格式化此摘要，並且願意接受建議。