(all_y)*1.2];

# Do regression on each pair and plot the input and results
figure;
haxs = NaN (1, 4);
for i_pair = 1:4
  x = data(i_pair).x;
  y = data(i_pair).y;
  # TODO: Port the anova and other characterizations from the R code
  # TODO: Do a linear regression and plot its line
  hax = subplot (2, 2, i_pair);
  haxs(i_pair) = hax;
  xlabel (sprintf (&quot;x%d&quot;, i_pair));
  ylabel (sprintf (&quot;y%d&quot;, i_pair));
  scatter (x, y, &quot;r&quot;);
endfor

# Fiddle with the plot axes parameters
lin