options ls=75 ps=55;

data fl_depr;
INPUT fluoride natural deprived t_decay;

district = fluoride;
if natural eq 1 then district = district + 1; /* 0 , 1, 2 */

LINES;
0  .  -23   1.4
0  .  -15   0.3
0  .  -15   1.4
0  .  -14   1.6
0  .  -13   1.65
0  .  -12   1.49
0  .  -10.5 1.25
0  .  -8    1.7
0  .  -6    1.1
0  .  -3    1.5
0  .  -1    2.1
0  .  -1    2.3
0  .  -0.5  1.6
0  .   1    1.1
0  .   1    1.7
0  .   1.5  2.2
0  .   2    2.3
0  .   4    3.2
0  .   5    2
0  .   5.5  2.35
0  .   9.5  2.3
0  .  10    2.1
0  .  11    2.6
0  .  12    2.4
0  .  13    2.5
0  .  15    2.9
0  .  15    2.3
0  .  19    2.7
0  .  20    2.5
0  .  20.5  2.9
0  .  22    2.8
0  .  23    2.9
0  .  23    3
0  .  23    3.1
0  .  33    4.2
0  .  36    3.4
0  .  37    3.3
0  .  40    2.65
0  .  41    4.2
0  .  41.5  4.3
0  .  48    3.3
1  0 -17.5  1.4
1  0  -5    0.55
1  0  -5    0.9
1  0  -4    0.6
1  0   1    1
1  0   2    1
1  0   3    1.3
1  0   3    1.1
1  0   5    1.05
1  0   5    0.35
1  0   8    0.7
1  0   9    1
1  0   9    1.25
1  0  10    2
1  0  11    1.8
1  0  11.5  1.5
1  0  12    1.7
1  0  18    1.2
1  0  20    1.3
1  0  22    1.3
1  0  22    1.8
1  0  22.5  1.2
1  0  26    1.2
1  0  30    1.3
1  0  32    1.4
1  0  32    2.5
1  0  39    1.8
1  0  41    1.3
1  0  42    2.2
1  0  42    1.5
1  0  44    1.4
1  0  53    1.7
1  0  56    1.8
1  1 -16    0.4
1  1  -12   0.6
1  1  -9    0.7
1  1  -8.5  0.4
1  1  -8.5  0.55
1  1  -2    0.7
1  1  -1    0.8
1  1   6    1
1  1  19    0.8
1  1  21    0.7
1  1  23    0.75
1  1  25    0.9
1  1  28    0.8
1  1  32    0.55
1  1  38    0.8
1  1  40    1
;
run;


proc plot data=fl_depr;
  plot t_decay*deprived=district;
run;