1### $Id: DNase.R,v 1.1 2003/12/11 07:16:04 ripley Exp $ 2### Data on calibration runs for an assay of the enzyme DNase. 3### cited in Davidian and Giltinan (1995) section 5.2.4, p.134 4"DNase" <- 5 structure(list( 6 Run = structure(ordered(c(4, 4, 4, 4, 4, 4, 4, 7 4, 4, 4, 4, 4, 4, 4, 4, 4, 10, 10, 10, 10, 10, 10, 10, 10, 10, 8 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 9 11, 11, 11, 11, 11, 11, 11, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 10 5, 5, 5, 5, 5, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 11 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 8, 8, 8, 8, 8, 12 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 13 6, 6, 6, 6, 6, 6, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 14 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 15 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2), levels=1:11), 16 class = c("ordered", "factor"), 17 .Label = c("10", "11", "9", "1", "4", "8", "5", "7", "6", "2", "3")), 18 conc = c(0.04882812, 0.04882812, 0.1953125, 19 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 1.5625, 20 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 21 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 22 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 23 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 24 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 25 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 26 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 27 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 28 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 29 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 30 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 31 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 32 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 33 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 34 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 35 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 36 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 37 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 38 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5, 0.04882812, 0.04882812, 39 0.1953125, 0.1953125, 0.390625, 0.390625, 0.78125, 0.78125, 1.5625, 40 1.5625, 3.125, 3.125, 6.25, 6.25, 12.5, 12.5), 41 density = c(0.017, 42 0.018, 0.121, 0.124, 0.206, 0.215, 0.377, 0.374, 0.614, 0.609, 43 1.019, 1.001, 1.334, 1.364, 1.73, 1.71, 0.045, 0.05, 0.137, 0.123, 44 0.225, 0.207, 0.401, 0.383, 0.672, 0.681, 1.116, 1.078, 1.554, 45 1.526, 1.932, 1.914, 0.07, 0.068, 0.173, 0.165, 0.277, 0.248, 46 0.434, 0.426, 0.703, 0.689, 1.067, 1.077, 1.629, 1.479, 2.003, 47 1.884, 0.011, 0.016, 0.118, 0.108, 0.2, 0.206, 0.364, 0.36, 0.62, 48 0.64, 0.979, 0.973, 1.424, 1.399, 1.74, 1.732, 0.035, 0.035, 49 0.132, 0.135, 0.224, 0.22, 0.385, 0.39, 0.658, 0.647, 1.06, 1.031, 50 1.425, 1.409, 1.75, 1.738, 0.086, 0.103, 0.191, 0.189, 0.272, 51 0.277, 0.44, 0.426, 0.686, 0.676, 1.062, 1.072, 1.424, 1.459, 52 1.768, 1.806, 0.094, 0.092, 0.182, 0.182, 0.282, 0.273, 0.444, 53 0.439, 0.686, 0.668, 1.052, 1.035, 1.409, 1.392, 1.759, 1.739, 54 0.054, 0.054, 0.152, 0.148, 0.226, 0.222, 0.392, 0.383, 0.658, 55 0.644, 1.043, 1.002, 1.466, 1.381, 1.743, 1.724, 0.032, 0.043, 56 0.142, 0.155, 0.239, 0.242, 0.42, 0.395, 0.624, 0.705, 1.046, 57 1.026, 1.398, 1.405, 1.693, 1.729, 0.052, 0.094, 0.164, 0.166, 58 0.259, 0.256, 0.439, 0.439, 0.69, 0.701, 1.042, 1.075, 1.34, 59 1.406, 1.699, 1.708, 0.047, 0.057, 0.159, 0.155, 0.246, 0.252, 60 0.427, 0.411, 0.704, 0.684, 0.994, 0.98, 1.421, 1.385, 1.715, 61 1.721)), 62row.names = 1:176, 63class = c("nfnGroupedData", "nfGroupedData", "groupedData", "data.frame"), 64formula = density ~ conc | Run, 65labels = list(x = "DNase concentration", y = "Optical density"), 66units = list(x = "(ng/ml)")) 67environment(attr(DNase, "formula")) <- emptyenv() 68