lattice: Trellis Graphics for R

"ethanol" <-
structure(list(NOx = c(3.741, 2.295, 1.498, 2.881, 0.76, 3.12, 
0.638, 1.17, 2.358, 0.606, 3.669, 1, 0.981, 1.192, 0.926, 1.59, 
1.806, 1.962, 4.028, 3.148, 1.836, 2.845, 1.013, 0.414, 0.812, 
0.374, 3.623, 1.869, 2.836, 3.567, 0.866, 1.369, 0.542, 2.739, 
1.2, 1.719, 3.423, 1.634, 1.021, 2.157, 3.361, 1.39, 1.947, 0.962, 
0.571, 2.219, 1.419, 3.519, 1.732, 3.206, 2.471, 1.777, 2.571, 
3.952, 3.931, 1.587, 1.397, 3.536, 2.202, 0.756, 1.62, 3.656, 
2.964, 3.76, 0.672, 3.677, 3.517, 3.29, 1.139, 0.727, 2.581, 
0.923, 1.527, 3.388, 2.085, 0.966, 3.488, 0.754, 0.797, 2.064, 
3.732, 0.586, 0.561, 0.563, 0.678, 0.37, 0.53, 1.9), C = c(12, 
12, 12, 12, 12, 9, 9, 9, 12, 12, 12, 12, 15, 18, 7.5, 12, 12, 
15, 15, 9, 9, 7.5, 7.5, 18, 18, 15, 15, 7.5, 7.5, 9, 15, 15, 
15, 15, 15, 9, 9, 7.5, 7.5, 7.5, 18, 18, 18, 18, 9, 9, 9, 9, 
7.5, 7.5, 7.5, 15, 18, 18, 15, 15, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 
7.5, 18, 18, 18, 12, 12, 9, 9, 9, 15, 15, 15, 15, 15, 7.5, 7.5, 
9, 7.5, 18, 18, 7.5, 9, 12, 15, 18, 18), E = c(0.907, 0.761, 
1.108, 1.016, 1.189, 1.001, 1.231, 1.123, 1.042, 1.215, 0.93, 
1.152, 1.138, 0.601, 0.696, 0.686, 1.072, 1.074, 0.934, 0.808, 
1.071, 1.009, 1.142, 1.229, 1.175, 0.568, 0.977, 0.767, 1.006, 
0.893, 1.152, 0.693, 1.232, 1.036, 1.125, 1.081, 0.868, 0.762, 
1.144, 1.045, 0.797, 1.115, 1.07, 1.219, 0.637, 0.733, 0.715, 
0.872, 0.765, 0.878, 0.811, 0.676, 1.045, 0.968, 0.846, 0.684, 
0.729, 0.911, 0.808, 1.168, 0.749, 0.892, 1.002, 0.812, 1.23, 
0.804, 0.813, 1.002, 0.696, 1.199, 1.03, 0.602, 0.694, 0.816, 
1.037, 1.181, 0.899, 1.227, 1.18, 0.795, 0.99, 1.201, 0.629, 
0.608, 0.584, 0.562, 0.535, 0.655)), class = "data.frame", row.names = c("1", 
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", 
"14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", 
"25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", 
"36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", 
"47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", 
"58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", 
"69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", 
"80", "81", "82", "83", "84", "85", "86", "87", "88"), .Names = c("NOx", 
"C", "E"))