Slide 77
Slide 77 text
dat <- data.frame(x = runif(N, 0, 10)) %>%
nest %>% list %>% rep(3) %>% bind_rows %>%
rowid_to_column("id") %>%
mutate(s = c(5, 10, 25)) %>%
mutate(data = map2(data, s,
~mutate(.x, y = a * x + b + rnorm(N, 0, .y)))) %>%
mutate(lm = map(data, ~lm(.$y ~ .$x)),
predict = map(lm, predict)) %>%
mutate(r_sq = map_dbl(lm, ~summary(.)$r.squared),
a = map_dbl(lm, ~.$coeff[2]),
b = map_dbl(lm, ~.$coeff[1]))
Nested data modeling