The R Project SVN R

### $Id: Pixel.R,v 1.1 1999/11/11 21:09:49 bates Exp $
### X-ray pixel intensity for several dogs over time
### Provided by Deborah Darien, Department of Medical Sciences, School of
### Veterinary Medicine, U. of Wisconsin - Madison.
"Pixel" <-
  structure(list(
  Dog = structure(factor(c(1, 1, 1, 1, 1, 1, 1, 
    3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 
    6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 10, 10, 
    2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 
    4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 
    8, 8, 8, 9, 9, 9, 9, 10, 10, 2, 2, 2), levels=1:10),
    class = "factor",
    .Label = c("1", "10", "2", "3", "4", "5", "6", "7", "8", "9")),
Side = structure(factor(c(2, 
  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
  2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), levels=1:2),
  class = "factor", .Label = c("L", "R")),
day = c(0, 1, 2, 4, 6, 10, 14, 0, 1, 2, 4, 6, 10, 14, 
  0, 4, 6, 8, 10, 14, 21, 0, 4, 6, 8, 10, 14, 21, 4, 6, 8, 10, 
  14, 4, 6, 8, 10, 14, 4, 6, 10, 14, 4, 6, 10, 14, 4, 8, 4, 6, 
  8, 0, 1, 2, 4, 6, 10, 14, 0, 1, 2, 4, 6, 10, 14, 0, 4, 6, 8, 
  10, 14, 21, 0, 4, 6, 8, 10, 14, 21, 4, 6, 8, 10, 14, 4, 6, 8, 
  10, 14, 4, 6, 10, 14, 4, 6, 10, 14, 4, 8, 4, 6, 8),
pixel = c(1045.8, 
  1044.5, 1042.9, 1050.4, 1045.2, 1038.9, 1039.8, 1041.8, 1045.6, 
  1051, 1054.1, 1052.7, 1062.0, 1050.8, 1039.8, 1082.9, 1091.1, 
  1088.5, 1099.4, 1086.8, 1073.6, 1034.1, 1090.7, 1126.4, 1108.8, 
  1122.5, 1100.7, 1102.2, 1124.4, 1128.2, 1124.1, 1126.7, 1109.9, 
  1094.2, 1084.1, 1093.9, 1082.7, 1076.1, 1151.3, 1160.3, 1139.7, 
  1133.2, 1077.5, 1085.9, 1100.8, 1083.4, 1068.2, 1064.0, 1095.0, 
  1094.7, 1081.4, 1040.9, 1054.1, 1051.9, 1054.1, 1051.3, 1051.4, 
  1046.1, 1053.6, 1052.9, 1053.6, 1057.1, 1064, 1065.3, 1053.4, 
  1043.9, 1069, 1082.5, 1077.7, 1073.5, 1076.4, 1053.8, 1039.9, 
  1114.1, 1111.5, 1100.3, 1124.0, 1113.2, 1077.0, 1114.7, 1117.6, 
  1128.2, 1130.5, 1102.2, 1106.0, 1121.5, 1115.5, 1105.2, 1092.8, 
  1129.7, 1136.2, 1108.5, 1094.3, 1098.8, 1103.6, 1106.8, 1105.0, 
  1097.2, 1099.5, 1132.3, 1154.5, 1161.1)),
row.names = 1:102,
class = c("nmGroupedData", "groupedData", "data.frame"),
formula = pixel ~ day | Dog/Side,
formulaList = list(Dog = ~Dog, Side = ~Side),
labels = list(x = "Time post injection", y = "Pixel intensity"),
units = list(x = "(days)"),
order.groups = list(Dog = TRUE, Side = TRUE),
FUN = function (x) max(x, na.rm = TRUE))