Как сделать GLM, где ответ зависит от ответа той же группы?

Question

У меня есть набор данных

y x1 x2 x3 
1 91 21 A1 
0 54 87 A1 
1 88 54 B1 
1 67 42 B1 
0 32 83 B1 
....... 

Я сначала хочу сделать GLM у ~ x1 + x2 на подмножестве, где x2> 50 (я могу сделать это с Подменю),

mdl1 <- glm(y~x1+x2, family=binomial, data=subset(df, x2>50)) 

я хотите сделать второй glm y ~ x1 + x2 + y (из x2> 50 и x3 - то же самое).

Например, вторая строка имеет x2 = 87 и x3 = A1, тогда это y = 0 будет иметь влияние на моделирование y для первой строки, которая находится в том же x3 = A1, и имеет x2 < 50 Аналогично, y из строк 3 и 5 будет влиять на строку 4.

Я не знаю, как включить это в модель.

mdl2 <- glm(y~x1+x2+???, family=binomial, data=subset(df, x2<=50)) 

Я заподозрить у, когда x2 < 50 зависит от других значений у одного и того же х3.

Answer 1

Я думаю, что вы могли бы соответствовать комбинированную GLM, который захватывает мысль, что «у, когда x2 < 50 зависит от других значений у одного и того же х3»:

mdl <- glm(y ~ x1 + x2 + x3:as.numeric(x2 < 50), family=binomial) 

Для x2> = 50, as.numeric(x2 < 50) равна нулю , поэтому модель имеет вид y ~ x1 + x2. Для x2 < 50, as.numeric(x2 < 50) - один, поэтому модель становится y ~ x1 + x2 + x3.

при таком подходе коэффициенты для х1 и х2 на основе всех данных, но эффекты x3 основаны только на данных с x2 < 50.