Describing the problem - Linear Example

Creating a matrix of means for each combination.

y_mean_matrix <-
  tapply(yx, list(x1c,x2c), mean) |>
  round(digits = 2)

y_mean_matrix

##           (0,0.2] (0.2,0.4] (0.4,0.6] (0.6,0.8] (0.8,1]
## (0,0.2]      0.19      0.42      0.63      0.80    0.99
## (0.2,0.4]    0.39      0.59      0.81      1.01    1.19
## (0.4,0.6]    0.59      0.81      0.98      1.20    1.44
## (0.6,0.8]    0.76      1.00      1.20      1.40    1.58
## (0.8,1]      1.01      1.22      1.38      1.58    1.77