
18.4
フォーミュラとモデルファミリー
325
18.4.3
交互作用(
2
連続変数)
2つの連続変数の等価モデルを調べます。最初は、先ほどの例と同じようになります。
mod1 <- lm(y ~ x1 + x2, data = sim4)
mod2 <- lm(y ~ x1 * x2, data = sim4)
grid <- sim4 %>%
data_grid(
x1 = seq_range(x1, 5),
x2 = seq_range(x2, 5)
) %>%
gather_predictions(mod1, mod2)
grid
#> # A tibble: 50
×
4
#> model x1 x2 pred
#> <chr> <dbl> <dbl> <dbl>
#> 1 mod1 -1.0 -1.0 0.996
#> 2 mod1 -1.0 -0
.5 -0.395
#> 3 mod1 -1.0 0.0 -1.786
#> 4 mod1 -1.0 0.5 -3.177
#> 5 mod1 -1.0 1.0 -4.569
#> 6 mod1 -0.5 -1.0 1.907
#> # ... with 44 more rows
data_grid()
内部での私の
seq_range()
の使い方に注意します。
x
の一意な値すべてを使う代わり
に、最小数と最大数の間の5つの格子値を私は等間隔に使います。これはそれほど重要なことではあ
りませんが、一般に役立つ技法です。
seq_range()
には他に3 つの役立つ引数があ ...