在阅读了 tidyverse 的更多内容后,我开始立即拟合许多线性模型,如中所述this https://blog.rstudio.org/2016/02/02/tidyr-0-4-0/。也就是说,我会按照以下方式做一些事情:
library(dplyr)
library(tidyr)
library(purrr)
df <- data.frame(y = rnorm(10),
x1 = runif(10),
x2 = runif(10))
df %>%
gather(covariate, value, x1:x2) %>%
group_by(covariate) %>%
nest() %>%
mutate(model = map(.x = data , .f = ~lm(y ~ value, data = .))) %>%
mutate(rsquared = map_dbl(.x = model, .f = ~summary(.)$r.squared))
问题是,当变量类型不同时,例如当一个是数字而一个是因子时,这种方法会失败,因为gather()
函数将强制整体value
向量化为因子。例如,
df <- data.frame(y = rnorm(10),
x1 = runif(10),
x3 = sample(c("a", "b", "c"), 10, replace = TRUE))
df %>%
gather(covariate, value, x1:x3) %>%
sapply(class)
随后是警告
Warning message:
attributes are not identical across measure variables; they will be dropped
y covariate value
"numeric" "character" "character"
and the value
列是一个字符,所以技巧nest()
将不再起作用,因为所有协变量都将作为因子放入。
我想知道是否有tidy这样做的方法。