Я подбираю линейные модели для нескольких разных периодов времени (в приведенном ниже примере используются страны). У меня два вопроса.
Вопрос 1. Есть ли способ передать столбец mods
в ggpredict
, не начиная с df1$mods
? что-то вроде df1 %>% select(mods)%>% map(., ggpredict)
. Я действительно пробовал это, но это не сработало. Это не имеет большого значения; но мне любопытно.
Вопрос 2: Есть ли способ автоматизировать получение имен для категорий группирующей переменной и сопоставление их с предсказанными значениями?
Спасибо!
library(tidyverse)
library(ggeffects)
#make a fake data frame
df1<-data.frame(country=rep(c("A", "B"), 100), var1=rnorm(200), var2=rnorm(200))
df1
df1 %>%
group_by(country) %>%
#in reality I have 10 to 12 variables and 1 grouping variable with 12 categories
nest(var1, var2) %>%
#in reality I'm also doing glms, but I don't think it matters
mutate(mods=map(data, function(x) lm(var1~var2, data=x))) ->out
out$mods %>%
map_df(., ggpredict, terms=c('var2 [0,0.5]')) %>%
#Is there a way to automate this line; by taking the values of the grouping variables
#somehow.
mutate(country=rep(c('A', 'B'), each=2))