Помогите, пожалуйста. Есть 2 группирующие переменные. Пол (м и ж) и курение(курит не курит). нужно провести регрессию стратифицированно, т.е. для 4х групп.
м курит
м не курит
ж курит
э не курит

сторою модельку
glm.linear <- glm(Emot ~ color)
Но как мне сделать, чтобы она построилась для 4 групп отдельно. Все что удалось нагуглить, так это то, что это делается через функцию lapply. А как именно я не поняла.

А в чем проблема "по взрослому" посчитать модель Emot ~ color + sex + smoke ? Или выбрать оптимальную модель из комбинации взаимодействий этих факторов.