Closed mwilson19 closed 4 years ago
I'm getting this, so I suppose it is fixed. Please use reprex to prepare example code, I had to load tidyverse and broom in addition to the code you shared to make it work:
library(tidyverse)
library(broom)
library(dslabs)
data(gapminder)
gap_nested <- gapminder %>%
group_by(country) %>%
nest()
gap_models <- gap_nested %>%
mutate(model = map(data, ~lm(formula = life_expectancy~year, data = .x)))
model_perf_nested <- gap_models %>%
mutate(fit = map(model, ~glance(.x)))
model_perf <- model_perf_nested %>%
unnest(fit)
model_perf %>%
slice_max(n = 5, order_by = r.squared)
#> # A tibble: 185 x 15
#> # Groups: country [185]
#> country data model r.squared adj.r.squared sigma statistic p.value df
#> <fct> <lis> <lis> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 Albania <tib… <lm> 0.963 0.963 0.772 1443. 3.72e-41 1
#> 2 Algeria <tib… <lm> 0.938 0.937 2.53 828. 7.84e-35 1
#> 3 Angola <tib… <lm> 0.989 0.989 0.698 5091. 6.69e-56 1
#> 4 Antigu… <tib… <lm> 0.938 0.937 0.977 828. 7.64e-35 1
#> 5 Argent… <tib… <lm> 0.983 0.982 0.479 3103. 4.58e-50 1
#> 6 Armenia <tib… <lm> 0.288 0.275 1.57 22.2 1.70e- 5 1
#> 7 Aruba <tib… <lm> 0.882 0.880 0.964 412. 3.28e-27 1
#> 8 Austra… <tib… <lm> 0.983 0.983 0.540 3240. 1.42e-50 1
#> 9 Austria <tib… <lm> 0.989 0.989 0.430 4949. 1.45e-55 1
#> 10 Azerba… <tib… <lm> 0.679 0.673 1.54 116. 3.48e-15 1
#> # … with 175 more rows, and 6 more variables: logLik <dbl>, AIC <dbl>,
#> # BIC <dbl>, deviance <dbl>, df.residual <int>, nobs <int>
Created on 2020-08-10 by the reprex package (v0.3.0.9001)
When nesting models then unnesting the slice function doesn't appear to be working, however on the last bit of code below add
as_tibble() %>%
and it works. Not sure why that is...Updated all packages and restarted R (tidyverse, tidymodels)