【问题标题】:Saved regression output of a nested regression保存的嵌套回归的回归输出
【发布时间】:2022-01-12 13:12:01
【问题描述】:

我正在尝试导出使用broommap 创建的嵌套回归的回归输出。

例如:

library(tidyverse)
library(broom)

year <- rep(2014:2015, length.out=10000)
group <- sample(c(0,1,2,3,4,5,6), replace=TRUE, size=10000)
female <- sample(c(0,1), replace=TRUE, size=10000)
smoker <- sample(c(0,1), replace=TRUE, size=10000)
dta <- tibble(year = year, group = group, female = female, smoker = smoker)

mods <- dta %>%
  nest(data = c(-year)) %>%
  mutate(model = map(data, ~ glm(smoker ~ female*group, data = .,
                                 family = binomial(link = "probit"))))

## export mods

## upload mods
  1. 如何保存模组?
  2. 保存后,如何将其加载到全局环境中?

【问题讨论】:

    标签: r tidyr broom


    【解决方案1】:

    你可以试试:

    saveRDS(mods, "/tmp/model.rds") # choose here the path to the wanted file
    
    mods <- readRDS("/tmp/model.rds")
    

    【讨论】:

    • 是的!我试过这个,但这真的很耗时。我基本上需要强制退出 Rstudios 才能摆脱它。还有其他选择吗?
    • 嗯,在我的电脑上需要 843 毫秒,不到 1 秒!所以,很难想象这么长的时间!但是,您可能需要考虑创建一个 ramdisk:直接在 RAM(内存)中保存和读取。
    • 另外,好像没有更快的办法:cran.r-project.org/web/packages/git2rdata/vignettes/…
    • 好的,谢谢!那我试试这个。
    猜你喜欢
    • 2021-05-25
    • 1970-01-01
    • 2018-08-29
    • 1970-01-01
    • 2016-02-28
    • 2020-08-18
    • 1970-01-01
    • 1970-01-01
    • 2019-10-11
    相关资源
    最近更新 更多