Всем привет! есть у кого идеи почему t1 %>% filter(

Question

Всем привет! есть у кого идеи почему t1 %>% filter(

) %>% arrange(client_code) %>%
mutate(pred=
lm(data=.,sum~dif/factor(client_code))$fitted.values
)

и
t1 %>% group_by(client_code) %>% nest() %>%
mutate(pred= lapply(data,function(df) lm(sum~dif,data=df)$fitted.values)) %>% unnest()

Дают разные результаты в предсказании?
При этом для не для всех а для отдельных ошибаются.
Оставил второй вариант - но хотелось бы разобраться конечно

#programming #r #russian

0

26.12.2020

6 ответов

22 просмотра

Ivan Struzhkov Автор вопроса

Stas
Предположу что сортировка не равна группировке, и ...

Звучит похоже конечно , хотя я так и не нащупал где проблема. Но мысли появились куда смотреть еще, спасибо

0

26.12.2020

Evgenii Kalenkovich

Вы в первом случай одну модель на всех фитите, а во втором - по отдельности на каждого клиента. Вот и предсказания разные. Или я что-то упускаю?

0

26.12.2020

Ivan Struzhkov Автор вопроса

Evgenii Kalenkovich
Вы в первом случай одну модель на всех фитите, а в...

предполагается что /factor(client_code) должен разделять как раз внутри lm()

0

26.12.2020

Grigoriy Demin

А что обозначает знак деления в формуле? Я поискал в документации - ничего не нашел. Если смотреть на результаты регрессии со знаком деления, то коэффициенты разные в зависимости от группы, но константа для всех уровней фактора одинаковая. То есть это не то же самое, что по отдельности делать. По смыслу должно быть *. Извиняюсь за аццкий однострочник: summary( lm(Sepal.Length ~ Sepal.Width*Species, data = iris)$fitted.values - unlist(lapply(lapply(split(iris, iris$Species), lm, formula = Sepal.Length ~ Sepal.Width), `[[`, 'fitted.values')) ) Разница не более 1e-15

0

27.12.2020

Ivan Struzhkov Автор вопроса

Grigoriy Demin
А что обозначает знак деления в формуле? Я поискал...

В синтетических датасетах все так. Но вот в практике - когда на одном и том же датасете делаю два подхода - вот через nest()

0

27.12.2020

Stas · Accepted Answer

Stas

Предположу что сортировка не равна группировке, и где-нибудь здесь кроется ошибка

0

26.12.2020

169 похожих чатов

Всем привет! есть у кого идеи почему t1 %>% filter(

6 ответов

Похожие вопросы