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

ну и по смыслу у вас же не 300 гипотез, гиптотеза одна, в данных есть просто переменная - город

Спасибо большое за ответ! Я опасалась вводить город как предиктор, потому что города очень разные, и я совсем не предполагаю, что зависимость А от В есть везде. Она есть, скажем, в 20 городах из 300, и я как раз хочу их вычленить. Модель с двумя факторами, насколько я понимаю, покажет, что зависимости А от В нет, раз она так редко встречается.

А уровень освещенности в рамках города варьируется?

Ну то есть, если вы модели по городам строите, то должен варьироваться

Можно попробовать так:

lm(a ~ b + city:b + city)

Первый коэффициент - в целом зависимость a от b, а значимые коэффиценты во втором блоке будут показывать, что зависимость a от b может варьироваться между городами

Да, варьируется!

lm(a ~ b + city:b + city)

Спасибо за совет! Кажется, это должно подойти!

Это известная проблема, про которую можно почитать, например, в этой книжке в начале третьей части. Как уже сказали выше, нужно, чтобы Вы четко сформулировали, что Вы хотите сказать:

1) согласно нашим данным про связь переменных y и x мы можем выделить группу городов А и группу городов Б (кластеризация или empirical bayes estimation)
2) нас интересует связь переменных y и x, какая-то такая общая, игнорирующая информацию про отдельные города (mixed effects model)

книжка, кстати одна из любимых, всем рекоммендую

Супер, спасибо! Меня интересует именно первый пункт.

Недавно её вспоминал. Её можно использовать, как учебник?

я бы считал ее полезным чтением

Просто никак не могу понять ее ЦА.

Скорее для расширения кругозора. Читать интересно, но большинство вещей с птичьего полета

А почему бы multilevel не попробовать?

Тут же задача именно на каких электродах/городах есть эффект, а не есть ли сам эффект

Мне кажется, да, все корректно, если не пугает либеральность fdr подхода в целом

Спасибо!

Коллеги, а у кого-нибудь есть опыт участия в соревнованиях на Kaggle на R?
Как это работает? Обучаешь модели на трейне, потом на тесте делаешь предсказания и заливаешь их на Kaggle? Или нужно саму модель залить?

модель - если выиграл, попросят показать