私はデータセットを持っています
y x1 x2 x3
1 91 21 A1
0 54 87 A1
1 88 54 B1
1 67 42 B1
0 32 83 B1
.......
最初に、x2>50 のサブセットで glm y~x1+x2 を実行したい (サブセット化でこれを実行できます)。
mdl1 <- glm(y~x1+x2, family=binomial, data=subset(df, x2>50))
次に、2 番目の glm y~x1+x2+ y (x2>50 および x3 は同じ) を実行します。
たとえば、2 行目に x2=87 と x3=A1 がある場合、y=0 は、同じ x3=A1 にあり、x2<50 である最初の行の y のモデリングに影響します。同様に、3 行目と 5 行目の y は 4 行目に影響を与えます。
これをモデルに含める方法がわかりません。
mdl2 <- glm(y~x1+x2+???, family=binomial, data=subset(df, x2<=50))
x2<50 が同じ x3 の他の y に依存している場合、y を疑っています。