"BD" <-
structure(list(KCl = c(30, 30, 20, 50, 25, 10, 40, 15, 50, 10, 
20, 45, 50, 35, 40, 35, 30, 15, 40, 30, 10, 50, 20, 10, 25, 45, 
30, 25, 45, 10, 45, 10, 30, 30, 20, 50, 10, 40, 15, 20, 45, 50, 
35, 35, 30, 15, 30, 10, 50, 20, 25, 45, 30, 25, 45, 10, 45, 10, 
30, 30, 25, 50, 10, 40, 40, 15, 35, 15, 10, 40, 10, 25, 35, 35, 
30, 25, 45, 50, 40, 20, 50, 30, 15, 10, 50, 15, 25, 45, 40), 
    MgCl2 = c(5, 5, 4, 7, 7, 4, 4, 6, 6, 7, 4, 7, 5, 6, 5, 7, 
    7, 6, 3, 4, 4, 4, 3, 3, 3, 4, 6, 3, 6, 5, 3, 6, 5, 5, 4, 
    7, 4, 4, 6, 4, 7, 5, 6, 7, 7, 6, 4, 4, 4, 3, 3, 4, 6, 3, 
    6, 5, 3, 6, 5, 5, 7, 6, 7, 5, 3, 6, 3, 5, 6, 6, 7, 7, 7, 
    5, 4, 6, 5, 6, 3, 7, 7, 4, 3, 4, 3, 3, 4, 4, 7), KPO4 = c(25, 
    25, 20, 20, 30, 25, 20, 45, 35, 25, 40, 30, 45, 20, 40, 45, 
    40, 25, 35, 45, 30, 40, 30, 20, 45, 30, 30, 35, 20, 40, 40, 
    35, 25, 25, 20, 20, 25, 20, 45, 40, 30, 45, 20, 45, 40, 25, 
    45, 30, 40, 30, 45, 30, 30, 35, 20, 40, 40, 35, 25, 25, 30, 
    35, 25, 40, 35, 40, 25, 20, 45, 25, 20, 35, 40, 30, 35, 40, 
    35, 30, 45, 40, 25, 20, 30, 45, 40, 30, 25, 25, 25), dNTP = c(625, 
    625, 1500, 250, 1500, 1250, 1250, 1250, 1500, 1250, 250, 
    1250, 1000, 1000, 1500, 750, 1000, 250, 250, 1500, 1000, 
    750, 750, 500, 500, 500, 250, 1000, 500, 500, 1250, 750, 
    625, 625, 1500, 250, 1250, 1250, 1250, 250, 1250, 1000, 1000, 
    750, 1000, 250, 1500, 1000, 750, 750, 500, 500, 250, 1000, 
    500, 500, 1250, 750, 625, 625, 1500, 1500, 1250, 1500, 250, 
    1500, 1500, 750, 250, 1500, 750, 1250, 500, 1250, 1250, 1000, 
    250, 500, 1000, 750, 1000, 250, 500, 1250, 750, 1000, 500, 
    750, 500), lnya = c(12.904207, 12.672946, 9.172639, 9.86786, 
    9.87817, 6.423247, 6.131226, 6.011267, 8.44247, 12.072541, 
    11.252859, 9.088173, 10.089967, 13.946539, 8.985946, 9.197255, 
    9.786954, 6.398595, 8.051978, 4.969813, 5.609472, 4.94876, 
    5.874931, 6.50279, 6.811244, 9.69892, 10.348173, 8.101678, 
    11.703546, 13.745088, 8.830543, 11.643954, 13.034624, 12.479909, 
    8.166216, 9.711116, 8.665613, 7.659171, 7.992945, 11.140411, 
    9.588777, 8.074026, 13.478638, 10.410305, 10.817776, 7.575585, 
    7.021084, 7.912057, 6.44254, 6.042633, 7.130899, 9.680344, 
    8.318742, 7.654443, 9.595603, 12.456831, 8.064636, 11.060369, 
    12.128111, 13.191889, 9.268609, 8.273847, 12.441145, 8.958025, 
    8.538955, 7.886081, 8.422883, 8.565983, 11.342137, 8.457443, 
    8.38936, 10.606585, 11.3679, 8.665613, 8.773385, 9.384294, 
    9.78132, 12.25009, 9.510445, 13.311329, 11.14331, 9.441452, 
    9.056023, 8.846497, 8.76873, 9.130214, 12.657148, 9.239899, 
    10.210972)), .Names = c("KCl", "MgCl2", "KPO4", "dNTP", "lnya"
), 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", "38", "39", 
"40", "43", "44", "45", "46", "48", "49", "50", "52", "53", "54", 
"55", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", 
"67", "68", "69", "70", "71", "72", "73", "75", "76", "77", "78", 
"79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", 
"90", "91", "92", "93", "94", "95", "96"))