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