R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.1.1  (2005-06-20), ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

> ### * <HEADER>
> ###
> attach(NULL, name = "CheckExEnv")
> assign(".CheckExEnv", as.environment(2), pos = length(search())) # base
> ## add some hooks to label plot pages for base and grid graphics
> setHook("plot.new", ".newplot.hook")
> setHook("persp", ".newplot.hook")
> setHook("grid.newpage", ".gridplot.hook")
> 
> assign("cleanEx",
+        function(env = .GlobalEnv) {
+ 	   rm(list = ls(envir = env, all.names = TRUE), envir = env)
+            RNGkind("default", "default")
+ 	   set.seed(1)
+    	   options(warn = 1)
+ 	   delayedAssign("T", stop("T used instead of TRUE"),
+ 		  assign.env = .CheckExEnv)
+ 	   delayedAssign("F", stop("F used instead of FALSE"),
+ 		  assign.env = .CheckExEnv)
+ 	   sch <- search()
+ 	   newitems <- sch[! sch %in% .oldSearch]
+ 	   for(item in rev(newitems))
+                eval(substitute(detach(item), list(item=item)))
+ 	   missitems <- .oldSearch[! .oldSearch %in% sch]
+ 	   if(length(missitems))
+ 	       warning("items ", paste(missitems, collapse=", "),
+ 		       " have been removed from the search path")
+        },
+        env = .CheckExEnv)
> assign("..nameEx", "__{must remake R-ex/*.R}__", env = .CheckExEnv) # for now
> assign("ptime", proc.time(), env = .CheckExEnv)
> grDevices::postscript("Devore6-Examples.ps")
> assign("par.postscript", graphics::par(no.readonly = TRUE), env = .CheckExEnv)
> options(contrasts = c(unordered = "contr.treatment", ordered = "contr.poly"))
> options(warn = 1)    
> library('Devore6')
> 
> assign(".oldSearch", search(), env = .CheckExEnv)
> assign(".oldNS", loadedNamespaces(), env = .CheckExEnv)
> cleanEx(); ..nameEx <- "ex01.11"
> 
> ### * ex01.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.11
> ### Title: data from exercise 1.11
> ### Aliases: ex01.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.11)
> str(ex01.11)
`data.frame':	40 obs. of  1 variable:
 $ Scores: int  74 89 80 93 64 67 72 70 66 85 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.12"
> 
> ### * ex01.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.12
> ### Title: data from exercise 1.12
> ### Aliases: ex01.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.12)
> str(ex01.12)
`data.frame':	36 obs. of  1 variable:
 $ SpecGrav: num  0.31 0.35 0.36 0.36 0.37 0.38 0.4 0.4 0.4 0.41 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.13"
> 
> ### * ex01.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.13
> ### Title: data from exercise 1.13
> ### Aliases: ex01.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.13)
> str(ex01.13)
`data.frame':	153 obs. of  1 variable:
 $ strength: int  122 127 130 131 132 133 134 134 135 135 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.14"
> 
> ### * ex01.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.14
> ### Title: data from exercise 1.14
> ### Aliases: ex01.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.14)
> str(ex01.14)
`data.frame':	129 obs. of  1 variable:
 $ Rate: num  4.6 12.3 7.1 7 4 9.2 6.7 6.9 11.5 5.1 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.15"
> 
> ### * ex01.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.15
> ### Title: data from exercise 1.15
> ### Aliases: ex01.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.15)
> str(ex01.15)
`data.frame':	37 obs. of  2 variables:
 $ Score: int  56 44 62 36 39 53 50 65 45 40 ...
 $ Type : Factor w/ 2 levels "Creamy","Crunchy": 1 1 1 1 1 1 1 1 1 1 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.17"
> 
> ### * ex01.17
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.17
> ### Title: data from exercise 1.17
> ### Aliases: ex01.17
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.17)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.18"
> 
> ### * ex01.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.18
> ### Title: data from exercise 1.18
> ### Aliases: ex01.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.19"
> 
> ### * ex01.19
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.19
> ### Title: data from exercise 1.19
> ### Aliases: ex01.19
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.19)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.20"
> 
> ### * ex01.20
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.20
> ### Title: data from exercise 1.20
> ### Aliases: ex01.20
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.20)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.21"
> 
> ### * ex01.21
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.21
> ### Title: data from exercise 1.21
> ### Aliases: ex01.21
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.21)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.23"
> 
> ### * ex01.23
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.23
> ### Title: data from exercise 1.23
> ### Aliases: ex01.23
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.23)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.24"
> 
> ### * ex01.24
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.24
> ### Title: data from exercise 1.24
> ### Aliases: ex01.24
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.24)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.25"
> 
> ### * ex01.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.25
> ### Title: data from exercise 1.25
> ### Aliases: ex01.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.27"
> 
> ### * ex01.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.27
> ### Title: data from exercise 1.27
> ### Aliases: ex01.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.27)
> str(ex01.27)
`data.frame':	50 obs. of  1 variable:
 $ lifetime: int  11 14 20 23 31 36 39 44 47 50 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.28"
> 
> ### * ex01.28
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.28
> ### Title: data from exercise 1.28
> ### Aliases: ex01.28
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.28)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.29"
> 
> ### * ex01.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.29
> ### Title: data from exercise 1.29
> ### Aliases: ex01.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.32"
> 
> ### * ex01.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.32
> ### Title: data from exercise 1.32
> ### Aliases: ex01.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.33"
> 
> ### * ex01.33
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.33
> ### Title: data from exercise 1.33
> ### Aliases: ex01.33
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.33)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.34"
> 
> ### * ex01.34
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.34
> ### Title: data from exercise 1.34
> ### Aliases: ex01.34
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.34)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.35"
> 
> ### * ex01.35
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.35
> ### Title: data from exercise 1.35
> ### Aliases: ex01.35
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.35)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.36"
> 
> ### * ex01.36
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.36
> ### Title: data from exercise 1.36
> ### Aliases: ex01.36
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.36)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.37"
> 
> ### * ex01.37
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.37
> ### Title: data from exercise 1.37
> ### Aliases: ex01.37
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.37)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.38"
> 
> ### * ex01.38
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.38
> ### Title: data from exercise 1.38
> ### Aliases: ex01.38
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.38)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.39"
> 
> ### * ex01.39
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.39
> ### Title: data from exercise 1.39
> ### Aliases: ex01.39
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.39)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.43"
> 
> ### * ex01.43
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.43
> ### Title: data from exercise 1.43
> ### Aliases: ex01.43
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.43)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.44"
> 
> ### * ex01.44
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.44
> ### Title: data from exercise 1.44
> ### Aliases: ex01.44
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.44)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.45"
> 
> ### * ex01.45
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.45
> ### Title: data from exercise 1.45
> ### Aliases: ex01.45
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.45)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.46"
> 
> ### * ex01.46
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.46
> ### Title: data from exercise 1.46
> ### Aliases: ex01.46
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.46)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.49"
> 
> ### * ex01.49
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.49
> ### Title: data from exercise 1.49
> ### Aliases: ex01.49
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.49)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.50"
> 
> ### * ex01.50
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.50
> ### Title: data from exercise 1.50
> ### Aliases: ex01.50
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.50)
> str(ex01.50)
`data.frame':	27 obs. of  1 variable:
 $ awards: int  37 60 75 115 135 140 149 150 238 290 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.51"
> 
> ### * ex01.51
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.51
> ### Title: data from exercise 1.51
> ### Aliases: ex01.51
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.51)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.54"
> 
> ### * ex01.54
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.54
> ### Title: data from exercise 1.54
> ### Aliases: ex01.54
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.54)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.56"
> 
> ### * ex01.56
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.56
> ### Title: data from exercise 1.56
> ### Aliases: ex01.56
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.56)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.59"
> 
> ### * ex01.59
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.59
> ### Title: data from exercise 1.59
> ### Aliases: ex01.59
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.59)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.60"
> 
> ### * ex01.60
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.60
> ### Title: data from exercise 1.60
> ### Aliases: ex01.60
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.60)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.63"
> 
> ### * ex01.63
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.63
> ### Title: data from exercise 1.63
> ### Aliases: ex01.63
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.63)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.64"
> 
> ### * ex01.64
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.64
> ### Title: data from exercise 1.64
> ### Aliases: ex01.64
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.64)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.65"
> 
> ### * ex01.65
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.65
> ### Title: data from exercise 1.65
> ### Aliases: ex01.65
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.65)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.67"
> 
> ### * ex01.67
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.67
> ### Title: data from exercise 1.67
> ### Aliases: ex01.67
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.67)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.70"
> 
> ### * ex01.70
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.70
> ### Title: data from exercise 1.70
> ### Aliases: ex01.70
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.70)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.72"
> 
> ### * ex01.72
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.72
> ### Title: data from exercise 1.72
> ### Aliases: ex01.72
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.72)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.73"
> 
> ### * ex01.73
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.73
> ### Title: data from exercise 1.73
> ### Aliases: ex01.73
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.73)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.75"
> 
> ### * ex01.75
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.75
> ### Title: data from exercise 1.75
> ### Aliases: ex01.75
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.75)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.77"
> 
> ### * ex01.77
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.77
> ### Title: data from exercise 1.77
> ### Aliases: ex01.77
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.77)
> str(ex01.77)
`data.frame':	46 obs. of  1 variable:
 $ Time: num  0.2 0.3 0.5 0.5 0.5 0.6 0.6 0.7 0.7 0.7 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.80"
> 
> ### * ex01.80
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.80
> ### Title: data from exercise 1.80
> ### Aliases: ex01.80
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.80)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex01.83"
> 
> ### * ex01.83
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex01.83
> ### Title: data from exercise 1.83
> ### Aliases: ex01.83
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex01.83)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.82"
> 
> ### * ex04.82
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.82
> ### Title: data from exercise 4.82
> ### Aliases: ex04.82
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.82)
> attach(ex04.82)
> boxplot(lifetime, ylab = "Lifetime (hr)",
+         main = "Bearing lifetimes from exercise 4.82",
+         col = "lightgray")
> ## Normal probability plot on the original time scale
> qqnorm(lifetime, ylab = "Lifetime (hr)", las = 1)
> qqline(lifetime)
> ## Try normal probability plot of the log(lifetime)
> qqnorm(log(lifetime), ylab = "log(lifetime) (log(hr))",
+        las = 1)
> qqline(log(lifetime))
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.83"
> 
> ### * ex04.83
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.83
> ### Title: data from exercise 4.83
> ### Aliases: ex04.83
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.83)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.84"
> 
> ### * ex04.84
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.84
> ### Title: data from exercise 4.84
> ### Aliases: ex04.84
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.84)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.86"
> 
> ### * ex04.86
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.86
> ### Title: data from exercise 4.86
> ### Aliases: ex04.86
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.86)
> attach(ex04.86)
> boxplot(loadlife, ylab = "Load-Life (million revs)",
+         main = "Bearing load-lifes from exercise 4.86",
+         col = "lightgray")
> ## Normal probability plot
> qqnorm(loadlife, ylab = "Load-life (million revs)")
> qqline(loadlife)
> ## Weibull probability plot
> plot(
+   log(-log(1 - (seq(along = loadlife) - 0.5)/length(loadlife))),
+   log(sort(loadlife)), xlab = "Theoretical Quantiles",
+   ylab = "log(load-life) (log(million revs))",
+   main = "Weibull Q-Q Plot", las = 1)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.88"
> 
> ### * ex04.88
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.88
> ### Title: data from exercise 4.88
> ### Aliases: ex04.88
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.88)
> attach(ex04.88)
> ## Normal probability plot
> qqnorm(preciptn, ylab = "Precipitation (in)",
+    main = "Precipitation during March in Minneapolis-St. Paul")
> qqline(preciptn)
> ## Normal probability plot on square root scale
> qqnorm(sqrt(preciptn),
+    ylab = expression(sqrt("Precipitation (in)")),
+    main =
+      "Precipitation during March in Minneapolis-St. Paul")
> qqline(sqrt(preciptn))
> detach()   
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.90"
> 
> ### * ex04.90
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.90
> ### Title: data from exercise 4.90
> ### Aliases: ex04.90
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.90)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex04.91"
> 
> ### * ex04.91
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex04.91
> ### Title: data from exercise 4.91
> ### Aliases: ex04.91
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex04.91)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.01"
> 
> ### * ex06.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.01
> ### Title: data from exercise 6.1
> ### Aliases: ex06.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.01)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.02"
> 
> ### * ex06.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.02
> ### Title: data from exercise 6.2
> ### Aliases: ex06.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.02)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.03"
> 
> ### * ex06.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.03
> ### Title: data from exercise 6.3
> ### Aliases: ex06.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.04"
> 
> ### * ex06.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.04
> ### Title: data from exercise 6.4
> ### Aliases: ex06.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.05"
> 
> ### * ex06.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.05
> ### Title: data from exercise 6.5
> ### Aliases: ex06.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.05)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.06"
> 
> ### * ex06.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.06
> ### Title: data from exercise 6.6
> ### Aliases: ex06.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.09"
> 
> ### * ex06.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.09
> ### Title: data from exercise 6.9
> ### Aliases: ex06.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.15"
> 
> ### * ex06.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.15
> ### Title: data from exercise 6.15
> ### Aliases: ex06.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex06.25"
> 
> ### * ex06.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex06.25
> ### Title: data from exercise 6.25
> ### Aliases: ex06.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex06.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.10"
> 
> ### * ex07.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.10
> ### Title: data from exercise 7.10
> ### Aliases: ex07.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.26"
> 
> ### * ex07.26
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.26
> ### Title: data from exercise 7.26
> ### Aliases: ex07.26
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.26)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.33"
> 
> ### * ex07.33
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.33
> ### Title: data from exercise 7.33
> ### Aliases: ex07.33
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.33)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.37"
> 
> ### * ex07.37
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.37
> ### Title: data from exercise 7.37
> ### Aliases: ex07.37
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.37)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.45"
> 
> ### * ex07.45
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.45
> ### Title: data from exercise 7.45
> ### Aliases: ex07.45
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.45)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.46"
> 
> ### * ex07.46
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.46
> ### Title: data from exercise 7.46
> ### Aliases: ex07.46
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.46)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.47"
> 
> ### * ex07.47
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.47
> ### Title: data from exercise 7.47
> ### Aliases: ex07.47
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.47)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.49"
> 
> ### * ex07.49
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.49
> ### Title: data from exercise 7.49
> ### Aliases: ex07.49
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.49)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.56"
> 
> ### * ex07.56
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.56
> ### Title: data from exercise 7.56
> ### Aliases: ex07.56
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.56)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex07.58"
> 
> ### * ex07.58
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex07.58
> ### Title: data from exercise 7.58
> ### Aliases: ex07.58
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex07.58)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.32"
> 
> ### * ex08.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.32
> ### Title: data from exercise 8.32
> ### Aliases: ex08.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.54"
> 
> ### * ex08.54
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.54
> ### Title: data from exercise 8.54
> ### Aliases: ex08.54
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.54)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.55"
> 
> ### * ex08.55
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.55
> ### Title: data from exercise 8.55
> ### Aliases: ex08.55
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.55)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.57"
> 
> ### * ex08.57
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.57
> ### Title: data from exercise 8.57
> ### Aliases: ex08.57
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.57)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.66"
> 
> ### * ex08.66
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.66
> ### Title: data from exercise 8.66
> ### Aliases: ex08.66
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.66)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.70"
> 
> ### * ex08.70
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.70
> ### Title: data from exercise 8.70
> ### Aliases: ex08.70
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.70)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex08.80"
> 
> ### * ex08.80
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex08.80
> ### Title: data from exercise 8.80
> ### Aliases: ex08.80
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex08.80)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.07"
> 
> ### * ex09.07
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.07
> ### Title: data from exercise 9.7
> ### Aliases: ex09.07
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.07)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.12"
> 
> ### * ex09.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.12
> ### Title: data from exercise 9.12
> ### Aliases: ex09.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.12)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.16"
> 
> ### * ex09.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.16
> ### Title: data from exercise 9.16
> ### Aliases: ex09.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.23"
> 
> ### * ex09.23
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.23
> ### Title: data from exercise 9.23
> ### Aliases: ex09.23
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.23)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.25"
> 
> ### * ex09.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.25
> ### Title: data from exercise 9.25
> ### Aliases: ex09.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.27"
> 
> ### * ex09.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.27
> ### Title: data from exercise 9.27
> ### Aliases: ex09.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.27)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.28"
> 
> ### * ex09.28
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.28
> ### Title: data from exercise 9.28
> ### Aliases: ex09.28
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.28)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.29"
> 
> ### * ex09.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.29
> ### Title: data from exercise 9.29
> ### Aliases: ex09.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.30"
> 
> ### * ex09.30
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.30
> ### Title: data from exercise 9.30
> ### Aliases: ex09.30
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.30)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.31"
> 
> ### * ex09.31
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.31
> ### Title: data from exercise 9.31
> ### Aliases: ex09.31
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.31)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.32"
> 
> ### * ex09.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.32
> ### Title: data from exercise 9.32
> ### Aliases: ex09.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.33"
> 
> ### * ex09.33
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.33
> ### Title: data from exercise 9.33
> ### Aliases: ex09.33
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.33)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.36"
> 
> ### * ex09.36
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.36
> ### Title: data from exercise 9.36
> ### Aliases: ex09.36
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.36)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.37"
> 
> ### * ex09.37
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.37
> ### Title: data from exercise 9.37
> ### Aliases: ex09.37
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.37)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.38"
> 
> ### * ex09.38
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.38
> ### Title: data from exercise 9.38
> ### Aliases: ex09.38
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.38)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.39"
> 
> ### * ex09.39
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.39
> ### Title: data from exercise 9.39
> ### Aliases: ex09.39
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.39)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.40"
> 
> ### * ex09.40
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.40
> ### Title: data from exercise 9.40
> ### Aliases: ex09.40
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.40)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.41"
> 
> ### * ex09.41
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.41
> ### Title: data from exercise 9.41
> ### Aliases: ex09.41
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.41)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.43"
> 
> ### * ex09.43
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.43
> ### Title: data from exercise 9.43
> ### Aliases: ex09.43
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.43)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.44"
> 
> ### * ex09.44
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.44
> ### Title: data from exercise 9.44
> ### Aliases: ex09.44
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.44)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.63"
> 
> ### * ex09.63
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.63
> ### Title: data from exercise 9.63
> ### Aliases: ex09.63
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.63)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.65"
> 
> ### * ex09.65
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.65
> ### Title: data from exercise 9.65
> ### Aliases: ex09.65
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.65)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.66"
> 
> ### * ex09.66
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.66
> ### Title: data from exercise 9.66
> ### Aliases: ex09.66
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.66)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.68"
> 
> ### * ex09.68
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.68
> ### Title: data from exercise 9.68
> ### Aliases: ex09.68
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.68)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.70"
> 
> ### * ex09.70
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.70
> ### Title: data from exercise 9.70
> ### Aliases: ex09.70
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.70)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.72"
> 
> ### * ex09.72
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.72
> ### Title: data from exercise 9.72
> ### Aliases: ex09.72
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.72)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.76"
> 
> ### * ex09.76
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.76
> ### Title: data from exercise 9.76
> ### Aliases: ex09.76
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.76)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.77"
> 
> ### * ex09.77
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.77
> ### Title: data from exercise 9.77
> ### Aliases: ex09.77
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.77)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.78"
> 
> ### * ex09.78
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.78
> ### Title: data from exercise 9.78
> ### Aliases: ex09.78
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.78)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.79"
> 
> ### * ex09.79
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.79
> ### Title: data from exercise 9.79
> ### Aliases: ex09.79
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.79)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.82"
> 
> ### * ex09.82
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.82
> ### Title: data from exercise 9.82
> ### Aliases: ex09.82
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.82)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.86"
> 
> ### * ex09.86
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.86
> ### Title: data from exercise 9.86
> ### Aliases: ex09.86
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.86)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.88"
> 
> ### * ex09.88
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.88
> ### Title: data from exercise 9.88
> ### Aliases: ex09.88
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.88)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex09.92"
> 
> ### * ex09.92
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex09.92
> ### Title: data from exercise 9.92
> ### Aliases: ex09.92
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex09.92)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.06"
> 
> ### * ex10.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.06
> ### Title: data from exercise 10.6
> ### Aliases: ex10.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.08"
> 
> ### * ex10.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.08
> ### Title: data from exercise 10.8
> ### Aliases: ex10.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.09"
> 
> ### * ex10.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.09
> ### Title: data from exercise 10.9
> ### Aliases: ex10.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.18"
> 
> ### * ex10.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.18
> ### Title: data from exercise 10.18
> ### Aliases: ex10.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.22"
> 
> ### * ex10.22
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.22
> ### Title: data from exercise 10.22
> ### Aliases: ex10.22
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.22)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.26"
> 
> ### * ex10.26
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.26
> ### Title: data from exercise 10.26
> ### Aliases: ex10.26
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.26)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.27"
> 
> ### * ex10.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.27
> ### Title: data from exercise 10.27
> ### Aliases: ex10.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.27)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.32"
> 
> ### * ex10.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.32
> ### Title: data from exercise 10.32
> ### Aliases: ex10.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.36"
> 
> ### * ex10.36
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.36
> ### Title: data from exercise 10.36
> ### Aliases: ex10.36
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.36)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.37"
> 
> ### * ex10.37
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.37
> ### Title: data from exercise 10.37
> ### Aliases: ex10.37
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.37)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.41"
> 
> ### * ex10.41
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.41
> ### Title: data from exercise 10.41
> ### Aliases: ex10.41
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.41)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.42"
> 
> ### * ex10.42
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.42
> ### Title: data from exercise 10.42
> ### Aliases: ex10.42
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.42)
> str(ex10.42)
`data.frame':	19 obs. of  2 variables:
 $ cff  : num  26.8 27.9 23.7 25 26.3 24.8 25.7 24.5 26.4 24.2 ...
 $ color: Factor w/ 3 levels "Blue","Brown",..: 2 2 2 2 2 2 2 2 3 3 ...
> boxplot(cff ~ color, ex10.42, horizontal = TRUE, las = 1,
+   xlab = "Critical Flicker Frequency (Hz)")
> fm1 <- aov(cff ~ color, data = ex10.42)
> summary(fm1)
            Df Sum Sq Mean Sq F value  Pr(>F)  
color        2 22.997  11.499  4.8023 0.02325 *
Residuals   16 38.310   2.394                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "ex10.44"
> 
> ### * ex10.44
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex10.44
> ### Title: data from exercise 10.44
> ### Aliases: ex10.44
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex10.44)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.02"
> 
> ### * ex11.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.02
> ### Title: data from exercise 11.2
> ### Aliases: ex11.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.02)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.03"
> 
> ### * ex11.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.03
> ### Title: data from exercise 11.3
> ### Aliases: ex11.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.04"
> 
> ### * ex11.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.04
> ### Title: data from exercise 11.4
> ### Aliases: ex11.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.05"
> 
> ### * ex11.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.05
> ### Title: data from exercise 11.5
> ### Aliases: ex11.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.05)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.08"
> 
> ### * ex11.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.08
> ### Title: data from exercise 11.8
> ### Aliases: ex11.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.09"
> 
> ### * ex11.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.09
> ### Title: data from exercise 11.9
> ### Aliases: ex11.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.10"
> 
> ### * ex11.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.10
> ### Title: data from exercise 11.10
> ### Aliases: ex11.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.10)
> xtabs(strength ~ Batch + Method, data = ex11.10)
     Method
Batch    A    B    C
   1  30.7 33.7 30.5
   2  29.1 30.6 32.6
   3  30.0 32.2 30.5
   4  31.9 34.6 33.5
   5  30.5 33.0 32.4
   6  26.9 29.3 27.8
   7  28.2 28.4 30.7
   8  32.4 32.4 33.6
   9  26.6 29.5 29.2
   10 28.6 29.4 33.2
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.15"
> 
> ### * ex11.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.15
> ### Title: data from exercise 11.15
> ### Aliases: ex11.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.16"
> 
> ### * ex11.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.16
> ### Title: data from exercise 11.16
> ### Aliases: ex11.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.17"
> 
> ### * ex11.17
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.17
> ### Title: data from exercise 11.17
> ### Aliases: ex11.17
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.17)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.18"
> 
> ### * ex11.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.18
> ### Title: data from exercise 11.18
> ### Aliases: ex11.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.20"
> 
> ### * ex11.20
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.20
> ### Title: data from exercise 11.20
> ### Aliases: ex11.20
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.20)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.29"
> 
> ### * ex11.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.29
> ### Title: data from exercise 11.29
> ### Aliases: ex11.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.31"
> 
> ### * ex11.31
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.31
> ### Title: data from exercise 11.31
> ### Aliases: ex11.31
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.31)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.34"
> 
> ### * ex11.34
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.34
> ### Title: data from exercise 11.34
> ### Aliases: ex11.34
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.34)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.35"
> 
> ### * ex11.35
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.35
> ### Title: data from exercise 11.35
> ### Aliases: ex11.35
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.35)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.39"
> 
> ### * ex11.39
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.39
> ### Title: data from exercise 11.39
> ### Aliases: ex11.39
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.39)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.40"
> 
> ### * ex11.40
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.40
> ### Title: data from exercise 11.40
> ### Aliases: ex11.40
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.40)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.42"
> 
> ### * ex11.42
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.42
> ### Title: data from exercise 11.42
> ### Aliases: ex11.42
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.42)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.43"
> 
> ### * ex11.43
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.43
> ### Title: data from exercise 11.43
> ### Aliases: ex11.43
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.43)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.48"
> 
> ### * ex11.48
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.48
> ### Title: data from exercise 11.48
> ### Aliases: ex11.48
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.48)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.50"
> 
> ### * ex11.50
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.50
> ### Title: data from exercise 11.50
> ### Aliases: ex11.50
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.50)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.52"
> 
> ### * ex11.52
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.52
> ### Title: data from exercise 11.52
> ### Aliases: ex11.52
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.52)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.53"
> 
> ### * ex11.53
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.53
> ### Title: data from exercise 11.53
> ### Aliases: ex11.53
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.53)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.54"
> 
> ### * ex11.54
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.54
> ### Title: data from exercise 11.54
> ### Aliases: ex11.54
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.54)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.55"
> 
> ### * ex11.55
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.55
> ### Title: data from exercise 11.55
> ### Aliases: ex11.55
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.55)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.56"
> 
> ### * ex11.56
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.56
> ### Title: data from exercise 11.56
> ### Aliases: ex11.56
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.56)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.57"
> 
> ### * ex11.57
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.57
> ### Title: data from exercise 11.57
> ### Aliases: ex11.57
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.57)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.59"
> 
> ### * ex11.59
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.59
> ### Title: data from exercise 11.59
> ### Aliases: ex11.59
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.59)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex11.61"
> 
> ### * ex11.61
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex11.61
> ### Title: data from exercise 11.61
> ### Aliases: ex11.61
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex11.61)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.01"
> 
> ### * ex12.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.01
> ### Title: data from exercise 12.1
> ### Aliases: ex12.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.01)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.02"
> 
> ### * ex12.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.02
> ### Title: data from exercise 12.2
> ### Aliases: ex12.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.02)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.03"
> 
> ### * ex12.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.03
> ### Title: data from exercise 12.3
> ### Aliases: ex12.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.04"
> 
> ### * ex12.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.04
> ### Title: data from exercise 12.4
> ### Aliases: ex12.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.05"
> 
> ### * ex12.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.05
> ### Title: data from exercise 12.5
> ### Aliases: ex12.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.05)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.13"
> 
> ### * ex12.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.13
> ### Title: data from exercise 12.13
> ### Aliases: ex12.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.13)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.15"
> 
> ### * ex12.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.15
> ### Title: data from exercise 12.15
> ### Aliases: ex12.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.16"
> 
> ### * ex12.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.16
> ### Title: data from exercise 12.16
> ### Aliases: ex12.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.19"
> 
> ### * ex12.19
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.19
> ### Title: data from exercise 12.19
> ### Aliases: ex12.19
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.19)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.20"
> 
> ### * ex12.20
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.20
> ### Title: data from exercise 12.20
> ### Aliases: ex12.20
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.20)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.21"
> 
> ### * ex12.21
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.21
> ### Title: data from exercise 12.21
> ### Aliases: ex12.21
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.21)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.24"
> 
> ### * ex12.24
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.24
> ### Title: data from exercise 12.24
> ### Aliases: ex12.24
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.24)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.29"
> 
> ### * ex12.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.29
> ### Title: data from exercise 12.29
> ### Aliases: ex12.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.36"
> 
> ### * ex12.36
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.36
> ### Title: data from exercise 12.36
> ### Aliases: ex12.36
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.36)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.37"
> 
> ### * ex12.37
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.37
> ### Title: data from exercise 12.37
> ### Aliases: ex12.37
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.37)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.46"
> 
> ### * ex12.46
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.46
> ### Title: data from exercise 12.46
> ### Aliases: ex12.46
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.46)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.50"
> 
> ### * ex12.50
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.50
> ### Title: data from exercise 12.50
> ### Aliases: ex12.50
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.50)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.52"
> 
> ### * ex12.52
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.52
> ### Title: data from exercise 12.52
> ### Aliases: ex12.52
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.52)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.54"
> 
> ### * ex12.54
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.54
> ### Title: data from exercise 12.54
> ### Aliases: ex12.54
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.54)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.55"
> 
> ### * ex12.55
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.55
> ### Title: data from exercise 12.55
> ### Aliases: ex12.55
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.55)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.58"
> 
> ### * ex12.58
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.58
> ### Title: data from exercise 12.58
> ### Aliases: ex12.58
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.58)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.59"
> 
> ### * ex12.59
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.59
> ### Title: data from exercise 12.59
> ### Aliases: ex12.59
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.59)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.61"
> 
> ### * ex12.61
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.61
> ### Title: data from exercise 12.61
> ### Aliases: ex12.61
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.61)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.62"
> 
> ### * ex12.62
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.62
> ### Title: data from exercise 12.62
> ### Aliases: ex12.62
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.62)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.63"
> 
> ### * ex12.63
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.63
> ### Title: data from exercise 12.63
> ### Aliases: ex12.63
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.63)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.65"
> 
> ### * ex12.65
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.65
> ### Title: data from exercise 12.65
> ### Aliases: ex12.65
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.65)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.68"
> 
> ### * ex12.68
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.68
> ### Title: data from exercise 12.68
> ### Aliases: ex12.68
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.68)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.69"
> 
> ### * ex12.69
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.69
> ### Title: data from exercise 12.69
> ### Aliases: ex12.69
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.69)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.71"
> 
> ### * ex12.71
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.71
> ### Title: data from exercise 12.71
> ### Aliases: ex12.71
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.71)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.72"
> 
> ### * ex12.72
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.72
> ### Title: data from exercise 12.72
> ### Aliases: ex12.72
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.72)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.73"
> 
> ### * ex12.73
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.73
> ### Title: data from exercise 12.73
> ### Aliases: ex12.73
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.73)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.75"
> 
> ### * ex12.75
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.75
> ### Title: data from exercise 12.75
> ### Aliases: ex12.75
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.75)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.82"
> 
> ### * ex12.82
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.82
> ### Title: data from exercise 12.82
> ### Aliases: ex12.82
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.82)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.83"
> 
> ### * ex12.83
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.83
> ### Title: data from exercise 12.83
> ### Aliases: ex12.83
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.83)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex12.84"
> 
> ### * ex12.84
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex12.84
> ### Title: data from exercise 12.84
> ### Aliases: ex12.84
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex12.84)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.02"
> 
> ### * ex13.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.02
> ### Title: data from exercise 13.2
> ### Aliases: ex13.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.02)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.04"
> 
> ### * ex13.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.04
> ### Title: data from exercise 13.4
> ### Aliases: ex13.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.05"
> 
> ### * ex13.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.05
> ### Title: data from exercise 13.5
> ### Aliases: ex13.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.05)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.06"
> 
> ### * ex13.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.06
> ### Title: data from exercise 13.6
> ### Aliases: ex13.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.07"
> 
> ### * ex13.07
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.07
> ### Title: data from exercise 13.7
> ### Aliases: ex13.07
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.07)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.08"
> 
> ### * ex13.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.08
> ### Title: data from exercise 13.8
> ### Aliases: ex13.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.09"
> 
> ### * ex13.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.09
> ### Title: data from exercise 13.9
> ### Aliases: ex13.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.14"
> 
> ### * ex13.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.14
> ### Title: data from exercise 13.14
> ### Aliases: ex13.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.15"
> 
> ### * ex13.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.15
> ### Title: data from exercise 13.15
> ### Aliases: ex13.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.16"
> 
> ### * ex13.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.16
> ### Title: data from exercise 13.16
> ### Aliases: ex13.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.17"
> 
> ### * ex13.17
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.17
> ### Title: data from exercise 13.17
> ### Aliases: ex13.17
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.17)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.18"
> 
> ### * ex13.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.18
> ### Title: data from exercise 13.18
> ### Aliases: ex13.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.19"
> 
> ### * ex13.19
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.19
> ### Title: data from exercise 13.19
> ### Aliases: ex13.19
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.19)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.21"
> 
> ### * ex13.21
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.21
> ### Title: data from exercise 13.21
> ### Aliases: ex13.21
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.21)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.24"
> 
> ### * ex13.24
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.24
> ### Title: data from exercise 13.24
> ### Aliases: ex13.24
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.24)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.25"
> 
> ### * ex13.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.25
> ### Title: data from exercise 13.25
> ### Aliases: ex13.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.27"
> 
> ### * ex13.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.27
> ### Title: data from exercise 13.27
> ### Aliases: ex13.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.27)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.29"
> 
> ### * ex13.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.29
> ### Title: data from exercise 13.29
> ### Aliases: ex13.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.30"
> 
> ### * ex13.30
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.30
> ### Title: data from exercise 13.30
> ### Aliases: ex13.30
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.30)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.31"
> 
> ### * ex13.31
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.31
> ### Title: data from exercise 13.31
> ### Aliases: ex13.31
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.31)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.32"
> 
> ### * ex13.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.32
> ### Title: data from exercise 13.32
> ### Aliases: ex13.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.33"
> 
> ### * ex13.33
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.33
> ### Title: data from exercise 13.33
> ### Aliases: ex13.33
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.33)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.34"
> 
> ### * ex13.34
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.34
> ### Title: data from exercise 13.34
> ### Aliases: ex13.34
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.34)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.35"
> 
> ### * ex13.35
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.35
> ### Title: data from exercise 13.35
> ### Aliases: ex13.35
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.35)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.47"
> 
> ### * ex13.47
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.47
> ### Title: data from exercise 13.47
> ### Aliases: ex13.47
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.47)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.48"
> 
> ### * ex13.48
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.48
> ### Title: data from exercise 13.48
> ### Aliases: ex13.48
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.48)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.49"
> 
> ### * ex13.49
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.49
> ### Title: data from exercise 13.49
> ### Aliases: ex13.49
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.49)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.51"
> 
> ### * ex13.51
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.51
> ### Title: data from exercise 13.51
> ### Aliases: ex13.51
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.51)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.52"
> 
> ### * ex13.52
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.52
> ### Title: data from exercise 13.52
> ### Aliases: ex13.52
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.52)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.53"
> 
> ### * ex13.53
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.53
> ### Title: data from exercise 13.53
> ### Aliases: ex13.53
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.53)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.54"
> 
> ### * ex13.54
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.54
> ### Title: data from exercise 13.54
> ### Aliases: ex13.54
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.54)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.55"
> 
> ### * ex13.55
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.55
> ### Title: data from exercise 13.55
> ### Aliases: ex13.55
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.55)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.68"
> 
> ### * ex13.68
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.68
> ### Title: data from exercise 13.68
> ### Aliases: ex13.68
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.68)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.69"
> 
> ### * ex13.69
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.69
> ### Title: data from exercise 13.69
> ### Aliases: ex13.69
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.69)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.70"
> 
> ### * ex13.70
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.70
> ### Title: data from exercise 13.70
> ### Aliases: ex13.70
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.70)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.71"
> 
> ### * ex13.71
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.71
> ### Title: data from exercise 13.71
> ### Aliases: ex13.71
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.71)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.73"
> 
> ### * ex13.73
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.73
> ### Title: data from exercise 13.73
> ### Aliases: ex13.73
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.73)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.74"
> 
> ### * ex13.74
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.74
> ### Title: data from exercise 13.74
> ### Aliases: ex13.74
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.74)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.75"
> 
> ### * ex13.75
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.75
> ### Title: data from exercise 13.75
> ### Aliases: ex13.75
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.75)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex13.76"
> 
> ### * ex13.76
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex13.76
> ### Title: data from exercise 13.76
> ### Aliases: ex13.76
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex13.76)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.09"
> 
> ### * ex14.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.09
> ### Title: data from exercise 14.9
> ### Aliases: ex14.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.11"
> 
> ### * ex14.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.11
> ### Title: data from exercise 14.11
> ### Aliases: ex14.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.11)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.12"
> 
> ### * ex14.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.12
> ### Title: data from exercise 14.12
> ### Aliases: ex14.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.12)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.13"
> 
> ### * ex14.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.13
> ### Title: data from exercise 14.13
> ### Aliases: ex14.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.13)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.14"
> 
> ### * ex14.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.14
> ### Title: data from exercise 14.14
> ### Aliases: ex14.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.15"
> 
> ### * ex14.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.15
> ### Title: data from exercise 14.15
> ### Aliases: ex14.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.16"
> 
> ### * ex14.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.16
> ### Title: data from exercise 14.16
> ### Aliases: ex14.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.17"
> 
> ### * ex14.17
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.17
> ### Title: data from exercise 14.17
> ### Aliases: ex14.17
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.17)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.18"
> 
> ### * ex14.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.18
> ### Title: data from exercise 14.18
> ### Aliases: ex14.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.20"
> 
> ### * ex14.20
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.20
> ### Title: data from exercise 14.20
> ### Aliases: ex14.20
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.20)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.21"
> 
> ### * ex14.21
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.21
> ### Title: data from exercise 14.21
> ### Aliases: ex14.21
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.21)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.22"
> 
> ### * ex14.22
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.22
> ### Title: data from exercise 14.22
> ### Aliases: ex14.22
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.22)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.23"
> 
> ### * ex14.23
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.23
> ### Title: data from exercise 14.23
> ### Aliases: ex14.23
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.23)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.26"
> 
> ### * ex14.26
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.26
> ### Title: data from exercise 14.26
> ### Aliases: ex14.26
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.26)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.27"
> 
> ### * ex14.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.27
> ### Title: data from exercise 14.27
> ### Aliases: ex14.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.27)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.28"
> 
> ### * ex14.28
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.28
> ### Title: data from exercise 14.28
> ### Aliases: ex14.28
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.28)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.29"
> 
> ### * ex14.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.29
> ### Title: data from exercise 14.29
> ### Aliases: ex14.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.30"
> 
> ### * ex14.30
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.30
> ### Title: data from exercise 14.30
> ### Aliases: ex14.30
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.30)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.31"
> 
> ### * ex14.31
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.31
> ### Title: data from exercise 14.31
> ### Aliases: ex14.31
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.31)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.32"
> 
> ### * ex14.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.32
> ### Title: data from exercise 14.32
> ### Aliases: ex14.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.38"
> 
> ### * ex14.38
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.38
> ### Title: data from exercise 14.38
> ### Aliases: ex14.38
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.38)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.40"
> 
> ### * ex14.40
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.40
> ### Title: data from exercise 14.40
> ### Aliases: ex14.40
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.40)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.41"
> 
> ### * ex14.41
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.41
> ### Title: data from exercise 14.41
> ### Aliases: ex14.41
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.41)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.42"
> 
> ### * ex14.42
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.42
> ### Title: data from exercise 14.42
> ### Aliases: ex14.42
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.42)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex14.44"
> 
> ### * ex14.44
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex14.44
> ### Title: data from exercise 14.44
> ### Aliases: ex14.44
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex14.44)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.01"
> 
> ### * ex15.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.01
> ### Title: data from exercise 15.1
> ### Aliases: ex15.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.01)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.03"
> 
> ### * ex15.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.03
> ### Title: data from exercise 15.3
> ### Aliases: ex15.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.04"
> 
> ### * ex15.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.04
> ### Title: data from exercise 15.4
> ### Aliases: ex15.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.05"
> 
> ### * ex15.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.05
> ### Title: data from exercise 15.5
> ### Aliases: ex15.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.05)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.08"
> 
> ### * ex15.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.08
> ### Title: data from exercise 15.8
> ### Aliases: ex15.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.10"
> 
> ### * ex15.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.10
> ### Title: data from exercise 15.10
> ### Aliases: ex15.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.11"
> 
> ### * ex15.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.11
> ### Title: data from exercise 15.11
> ### Aliases: ex15.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.11)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.12"
> 
> ### * ex15.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.12
> ### Title: data from exercise 15.12
> ### Aliases: ex15.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.12)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.13"
> 
> ### * ex15.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.13
> ### Title: data from exercise 15.13
> ### Aliases: ex15.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.13)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.14"
> 
> ### * ex15.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.14
> ### Title: data from exercise 15.14
> ### Aliases: ex15.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.15"
> 
> ### * ex15.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.15
> ### Title: data from exercise 15.15
> ### Aliases: ex15.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.23"
> 
> ### * ex15.23
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.23
> ### Title: data from exercise 15.23
> ### Aliases: ex15.23
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.23)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.24"
> 
> ### * ex15.24
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.24
> ### Title: data from exercise 15.24
> ### Aliases: ex15.24
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.24)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.25"
> 
> ### * ex15.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.25
> ### Title: data from exercise 15.25
> ### Aliases: ex15.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.26"
> 
> ### * ex15.26
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.26
> ### Title: data from exercise 15.26
> ### Aliases: ex15.26
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.26)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.27"
> 
> ### * ex15.27
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.27
> ### Title: data from exercise 15.27
> ### Aliases: ex15.27
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.27)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.28"
> 
> ### * ex15.28
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.28
> ### Title: data from exercise 15.28
> ### Aliases: ex15.28
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.28)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.29"
> 
> ### * ex15.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.29
> ### Title: data from exercise 15.29
> ### Aliases: ex15.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.29)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.30"
> 
> ### * ex15.30
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.30
> ### Title: data from exercise 15.30
> ### Aliases: ex15.30
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.30)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.32"
> 
> ### * ex15.32
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.32
> ### Title: data from exercise 15.32
> ### Aliases: ex15.32
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.32)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.33"
> 
> ### * ex15.33
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.33
> ### Title: data from exercise 15.33
> ### Aliases: ex15.33
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.33)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex15.35"
> 
> ### * ex15.35
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex15.35
> ### Title: data from exercise 15.35
> ### Aliases: ex15.35
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex15.35)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.06"
> 
> ### * ex16.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.06
> ### Title: data from exercise 16.6
> ### Aliases: ex16.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.09"
> 
> ### * ex16.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.09
> ### Title: data from exercise 16.9
> ### Aliases: ex16.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.14"
> 
> ### * ex16.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.14
> ### Title: data from exercise 16.14
> ### Aliases: ex16.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.25"
> 
> ### * ex16.25
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.25
> ### Title: data from exercise 16.25
> ### Aliases: ex16.25
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.25)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.41"
> 
> ### * ex16.41
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.41
> ### Title: data from exercise 16.41
> ### Aliases: ex16.41
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.41)
> 
> 
> 
> cleanEx(); ..nameEx <- "ex16.43"
> 
> ### * ex16.43
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: ex16.43
> ### Title: data from exercise 16.43
> ### Aliases: ex16.43
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(ex16.43)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.01"
> 
> ### * xmp01.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.01
> ### Title: data from Example 1.1
> ### Aliases: xmp01.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.01)
> attach(xmp01.01)
> summary(temp)             # summary statistics
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  31.00   59.50   67.50   65.86   75.00   84.00 
> stem(temp)

  The decimal point is 1 digit(s) to the right of the |

  3 | 1
  4 | 059
  5 | 23788
  6 | 01136777789
  7 | 000023556689
  8 | 0134

> hist(temp, xlab = "Temperature (deg. F)")
> rug(temp)
> hist(temp, xlab = "Temperature (deg. F)",
+      prob = TRUE, col = "lightgray")
> lines(density(temp), col = "blue")
> rug(temp)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.02"
> 
> ### * xmp01.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.02
> ### Title: data from Example 1.2
> ### Aliases: xmp01.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.02)
> attach(xmp01.02)
> hist(strength, xlab = "Flexural strength (MPa)",
+      col = "lightgray")
> rug(strength)
> summary(strength)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  5.900   7.000   7.700   8.141   8.850  11.800 
> boxplot(strength, col = "lightgray", notch = TRUE,
+         ylab = "Flexural strength (MPa)",
+         main = "Boxplot of strength",
+         sub =
+ "Notches show a 95% confidence interval on the median strength")
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.05"
> 
> ### * xmp01.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.05
> ### Title: data from Example 1.5
> ### Aliases: xmp01.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.05)
> str(xmp01.05)
`data.frame':	140 obs. of  1 variable:
 $ bingePct: int  4 11 13 14 15 16 17 18 18 18 ...
> stem(xmp01.05$bingePct)

  The decimal point is 1 digit(s) to the right of the |

  0 | 4
  0 | 
  1 | 134
  1 | 5678889
  2 | 12234
  2 | 56666777889999
  3 | 0112233344
  3 | 555666677777888899999
  4 | 11122222334444
  4 | 5566666677788888999
  5 | 001112222334
  5 | 55666667777888899
  6 | 011112444
  6 | 55666778

> stem(xmp01.05$bingePct, scale = 0.5)  # compare to Figure 1.4, p. 12

  The decimal point is 1 digit(s) to the right of the |

  0 | 4
  1 | 1345678889
  2 | 1223456666777889999
  3 | 0112233344555666677777888899999
  4 | 111222223344445566666677788888999
  5 | 00111222233455666667777888899
  6 | 01111244455666778

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.06"
> 
> ### * xmp01.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.06
> ### Title: data from Example 1.6
> ### Aliases: xmp01.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.06)
> attach(xmp01.06)
> summary(yardage)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
   6433    6674    6872    6845    7043    7280 
> stem(yardage)

  The decimal point is 2 digit(s) to the right of the |

  64 | 3467
  65 | 1338
  66 | 119
  67 | 015779
  68 | 05779
  69 | 0034
  70 | 112455
  71 | 113777
  72 | 18

> hist(yardage, col = "lightgray",
+      xlab = "Golf course yardages")
> rug(yardage)
> qqnorm(yardage, las = 1, ylab = "Golf course yardages")
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.10"
> 
> ### * xmp01.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.10
> ### Title: data from Example 1.10
> ### Aliases: xmp01.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.10)
> attach(xmp01.10)
> hist(consump, col = "lightgray",
+      xlab = "Adjusted energy consumption")
> rug(consump)
> hist(consump, col = "lightgray",
+      xlab = "Adjusted energy consumption",
+      prob = TRUE)
> lines(density(consump), col = "blue")
> rug(consump)
> summary(consump)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  2.970   7.948   9.835  10.040  12.050  18.260 
> # Make a histogram like Fig 1.9, p. 19
> hist(consump, breaks = 1 + 2*(0:9),
+      xlab = "BTUN", prob = TRUE, col = "lightgray")
> rug(consump)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.11"
> 
> ### * xmp01.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.11
> ### Title: data from Example 1.11
> ### Aliases: xmp01.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.11)
> attach(xmp01.11)
> hist(strength, xlab = "Bond strength", col = "lightgray")
> rug(strength)
> hist(log(strength), xlab = "log(bond strength)", col = "lightgray")
> rug(log(strength))
> ## Create a histogram like Fig 1.11, page 20
> hist(strength, breaks = c(2,4,6,8,12,20,30), prob = TRUE,
+      col = "lightgray", xlab = "Bond strength")
> rug(strength)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.13"
> 
> ### * xmp01.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.13
> ### Title: data from Example 1.13
> ### Aliases: xmp01.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.13)
> str(xmp01.13)
`data.frame':	21 obs. of  1 variable:
 $ crackLength: num  16.1 9.6 24.9 20.4 12.7 21.2 30.2 25.8 18.5 10.3 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.14"
> 
> ### * xmp01.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.14
> ### Title: data from Example 1.14
> ### Aliases: xmp01.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.14)
> str(xmp01.14)
`data.frame':	12 obs. of  1 variable:
 $ concentration: num  15.2 9.3 7.6 11.9 10.4 9.7 20.4 9.4 11.5 16.2 ...
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.15"
> 
> ### * xmp01.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.15
> ### Title: data from Example 1.15
> ### Aliases: xmp01.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.15)
> attach(xmp01.15)
> summary(lifetime)              # produces mean, median, etc.
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  612.0   894.3  1010.0   965.0  1086.0  1201.0 
> mean(lifetime, trim = 0.1)     # 10% trimmed mean
[1] 979.125
> mean(lifetime, trim = 0.2)     # 20% trimmed mean
[1] 999.9167
> dotchart(lifetime)
> hist(lifetime)                 # display a histogram
> rug(lifetime)                  # add the data
> abline(v = median(lifetime), col = 2, lty = 2)
> abline(v = mean(lifetime), col = 3, lty = 2)
> abline(v = mean(lifetime, trim = 0.1), col = 4, lty = 2)
> abline(v = mean(lifetime, trim = 0.2), col = 5, lty = 2)
> legend(600, 6,
+   c("median", "mean", "10% trimmed mean", "20% trimmed mean"),
+   col = 2:5, lty = 2)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.16"
> 
> ### * xmp01.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.16
> ### Title: data from Example 1.16
> ### Aliases: xmp01.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.18"
> 
> ### * xmp01.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.18
> ### Title: data from Example 1.18
> ### Aliases: xmp01.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp01.19"
> 
> ### * xmp01.19
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp01.19
> ### Title: data from Example 1.19
> ### Aliases: xmp01.19
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp01.19)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp04.28"
> 
> ### * xmp04.28
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp04.28
> ### Title: data from Example 4.28
> ### Aliases: xmp04.28
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp04.28)
> attach(xmp04.28)
> qqnorm(meas.err)   # compare to Figure 4.31, p. 188
> qqline(meas.err)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp04.29"
> 
> ### * xmp04.29
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp04.29
> ### Title: data from Example 4.29
> ### Aliases: xmp04.29
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp04.29)
> attach(xmp04.29)
> ## compare to Figure 4.33, page 190
> qqp <- qqnorm(Voltage)
> qqline(Voltage)
> detach()
> ## compare quantiles with those given in book
> cbind(qqp, xmp04.29)
             x     y Voltage z.percentile
1  -1.95996398 24.46   24.46        -1.96
2  -1.43953147 25.61   25.61        -1.44
3  -1.15034938 26.25   26.25        -1.15
4  -0.93458929 26.42   26.42        -0.93
5  -0.75541503 26.66   26.66        -0.76
6  -0.59776013 27.15   27.15        -0.60
7  -0.45376219 27.31   27.31        -0.45
8  -0.31863936 27.54   27.54        -0.32
9  -0.18911843 27.74   27.74        -0.19
10 -0.06270678 27.94   27.94        -0.06
11  0.06270678 27.98   27.98         0.06
12  0.18911843 28.04   28.04         0.19
13  0.31863936 28.28   28.28         0.32
14  0.45376219 28.49   28.49         0.45
15  0.59776013 28.50   28.50         0.60
16  0.75541503 28.87   28.87         0.76
17  0.93458929 29.11   29.11         0.93
18  1.15034938 29.13   29.13         1.15
19  1.43953147 29.50   29.50         1.44
20  1.95996398 30.88   30.88         1.96
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp04.30"
> 
> ### * xmp04.30
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp04.30
> ### Title: data from Example 4.30
> ### Aliases: xmp04.30
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp04.30)
> attach(xmp04.30)
> ## Try normal probability plot first
> qqnorm(lifetime, ylab = "Lifetime (hr)")
> qqline(lifetime)
> ## Weibull probability plot, compare Figure 4.36, p. 194
> plot(log(-log(1 - seq(0.05, 0.95, 0.1))),
+      log(sort(lifetime)), xlab = "Theoretical Quantiles",
+      ylab = "log(Lifetime) (log(hr))",
+      main = "Weibull Q-Q Plot", las = 1)
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp06.02"
> 
> ### * xmp06.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp06.02
> ### Title: data from Example 6.2
> ### Aliases: xmp06.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp06.02)
> summary(xmp06.02)                # gives mean, median, etc.
    Voltage     
 Min.   :24.46  
 1st Qu.:27.03  
 Median :27.96  
 Mean   :27.79  
 3rd Qu.:28.59  
 Max.   :30.88  
> attach(xmp06.02)
> mean(range(Voltage))             # average of the extremes
[1] 27.67
> mean(Voltage, trim = 0.1)        # trimmed mean
[1] 27.83812
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp06.03"
> 
> ### * xmp06.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp06.03
> ### Title: data from Example 6.3
> ### Aliases: xmp06.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp06.03)
> attach(xmp06.03)
> stem(Strength)

  The decimal point is at the |

  43 | 1
  43 | 89
  44 | 022
  44 | 67

> var(Strength)  # usual (unbiased) estimate of sigma^2
[1] 0.25125
> ## alternative estimate of sigma^2 with n in denominator
> sum((Strength - mean(Strength))^2)/length(Strength)
[1] 0.2198438
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp06.12"
> 
> ### * xmp06.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp06.12
> ### Title: data from Example 6.12
> ### Aliases: xmp06.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp06.12)
> attach(xmp06.12)
> gamma.MoM <- function(x) {
+    ## calculate method of moments estimates for gamma distribution
+    xbar <- mean(x)
+    mnSqDev <- mean((x - xbar)^2)
+    c(alpha = xbar^2/mnSqDev, beta = mnSqDev/xbar)
+ }
> ## method of moments estimates
> print(surv.MoM <- gamma.MoM(Survival))
   alpha     beta 
10.57725 10.72585 
> ## evaluating the negative log-likelihood
> gammaLlik <- function(x) {
+    ## argument x is a vector of shape (alpha) and scale (beta)
+    -sum(dgamma(Survival, shape = x[1], scale = x[2], log = TRUE))
+ }
> ## maximum likelihood estimates - use MoM estimates as starting value
> MLE <- optim(par = surv.MoM, gammaLlik)
> print(MLE)
$par
    alpha      beta 
 8.792303 12.904491 

$value
[1] 100.4799

$counts
function gradient 
      45       NA 

$convergence
[1] 0

$message
NULL

> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp06.13"
> 
> ### * xmp06.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp06.13
> ### Title: data from Example 6.13
> ### Aliases: xmp06.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp06.13)
> attach(xmp06.13)
> table(goals)      # compare to frequency table on p. 267
goals
 0  1  2  3  4  5  6  7  8  9 10 
29 71 82 89 65 45 24  7  4  1  3 
> hist(goals, breaks = 0:12 - 0.5, las = 1, col = "lightgray")
> negBinom.MoM <- function(x) {
+    ## method of moments estimates for negative binomial distribution
+    xbar <- mean(x)
+    mnSqDev <- mean((x - xbar)^2)
+    c(p = xbar/mnSqDev, r = xbar^2/(mnSqDev - xbar))
+ }
> print(goals.MoM <- negBinom.MoM(goals))
       p        r 
 0.84367 16.07454 
> ## MLE's
> optim(goals.MoM, function(x)
+   -sum(dnbinom(goals, p = x[1], size = x[2], log = TRUE)))
Warning in dnbinom(x, size, prob, log) : NaNs produced
Warning in dnbinom(x, size, prob, log) : NaNs produced
Warning in dnbinom(x, size, prob, log) : NaNs produced
Warning in dnbinom(x, size, prob, log) : NaNs produced
Warning in dnbinom(x, size, prob, log) : NaNs produced
Warning in dnbinom(x, size, prob, log) : NaNs produced
$par
         p          r 
 0.8410075 15.7549316 

$value
[1] 838.8255

$counts
function gradient 
      59       NA 

$convergence
[1] 0

$message
NULL

> ## would have been better to use a transformation of p
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp07.06"
> 
> ### * xmp07.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp07.06
> ### Title: data from Example 7.6
> ### Aliases: xmp07.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp07.06)
> boxplot(xmp07.06,
+   main = "AC Breakdown Voltage (kV)")
> # t.test gives a 95% confidence interval on the mean
> t.test(xmp07.06)   

	One Sample t-test

data:  xmp07.06 
t = 72.4631, df = 47, p-value < 2.2e-16
alternative hypothesis: true mean is not equal to 0 
95 percent confidence interval:
 53.18950 56.22716 
sample estimates:
mean of x 
 54.70833 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp07.11"
> 
> ### * xmp07.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp07.11
> ### Title: data from Example 7.11
> ### Aliases: xmp07.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp07.11)
> boxplot(xmp07.11)
> with(xmp07.11, qqnorm(Elasticity))
> with(xmp07.11, qqline(Elasticity))
> with(xmp07.11, t.test(Elasticity))

	One Sample t-test

data:  Elasticity 
t = 28.2778, df = 15, p-value = 1.989e-14
alternative hypothesis: true mean is not equal to 0 
95 percent confidence interval:
 13437.11 15627.89 
sample estimates:
mean of x 
  14532.5 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp07.15"
> 
> ### * xmp07.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp07.15
> ### Title: data from Example 7.15
> ### Aliases: xmp07.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp07.15)
> boxplot(xmp07.15, main = "Breakdown voltage")
> with(xmp07.15, qqnorm(voltage, main = "Breakdown voltage"))
> with(xmp07.15, qqline(voltage))
> attach(xmp07.15)
> var(voltage) * (length(voltage) - 1)/
+   qchisq(c(0.975, 0.025), df = length(voltage) - 1)
[1]  76171.31 318079.76
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp08.08"
> 
> ### * xmp08.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp08.08
> ### Title: data from Example 8.8
> ### Aliases: xmp08.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp08.08)
> boxplot(xmp08.08, main = "DCP readings")
> attach(xmp08.08)
> hist(DCP, breaks = 8, prob = TRUE)
> rug(DCP)
> lines(density(DCP), col = "blue")
> t.test(DCP, alt = "less", mu = 30)

	One Sample t-test

data:  DCP 
t = -0.7282, df = 51, p-value = 0.2349
alternative hypothesis: true mean is less than 30 
95 percent confidence interval:
     -Inf 31.61088 
sample estimates:
mean of x 
 28.76154 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp08.09"
> 
> ### * xmp08.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp08.09
> ### Title: data from Example 8.9
> ### Aliases: xmp08.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp08.09)
> str(xmp08.09)
`data.frame':	5 obs. of  1 variable:
 $ MAWL: num  25.8 36.6 26.3 21.8 27.2
> with(xmp08.09, t.test(MAWL, mu = 25, alt = "greater"))

	One Sample t-test

data:  MAWL 
t = 1.0382, df = 4, p-value = 0.1789
alternative hypothesis: true mean is greater than 25 
95 percent confidence interval:
 22.32433      Inf 
sample estimates:
mean of x 
    27.54 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp09.04"
> 
> ### * xmp09.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp09.04
> ### Title: data from Example 9.4
> ### Aliases: xmp09.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp09.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp09.06"
> 
> ### * xmp09.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp09.06
> ### Title: data from Example 9.6
> ### Aliases: xmp09.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp09.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp09.08"
> 
> ### * xmp09.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp09.08
> ### Title: data from Example 9.8
> ### Aliases: xmp09.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp09.08)
> boxplot(xmp09.08, main = "Boxplot of data from Example 9.8")
> attach(xmp09.08)
> boxplot(bottom-surface, main = "Boxplot of differences from Example 9.8")
> t.test(bottom, surface, alt = "greater", paired = TRUE)

	Paired t-test

data:  bottom and surface 
t = 3.6998, df = 5, p-value = 0.007001
alternative hypothesis: true difference in means is greater than 0 
95 percent confidence interval:
 0.04174206        Inf 
sample estimates:
mean of the differences 
             0.09166667 

> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp09.09"
> 
> ### * xmp09.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp09.09
> ### Title: data from Example 9.9
> ### Aliases: xmp09.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp09.09)
> boxplot(xmp09.09[, c("Before", "After")],
+    main = "Data from Example 9.9")
> attach(xmp09.09)
> boxplot(Difference, main = "Differences in Example 9.9")
> qqnorm(Difference,
+    main = "Normal probability plot (compare Figure 9.5, p. 377)")
> t.test(Difference)

	One Sample t-test

data:  Difference 
t = 3.2791, df = 15, p-value = 0.005072
alternative hypothesis: true mean is not equal to 0 
95 percent confidence interval:
  2.362371 11.137629 
sample estimates:
mean of x 
     6.75 

> t.test(Before, After, paired = TRUE)   # same test

	Paired t-test

data:  Before and After 
t = 3.2791, df = 15, p-value = 0.005072
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
  2.362371 11.137629 
sample estimates:
mean of the differences 
                   6.75 

> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp09.10"
> 
> ### * xmp09.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp09.10
> ### Title: data from Example 9.10
> ### Aliases: xmp09.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp09.10)
> boxplot(xmp09.10[, c("t1.min", "t4.wks")],
+      main = "Data from Example 9.10")
> attach(xmp09.10)
> ## compare to Figure 9.7, page 379
> qqnorm(Difference, main = "Differences from Example 9.10",
+        ylab = "Difference in modulus of elasticity")
> qqline(Difference)
> t.test(Difference, conf = 0.99)

	One Sample t-test

data:  Difference 
t = 20.7266, df = 15, p-value = 1.88e-12
alternative hypothesis: true mean is not equal to 0 
99 percent confidence interval:
 2260.917 3010.333 
sample estimates:
mean of x 
 2635.625 

> t.test(t1.min, t4.wks, paired = TRUE, conf = 0.99) # same thing

	Paired t-test

data:  t1.min and t4.wks 
t = 20.7266, df = 15, p-value = 1.88e-12
alternative hypothesis: true difference in means is not equal to 0 
99 percent confidence interval:
 2260.917 3010.333 
sample estimates:
mean of the differences 
               2635.625 

> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp10.01"
> 
> ### * xmp10.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp10.01
> ### Title: data from Example 10.1
> ### Aliases: xmp10.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp10.01)
> boxplot(strength ~ type, data = xmp10.01, 
+     main = "Data from Example (compare Figure 10.1, p. 405)")
> fm1 <- lm( strength ~ type, data = xmp10.01 )  # fit anova model
> qqnorm(resid(fm1), main = "Compare to Figure 10.2, p. 407")
> anova(fm1)     # compare results in Example 10.2, p. 409
Analysis of Variance Table

Response: strength
          Df Sum Sq Mean Sq F value    Pr(>F)    
type       3 127375   42458  25.094 5.525e-07 ***
Residuals 20  33839    1692                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp10.03"
> 
> ### * xmp10.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp10.03
> ### Title: data from Example 10.3
> ### Aliases: xmp10.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp10.03)
> xmp10.03$Mixture <- factor(xmp10.03$Mixture)
> plot(Soiling ~ Mixture, data = xmp10.03, col = "lightgray",
+      main = "Data from Example 10.3")
> summary(xmp10.03)          # check ranges and balance
    Soiling       Mixture
 Min.   :0.5600   1:5    
 1st Qu.:0.7500   2:5    
 Median :0.9100   3:5    
 Mean   :0.8833          
 3rd Qu.:1.0300          
 Max.   :1.1200          
> fm1 <- lm(Soiling ~ Mixture, data = xmp10.03)
> anova(fm1)                 # compare to table shown on p. 412
Analysis of Variance Table

Response: Soiling
          Df  Sum Sq Mean Sq F value Pr(>F)
Mixture    2 0.06085 0.03043  0.9866 0.4012
Residuals 12 0.37008 0.03084               
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp10.05"
> 
> ### * xmp10.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp10.05
> ### Title: data from Example 10.5
> ### Aliases: xmp10.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp10.05)
> plot(REMtime ~ ethanol, data = xmp10.05,
+      xlab = "Ethanol concentration administered (g/kg)",
+      ylab = "Amount of REM sleep during a 24 hour period")
> fm1 <- lm(REMtime ~ factor(ethanol), data = xmp10.05)
> anova(fm1)      # compare with Table 10.4, p. 417
Analysis of Variance Table

Response: REMtime
                Df Sum Sq Mean Sq F value    Pr(>F)    
factor(ethanol)  3 5882.4  1960.8  21.092 8.325e-06 ***
Residuals       16 1487.4    93.0                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> summary(fm1)    # differences with baseline (0 g/kg)

Call:
lm(formula = REMtime ~ factor(ethanol), data = xmp10.05)

Residuals:
    Min      1Q  Median      3Q     Max 
-11.440  -7.660  -2.390   8.615  12.540 

Coefficients:
                 Estimate Std. Error t value Pr(>|t|)    
(Intercept)        79.280      4.312  18.386 3.49e-12 ***
factor(ethanol)1  -17.740      6.098  -2.909   0.0102 *  
factor(ethanol)2  -31.360      6.098  -5.143 9.83e-05 ***
factor(ethanol)4  -46.520      6.098  -7.629 1.02e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 9.642 on 16 degrees of freedom
Multiple R-Squared: 0.7982,	Adjusted R-squared: 0.7603 
F-statistic: 21.09 on 3 and 16 DF,  p-value: 8.325e-06 

> ##  more appropriate to use an ordered factor
> fm2 <- lm(REMtime ~ ordered(ethanol), data = xmp10.05)
> anova(fm2)      # same as above
Analysis of Variance Table

Response: REMtime
                 Df Sum Sq Mean Sq F value    Pr(>F)    
ordered(ethanol)  3 5882.4  1960.8  21.092 8.325e-06 ***
Residuals        16 1487.4    93.0                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> summary(fm2)    # polynomial contrasts

Call:
lm(formula = REMtime ~ ordered(ethanol), data = xmp10.05)

Residuals:
    Min      1Q  Median      3Q     Max 
-11.440  -7.660  -2.390   8.615  12.540 

Coefficients:
                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)          55.375      2.156  25.685 1.96e-14 ***
ordered(ethanol).L  -34.252      4.312  -7.944 6.08e-07 ***
ordered(ethanol).Q    1.290      4.312   0.299    0.769    
ordered(ethanol).C   -1.266      4.312  -0.294    0.773    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 9.642 on 16 degrees of freedom
Multiple R-Squared: 0.7982,	Adjusted R-squared: 0.7603 
F-statistic: 21.09 on 3 and 16 DF,  p-value: 8.325e-06 

> ##  best model uses square root of ethanol concentration
> plot(REMtime ~ sqrt(ethanol), data = xmp10.05,
+      xlab = expression(sqrt(
+         plain("Ethanol concentration administered (g/kg)"))),
+      ylab = "Amount of REM sleep during a 24 hour period")
> fm3 <- lm(REMtime ~ sqrt(ethanol), data = xmp10.05)
> summary(fm3)

Call:
lm(formula = REMtime ~ sqrt(ethanol), data = xmp10.05)

Residuals:
    Min      1Q  Median      3Q     Max 
-13.093  -7.903  -3.359   8.700  13.712 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)      81.09       3.81  21.283 3.28e-14 ***
sqrt(ethanol)   -23.30       2.88  -8.091 2.08e-07 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 9.397 on 18 degrees of freedom
Multiple R-Squared: 0.7843,	Adjusted R-squared: 0.7724 
F-statistic: 65.47 on 1 and 18 DF,  p-value: 2.084e-07 

> abline(fm3)
> anova(fm3, fm1)   # lack of fit test
Analysis of Variance Table

Model 1: REMtime ~ sqrt(ethanol)
Model 2: REMtime ~ factor(ethanol)
  Res.Df    RSS Df Sum of Sq      F Pr(>F)
1     18 1589.3                           
2     16 1487.4  2     101.9 0.5481 0.5885
> opar <- par(mfrow = c(2,2))
> plot(fm3, main = "Continuous fit to data in Example 10.5")
> par(opar)
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp10.08"
> 
> ### * xmp10.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp10.08
> ### Title: data from Example 10.8
> ### Aliases: xmp10.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp10.08)
> str(xmp10.08)
`data.frame':	22 obs. of  2 variables:
 $ elastic: num  45.5 45.3 45.4 44.4 44.6 43.9 44.6 44 44.2 43.9 ...
 $ type   : Factor w/ 3 levels "Die","Permanent",..: 2 2 2 2 2 2 2 2 1 1 ...
> fm1 <- aov(elastic ~ type, data = xmp10.08)
> anova(fm1)
Analysis of Variance Table

Response: elastic
          Df Sum Sq Mean Sq F value    Pr(>F)    
type       2 7.9301  3.9650  12.565 0.0003336 ***
Residuals 19 5.9958  0.3156                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp10.10"
> 
> ### * xmp10.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp10.10
> ### Title: data from Example 10.10
> ### Aliases: xmp10.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp10.10)
> xmp10.10$Rail <- factor(xmp10.10$Rail)
> boxplot(travel ~ Rail, xmp10.10, col = "lightgray",
+   xlab = "Rail", ylab = "Zero-force travel time (microsec)",
+   main = "Travel times in rails, from example 10.10")
> fm1 <- lm(travel ~ Rail, data = xmp10.10)
> anova(fm1)
Analysis of Variance Table

Response: travel
          Df Sum Sq Mean Sq F value    Pr(>F)    
Rail       5 9310.5  1862.1  115.18 1.033e-09 ***
Residuals 12  194.0    16.2                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.01"
> 
> ### * xmp11.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.01
> ### Title: data from Example 11.1
> ### Aliases: xmp11.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.01)
> xmp11.01$brand <- factor(xmp11.01$brand)
> xmp11.01$treatment <- factor(xmp11.01$treatment)
> plot(strength ~ treatment, data = xmp11.01, col = "lightgray",
+      main = "Interaction plot for Example 11.01",
+      xlab = "Washing treatment")
> lines(strength ~ as.integer(treatment), data = xmp11.01,
+       subset = brand == 1, col = 4, type = "b")
> lines(strength ~ as.integer(treatment), data = xmp11.01,
+       subset = brand == 2, col = 2, type = "b")
> lines(strength ~ as.integer(treatment), data = xmp11.01,
+       subset = brand == 3, col = 3, type = "b")
> legend(3, 0.9, paste("Brand", 1:3), col = c(4, 2, 3), lty = 1)
> fm1 <- lm(strength ~ brand + treatment, data = xmp11.01)
> anova(fm1)      # compare to table 11.1, page 439
Analysis of Variance Table

Response: strength
          Df  Sum Sq Mean Sq F value   Pr(>F)   
brand      2 0.12822 0.06411  4.4323 0.065765 . 
treatment  3 0.47969 0.15990 11.0549 0.007399 **
Residuals  6 0.08678 0.01446                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.05"
> 
> ### * xmp11.05
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.05
> ### Title: data from Example 11.5
> ### Aliases: xmp11.05
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.05)
> plot(power ~ humid, data = xmp11.05, col = "lightgray",
+      xlab = "Level of humidity",
+      ylab = "Estimated annual power consumption (kwh)",
+      main = "Data from Example 11.5")
> lines(power ~ as.integer(humid), data = xmp11.05,
+         subset = brand == 1, col = 2, type = "b")
> lines(power ~ as.integer(humid), data = xmp11.05,
+         subset = brand == 2, col = 3, type = "b")
> lines(power ~ as.integer(humid), data = xmp11.05,
+         subset = brand == 3, col = 4, type = "b")
> lines(power ~ as.integer(humid), data = xmp11.05,
+         subset = brand == 4, col = 5, type = "b")
> lines(power ~ as.integer(humid), data = xmp11.05,
+         subset = brand == 5, col = 6, type = "b")
> legend(0.6, 1010, paste("Brand", 1:5), col = 1 + (1:5),
+        lty = 1)
> fm1 <- lm(power ~ humid + brand, data = xmp11.05)
> anova(fm1)    # compare with Table 11.3, page 442
Analysis of Variance Table

Response: power
          Df Sum Sq Mean Sq F value    Pr(>F)    
humid      3 116218   38739 278.199 2.364e-11 ***
brand      4  53231   13308  95.567 5.419e-09 ***
Residuals 12   1671     139                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> summary(fm1)

Call:
lm(formula = power ~ humid + brand, data = xmp11.05)

Residuals:
    Min      1Q  Median      3Q     Max 
-15.650  -4.813  -1.275   7.412  19.150 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 797.5000     5.9002 135.165  < 2e-16 ***
humid.L     151.4712     5.2773  28.702 1.99e-12 ***
humid.Q     -17.3000     5.2773  -3.278 0.006602 ** 
humid.C       0.8497     5.2773   0.161 0.874764    
brand2       46.0000     8.3442   5.513 0.000133 ***
brand3       41.5000     8.3442   4.974 0.000323 ***
brand4      116.5000     8.3442  13.962 8.80e-09 ***
brand5      139.7500     8.3442  16.748 1.09e-09 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 11.8 on 12 degrees of freedom
Multiple R-Squared: 0.9902,	Adjusted R-squared: 0.9845 
F-statistic: 173.8 on 7 and 12 DF,  p-value: 4.233e-11 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.06"
> 
> ### * xmp11.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.06
> ### Title: data from Example 11.6
> ### Aliases: xmp11.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.06)
> plot(Resp ~ Stimulus, data = xmp11.06, col = "lightgray",
+      main = "Data from Example 11.6",
+      ylab = "Mean number of responses")
> for (i in seq(along = levels(xmp11.06$Subject))) {
+     attach(xmp11.06[ xmp11.06$Subject == i, ])
+     lines(Resp ~ as.integer(Stimulus), col = i+1, type = "b")
+ }
> legend(0.8, 95, paste("Subject", levels(xmp11.06$Subject)),
+        col = 1 + seq(along = levels(xmp11.06$Subject)),
+        lty = 1)
> fm1 <- lm(Resp ~ Stimulus + Subject, data = xmp11.06)
> anova(fm1)     # compare to Table 11.5, page 443
Analysis of Variance Table

Response: Resp
          Df  Sum Sq Mean Sq F value    Pr(>F)    
Stimulus   5  1428.3   285.7  5.4888   0.00455 ** 
Subject    3 13444.6  4481.5 86.1112 1.110e-09 ***
Residuals 15   780.7    52.0                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> attach(xmp11.06)
> means <- sort(tapply(Resp, Stimulus, mean))
> means
    L1     L2      T  L1+L2   L1+T   L2+T 
24.825 27.875 29.100 40.350 41.225 45.050 
> diff(means)    # successive differences
    L2      T  L1+L2   L1+T   L2+T 
 3.050  1.225 11.250  0.875  3.825 
> qtukey(0.95, nmeans = 6, df = 15)  #for Tukey comparisons
[1] 4.594735
> detach()
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.07"
> 
> ### * xmp11.07
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.07
> ### Title: data from Example 11.7
> ### Aliases: xmp11.07
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.07)
> plot(Yield ~ Density, data = xmp11.07, col = "lightgray",
+   main = "Data from Example 11.7, page 450",
+   xlab = "Density (plants/hectare)")
> means <- sapply(split(xmp11.07, xmp11.07$Density),
+   function(x) tapply(x$Yield, x$Variety, mean))
> round(means, 2)
      1     2     3     4
1  9.20 12.43 12.90 10.80
2  8.93 12.63 14.50 12.77
3 16.30 18.10 19.93 18.17
> lines(1:4, means[1, ], col = 4, type = "b")
> lines(1:4, means[2, ], col = 2, type = "b")
> lines(1:4, means[3, ], col = 3, type = "b")
> legend(0.4, 21.2, levels(xmp11.07$Variety), lty = 2,
+   col = c(4, 2, 3))
> fm1 <- lm(Yield ~ Variety * Density, data = xmp11.07)
> anova(fm1)    # compare with Table 11.7, page 452
Analysis of Variance Table

Response: Yield
                Df Sum Sq Mean Sq  F value    Pr(>F)    
Variety          2 327.60  163.80 103.3430 1.608e-12 ***
Density          3  86.69   28.90  18.2306 2.212e-06 ***
Variety:Density  6   8.03    1.34   0.8445    0.5484    
Residuals       24  38.04    1.58                       
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> fm2 <- update(fm1, . ~ Variety + Density)  # additive model
> anova(fm2)
Analysis of Variance Table

Response: Yield
          Df Sum Sq Mean Sq F value    Pr(>F)    
Variety    2 327.60  163.80 106.659 2.313e-14 ***
Density    3  86.69   28.90  18.816 4.690e-07 ***
Residuals 30  46.07    1.54                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> sort(tapply(xmp11.07$Yield, xmp11.07$Variety, mean))
       1        2        3 
11.33333 12.20833 18.12500 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.11"
> 
> ### * xmp11.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.11
> ### Title: data from Example 11.11
> ### Aliases: xmp11.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.11)
> coplot(Tempr ~ as.integer(Period) | Strain * Coat,
+        data = xmp11.11, show.given = FALSE)
> coplot(Tempr ~ as.integer(Strain) | Period * Coat,
+        data = xmp11.11, show.given = FALSE)
> coplot(Tempr ~ as.integer(Coat) | Period * Strain,
+        data = xmp11.11, show.given = FALSE)
> fm1 <- lm(Tempr ~ Period * Strain * Coat, xmp11.11)
> anova(fm1)  # compare with Table 11.8, page 461
Analysis of Variance Table

Response: Tempr
                   Df Sum Sq Mean Sq  F value    Pr(>F)    
Period              3  0.490   0.163   4.1398 0.0095807 ** 
Strain              1  6.458   6.458 163.5910 < 2.2e-16 ***
Coat                3 48.934  16.311 413.1671 < 2.2e-16 ***
Period:Strain       3  0.020   0.007   0.1715 0.9152673    
Period:Coat         9  1.608   0.179   4.5268 0.0001254 ***
Strain:Coat         3  0.876   0.292   7.3975 0.0002489 ***
Period:Strain:Coat  9  0.250   0.028   0.7039 0.7029981    
Residuals          64  2.527   0.039                       
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.12"
> 
> ### * xmp11.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.12
> ### Title: data from Example 11.12
> ### Aliases: xmp11.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.12)
> xmp11.12$row <- ordered(xmp11.12$row)
> xmp11.12$column <- ordered(xmp11.12$column)
> xmp11.12$humidity <- ordered(xmp11.12$humidity)
> attach(xmp11.12)  # to check the design
> table(row, column)
   column
row 1 2 3 4 5 6
  1 1 1 1 1 1 1
  2 1 1 1 1 1 1
  3 1 1 1 1 1 1
  4 1 1 1 1 1 1
  5 1 1 1 1 1 1
  6 1 1 1 1 1 1
> table(row, humidity)
   humidity
row 25% 37% 50% 62% 75% 87%
  1   1   1   1   1   1   1
  2   1   1   1   1   1   1
  3   1   1   1   1   1   1
  4   1   1   1   1   1   1
  5   1   1   1   1   1   1
  6   1   1   1   1   1   1
> table(humidity, column)
        column
humidity 1 2 3 4 5 6
     25% 1 1 1 1 1 1
     37% 1 1 1 1 1 1
     50% 1 1 1 1 1 1
     62% 1 1 1 1 1 1
     75% 1 1 1 1 1 1
     87% 1 1 1 1 1 1
> detach()
> fm1 <- lm(abrasion ~ row + column + humidity, xmp11.12)
> anova(fm1)    # compare with Table 11.9, page 464
Analysis of Variance Table

Response: abrasion
          Df  Sum Sq Mean Sq F value   Pr(>F)    
row        5  2.1897  0.4379  2.5106  0.06407 .  
column     5  2.5743  0.5149  2.9516  0.03731 *  
humidity   5 23.5301  4.7060 26.9789 3.03e-08 ***
Residuals 20  3.4887  0.1744                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.13"
> 
> ### * xmp11.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.13
> ### Title: data from Example 11.13
> ### Aliases: xmp11.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.13)
> fm1 <- lm(Strength ~ age * tempture * soil, xmp11.13)
> anova(fm1)   # compare with Table 11.12, page 471
Analysis of Variance Table

Response: Strength
                  Df Sum Sq Mean Sq  F value    Pr(>F)    
age                1 252255  252255 117.9210 4.573e-06 ***
tempture           1  28985   28985  13.5496   0.00621 ** 
soil               1   2328    2328   1.0883   0.32735    
age:tempture       1   3393    3393   1.5861   0.24338    
age:soil           1   1425    1425   0.6662   0.43802    
tempture:soil      1    315     315   0.1473   0.71115    
age:tempture:soil  1   3335    3335   1.5590   0.24711    
Residuals          8  17114    2139                       
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp11.16"
> 
> ### * xmp11.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp11.16
> ### Title: data from Example 11.16
> ### Aliases: xmp11.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp11.16)
> ## leave -1/+1 encoding for experimental factors, convert block
> fm1 <- aov(strength ~ tempture * gas * conc + block,
+            data = xmp11.16)
> summary(fm1)     # anova table
              Df  Sum Sq Mean Sq F value    Pr(>F)    
tempture       1 12996.0 12996.0 39.8243 0.0007389 ***
gas            1   702.2   702.2  2.1519 0.1927565    
conc           1  2756.2  2756.2  8.4461 0.0271136 *  
block          3  5204.0  1734.7  5.3156 0.0398310 *  
tempture:gas   1   210.3   210.3  0.6443 0.4527783    
tempture:conc  1    30.2    30.2  0.0927 0.7710602    
gas:conc       1    25.0    25.0  0.0766 0.7912347    
Residuals      6  1958.0   326.3                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.01"
> 
> ### * xmp12.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.01
> ### Title: data from Example 12.1
> ### Aliases: xmp12.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.01)
> plot(OSA ~ palprebal, data = xmp12.01,
+      xlab = "Palprebal fissure width (cm)",
+      ylab = expression(paste(plain("Occular surface area (cm")^2,
+                plain(")"))),
+      main = "Data from Example 12.1, page 490", las = 1)
> summary(xmp12.01)
   palprebal          OSA       
 Min.   :0.400   Min.   :0.880  
 1st Qu.:0.750   1st Qu.:1.755  
 Median :1.135   Median :3.085  
 Mean   :1.051   Mean   :2.840  
 3rd Qu.:1.363   3rd Qu.:3.808  
 Max.   :1.600   Max.   :4.920  
> fm1 <- lm(OSA ~ palprebal, data = xmp12.01)
> summary(fm1)

Call:
lm(formula = OSA ~ palprebal, data = xmp12.01)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.60942 -0.19875 -0.01902  0.21727  0.66378 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  -0.3977     0.1680  -2.367   0.0251 *  
palprebal     3.0800     0.1506  20.453   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.308 on 28 degrees of freedom
Multiple R-Squared: 0.9373,	Adjusted R-squared: 0.935 
F-statistic: 418.3 on 1 and 28 DF,  p-value: < 2.2e-16 

> abline(fm1)
> opar <- par(mfrow = c(2,2))
> plot(fm1)
> par(opar)
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp12.02"
> 
> ### * xmp12.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.02
> ### Title: data from Example 12.2
> ### Aliases: xmp12.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.02)
> plot(dieback ~ soil.pH, data = xmp12.02,
+      xlab = "soil pH", ylab = "mean crown dieback (%)",
+      main = "Data from Example 12.2, page 491")
> fm1 <- lm(dieback ~ soil.pH, data = xmp12.02)
> abline(fm1)
> summary(fm1)

Call:
lm(formula = dieback ~ soil.pH, data = xmp12.02)

Residuals:
    Min      1Q  Median      3Q     Max 
-4.6253 -1.5520 -0.4994  2.0188  5.6876 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   31.040      5.445   5.701  2.6e-05 ***
soil.pH       -5.792      1.363  -4.250 0.000541 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 2.981 on 17 degrees of freedom
Multiple R-Squared: 0.5151,	Adjusted R-squared: 0.4866 
F-statistic: 18.06 on 1 and 17 DF,  p-value: 0.0005406 

> opar <- par(mfrow = c(2,2))
> plot(fm1)
> par(opar)
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp12.04"
> 
> ### * xmp12.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.04
> ### Title: data from Example 12.4
> ### Aliases: xmp12.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.04)
> plot(porosity ~ weight, data = xmp12.04,
+      xlab = "Unit weight (pcf) of concrete specimen",
+      ylab = "Porosity (%)",
+      main = "Data from Example 12.4, page 500")
> fm1 <- lm(porosity ~ weight, data = xmp12.04)
> abline(fm1)
> summary(fm1)

Call:
lm(formula = porosity ~ weight, data = xmp12.04)

Residuals:
    Min      1Q  Median      3Q     Max 
-1.7754 -0.5727 -0.1325  0.6034  1.6818 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 118.90992    4.49912   26.43 1.10e-12 ***
weight       -0.90473    0.04109  -22.02 1.12e-11 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.938 on 13 degrees of freedom
Multiple R-Squared: 0.9739,	Adjusted R-squared: 0.9719 
F-statistic: 484.8 on 1 and 13 DF,  p-value: 1.125e-11 

> opar <- par(mfrow = c(2,2))
> plot(fm1)
> par(opar)
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp12.06"
> 
> ### * xmp12.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.06
> ### Title: data from Example 12.6
> ### Aliases: xmp12.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.06)
> plot(lead ~ traffic, data = xmp12.06,
+      xlab = "Traffic flow (1000's of cars per 24 hours)",
+      ylab = expression(paste(plain("Lead content of tree bark ("),
+        mu,plain("g/g dry wt)"))),
+      main = "Data from Example 12.6, page 503", las = 1)
> fm1 <- lm(lead ~ traffic, data = xmp12.06)
> abline(fm1)
> summary(fm1)

Call:
lm(formula = lead ~ traffic, data = xmp12.06)

Residuals:
    Min      1Q  Median      3Q     Max 
-128.43  -63.13   24.52   69.32  125.72 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  -12.842     72.143  -0.178    0.863    
traffic       36.184      3.693   9.798 4.24e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 92.19 on 9 degrees of freedom
Multiple R-Squared: 0.9143,	Adjusted R-squared: 0.9048 
F-statistic: 96.01 on 1 and 9 DF,  p-value: 4.239e-06 

> opar <- par(mfrow = c(2, 2))
> plot(fm1)
> par(opar)
> ## compare to table on page 503
> cbind(xmp12.06, yhat = fitted(fm1), resid = resid(fm1))
   traffic lead      yhat      resid
1      8.3  227  287.4844  -60.48439
2      8.3  312  287.4844   24.51561
3     12.1  362  424.9830  -62.98301
4     12.1  521  424.9830   96.01699
5     17.0  640  602.2839   37.71613
6     17.0  539  602.2839  -63.28387
7     17.0  728  602.2839  125.71613
8     24.3  945  866.4260   78.57404
9     24.3  738  866.4260 -128.42596
10    24.3  759  866.4260 -107.42596
11    33.6 1263 1202.9357   60.06426
> anova(fm1)
Analysis of Variance Table

Response: lead
          Df Sum Sq Mean Sq F value    Pr(>F)    
traffic    1 815966  815966  96.005 4.239e-06 ***
Residuals  9  76493    8499                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp12.08"
> 
> ### * xmp12.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.08
> ### Title: data from Example 12.8
> ### Aliases: xmp12.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.08)
> plot(attenuat ~ strength, data = xmp12.08,
+      xlab = "Fracture strength (% of ultimate tensile strength)",
+      ylab = "Attenuation (neper/cm)",
+      main = "Data from Example 12.8, page 504")
> fm1 <- lm(attenuat ~ strength, data = xmp12.08)
> abline(fm1)
> summary(fm1)

Call:
lm(formula = attenuat ~ strength, data = xmp12.08)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.23808 -0.07052 -0.02065  0.09414  0.30869 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  3.620906   0.099490   36.40 1.19e-13 ***
strength    -0.014711   0.001436  -10.24 2.76e-07 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.1479 on 12 degrees of freedom
Multiple R-Squared: 0.8974,	Adjusted R-squared: 0.8888 
F-statistic: 104.9 on 1 and 12 DF,  p-value: 2.762e-07 

> opar <- par(mfrow = c(2, 2))
> plot(fm1)
> par(opar)
> anova(fm1)
Analysis of Variance Table

Response: attenuat
          Df  Sum Sq Mean Sq F value    Pr(>F)    
strength   1 2.29469 2.29469  104.92 2.762e-07 ***
Residuals 12 0.26246 0.02187                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> graphics::par(get("par.postscript", env = .CheckExEnv))
> cleanEx(); ..nameEx <- "xmp12.10"
> 
> ### * xmp12.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.10
> ### Title: data from Example 12.10
> ### Aliases: xmp12.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.11"
> 
> ### * xmp12.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.11
> ### Title: data from Example 12.11
> ### Aliases: xmp12.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.11)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.12"
> 
> ### * xmp12.12
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.12
> ### Title: data from Example 12.12
> ### Aliases: xmp12.12
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.12)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.13"
> 
> ### * xmp12.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.13
> ### Title: data from Example 12.13
> ### Aliases: xmp12.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.13)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.14"
> 
> ### * xmp12.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.14
> ### Title: data from Example 12.14
> ### Aliases: xmp12.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.15"
> 
> ### * xmp12.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.15
> ### Title: data from Example 12.15
> ### Aliases: xmp12.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp12.16"
> 
> ### * xmp12.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp12.16
> ### Title: data from Example 12.16
> ### Aliases: xmp12.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp12.16)
> plot(y ~ x, data = xmp12.16, xlab = "Ozone concentration (ppm)",
+      ylab = "Secondary carbon concentration")
> cor(xmp12.16$x, xmp12.16$y)
[1] 0.7155355
> cor.test(~ x + y, data = xmp12.16)

	Pearson's product-moment correlation

data:  x and y 
t = 3.8325, df = 14, p-value = 0.001830
alternative hypothesis: true correlation is not equal to 0 
95 percent confidence interval:
 0.3406608 0.8941061 
sample estimates:
      cor 
0.7155355 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.01"
> 
> ### * xmp13.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.01
> ### Title: data from Example 13.1
> ### Aliases: xmp13.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.01)
> plot(emission ~ rate, data = xmp13.01,
+      xlab = "Burner area liberation rate",
+      ylab = expression(plain("NO")["x"]*
+                  plain("emissions")), las = 1,
+      main = "Data from Example 13.1, page 545")
> fm1 <- lm(emission ~ rate, data = xmp13.01)
> abline(fm1)                   # plot 1, Figure 13.1
> if (require(MASS)) {
+    sres <- stdres(fm1)
+    plot(sres ~ fitted(fm1), ylab = "Standardized residuals",
+        xlab = "Fitted values")  # plot 2, Figure 13.1
+    abline(h = 0, lty = 2, lwd = 0)   # horizontal reference
+ }
Loading required package: MASS
> plot(resid(fm1) ~ fitted(fm1), ylab = "Residuals",
+      xlab = "Fitted values")  # alternative plot 2
> abline(h = 0, lty = 2, lwd = 0)   # horizontal reference
> plot(fitted(fm1) ~ emission, data = xmp13.01)
> abline(0, 1)                  # plot 3, Figure 13.1
> if (require(MASS)) {
+    plot(sres ~ rate, data = xmp13.01) # plot 4
+    abline(h = 0, lty = 2, lwd = 0)
+    qqnorm(sres)               # plot 5
+ } else {
+    plot(resid(fm1) ~ rate, data = xmp13.01) # plot 4
+    qqnorm(resid(fm1))         # plot 5
+ }
> ## The residuals versus fitted plot and the normal
> ## probability plot of the standardized residuals 
> plot(fm1, which = 1:2)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.03"
> 
> ### * xmp13.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.03
> ### Title: data from Example 13.3
> ### Aliases: xmp13.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.03)
> plot(ToolLife ~ time, data = xmp13.03)
> plot(ToolLife ~ time, data = xmp13.03, log = "xy")
> fm1 <- lm(log(ToolLife) ~ I(log(time)), data = xmp13.03)
> summary(fm1)

Call:
lm(formula = log(ToolLife) ~ I(log(time)), data = xmp13.03)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.27312 -0.17695 -0.05837  0.15764  0.69134 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)   35.6693     3.0731   11.61 3.99e-07 ***
I(log(time))  -5.3997     0.4954  -10.90 7.18e-07 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.2846 on 10 degrees of freedom
Multiple R-Squared: 0.9224,	Adjusted R-squared: 0.9146 
F-statistic: 118.8 on 1 and 10 DF,  p-value: 7.179e-07 

> plot(fm1, which = 1)     # plot of residuals versus fitted values
> plot(exp(fitted(fm1)) ~ xmp13.03$ToolLife,
+    xlab = "y", ylab = expression(hat("y")),
+    main = "Compare to Figure 13.4, page 555")
> abline(0, 1)       # reference line
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.04"
> 
> ### * xmp13.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.04
> ### Title: data from Example 13.4
> ### Aliases: xmp13.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.04)
> plot(Ethylene ~ time, data = xmp13.04,
+      xlab = "Exposure time (min)",
+      ylab = "Ethylene content (nL/g dry wt)",
+      main = "Compare to Figure 13.5 (a), page 556")
> fm1 <- lm(Ethylene ~ time, data = xmp13.04)
> abline(fm1)
> plot(resid(fm1) ~ xmp13.04$time)
> abline(h = 0, lty = 2)
> title(main = "Compare to Figure 13.5 (b), page 556")
> title(sub = "Using raw residuals instead of standardized")
> fm2 <- lm(log(Ethylene) ~ time, data = xmp13.04)
> plot(resid(fm2) ~ xmp13.04$time)
> abline(h = 0, lty = 2)
> title(main = "Compare to Figure 13.6 (a), page 557")
> title(sub = "Using raw residuals instead of standardized")
> summary(fm2)

Call:
lm(formula = log(Ethylene) ~ time, data = xmp13.04)

Residuals:
       Min         1Q     Median         3Q        Max 
-0.1475375 -0.0332298  0.0004249  0.0398234  0.1354295 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)  5.9404951  0.0443816   133.8 3.68e-16 ***
time        -0.0323287  0.0007501   -43.1 9.73e-12 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.07797 on 9 degrees of freedom
Multiple R-Squared: 0.9952,	Adjusted R-squared: 0.9946 
F-statistic:  1857 on 1 and 9 DF,  p-value: 9.734e-12 

> plot(exp(fitted(fm2)) ~ xmp13.04$Ethylene)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.06"
> 
> ### * xmp13.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.06
> ### Title: data from Example 13.6
> ### Aliases: xmp13.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.06)
> fm1 <- glm(Failure ~ Temperature, 
+     data = xmp13.06, family = "binomial")
> ## results are different from JMP results in Figure 13.9
> summary(fm1)

Call:
glm(formula = Failure ~ Temperature, family = "binomial", data = xmp13.06)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-1.12489  -0.72335  -0.40139  -0.04056   2.22441  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept) 11.74641    6.02142   1.951   0.0511 .
Temperature -0.18843    0.08909  -2.115   0.0344 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 26.992  on 23  degrees of freedom
Residual deviance: 20.482  on 22  degrees of freedom
AIC: 24.482

Number of Fisher Scoring iterations: 5

> temp <- seq(55, 85, len = 101)   # for doing the prediction
> plot(
+   predict(fm1, new = list(Temperature = temp), type = "resp") ~ temp,
+   main = "Compare with Figure 13.8, page 560", type = "l",
+   ylab = "P(F)")
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.07"
> 
> ### * xmp13.07
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.07
> ### Title: data from Example 13.6
> ### Aliases: xmp13.07
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.07)
> plot(yield ~ days, data = xmp13.07,
+    main = "Compare to Figure 13.11(a), page 579")
> fm1 <- lm(yield ~ days + I(days^2), data = xmp13.07)
> summary(fm1)

Call:
lm(formula = yield ~ days + I(days^2), data = xmp13.07)

Residuals:
    Min      1Q  Median      3Q     Max 
-303.96 -118.11   13.86  115.67  319.06 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept) -1070.3977   617.2527  -1.734    0.107    
days          293.4829    42.1776   6.958 9.94e-06 ***
I(days^2)      -4.5358     0.6744  -6.726 1.41e-05 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 203.9 on 13 degrees of freedom
Multiple R-Squared: 0.7942,	Adjusted R-squared: 0.7625 
F-statistic: 25.08 on 2 and 13 DF,  p-value: 3.452e-05 

> anova(fm1)
Analysis of Variance Table

Response: yield
          Df  Sum Sq Mean Sq F value    Pr(>F)    
days       1  204526  204526  4.9202   0.04497 *  
I(days^2)  1 1880253 1880253 45.2328 1.414e-05 ***
Residuals 13  540388   41568                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> predict(fm1, list(days = 25), interval = "conf")
        fit      lwr      upr
[1,] 3431.8 3282.275 3581.325
> predict(fm1, list(days = 25), interval = "pred")
        fit      lwr      upr
[1,] 3431.8 2966.649 3896.951
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.09"
> 
> ### * xmp13.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.09
> ### Title: data from Example 13.9
> ### Aliases: xmp13.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.09)
> plot(strength ~ tempture, data = xmp13.09)
> fm1 <- lm(strength ~ tempture + I(tempture^2), data = xmp13.09)
> summary(fm1)

Call:
lm(formula = strength ~ tempture + I(tempture^2), data = xmp13.09)

Residuals:
      1       2       3       4       5       6       7       8 
-22.542  -2.187  50.317  35.935 -17.485 -38.948 -47.497  42.407 

Coefficients:
                Estimate Std. Error t value Pr(>|t|)  
(Intercept)   -2.622e+04  1.191e+04  -2.201   0.0790 .
tempture       1.892e+02  8.024e+01   2.358   0.0649 .
I(tempture^2) -3.312e-01  1.350e-01  -2.453   0.0577 .
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 45.2 on 5 degrees of freedom
Multiple R-Squared: 0.8731,	Adjusted R-squared: 0.8224 
F-statistic:  17.2 on 2 and 5 DF,  p-value: 0.005734 

> xmp13.09$Tcentered <- scale(xmp13.09$tempture, scale = FALSE)
> fm2 <- lm(strength ~ Tcentered + I(Tcentered^2), data = xmp13.09)
> summary(fm2)

Call:
lm(formula = strength ~ Tcentered + I(Tcentered^2), data = xmp13.09)

Residuals:
      1       2       3       4       5       6       7       8 
-22.542  -2.187  50.317  35.935 -17.485 -38.948 -47.497  42.407 

Coefficients:
               Estimate Std. Error t value Pr(>|t|)    
(Intercept)    759.3948    23.2009  32.731    5e-07 ***
Tcentered       -7.6073     1.4315  -5.314  0.00315 ** 
I(Tcentered^2)  -0.3312     0.1350  -2.453  0.05769 .  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 45.2 on 5 degrees of freedom
Multiple R-Squared: 0.8731,	Adjusted R-squared: 0.8224 
F-statistic:  17.2 on 2 and 5 DF,  p-value: 0.005734 

> ## another approach using orthogonal polynomials
> fm3 <- lm(strength ~ poly(tempture, 2), data = xmp13.09)
> summary(fm3)

Call:
lm(formula = strength ~ poly(tempture, 2), data = xmp13.09)

Residuals:
      1       2       3       4       5       6       7       8 
-22.542  -2.187  50.317  35.935 -17.485 -38.948 -47.497  42.407 

Coefficients:
                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)          718.12      15.98  44.942 1.03e-07 ***
poly(tempture, 2)1  -240.81      45.20  -5.328  0.00312 ** 
poly(tempture, 2)2  -110.88      45.20  -2.453  0.05769 .  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 45.2 on 5 degrees of freedom
Multiple R-Squared: 0.8731,	Adjusted R-squared: 0.8224 
F-statistic:  17.2 on 2 and 5 DF,  p-value: 0.005734 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.10"
> 
> ### * xmp13.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.10
> ### Title: data from Example 13.10
> ### Aliases: xmp13.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.11"
> 
> ### * xmp13.11
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.11
> ### Title: data from Example 13.11
> ### Aliases: xmp13.11
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.11)
> fm1 <- lm(Strength ~ Force + Power + Temperature + Time,
+      data = xmp13.11)
> summary(fm1)

Call:
lm(formula = Strength ~ Force + Power + Temperature + Time, data = xmp13.11)

Residuals:
     Min       1Q   Median       3Q      Max 
-11.0900  -1.7608  -0.3067   2.4392   7.5933 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) -37.47667   13.09964  -2.861  0.00841 ** 
Force         0.21167    0.21057   1.005  0.32444    
Power         0.49833    0.07019   7.100 1.93e-07 ***
Temperature   0.12967    0.04211   3.079  0.00499 ** 
Time          0.25833    0.21057   1.227  0.23132    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 5.158 on 25 degrees of freedom
Multiple R-Squared: 0.714,	Adjusted R-squared: 0.6682 
F-statistic:  15.6 on 4 and 25 DF,  p-value: 1.592e-06 

> anova(fm1)
Analysis of Variance Table

Response: Strength
            Df  Sum Sq Mean Sq F value    Pr(>F)    
Force        1   26.88   26.88  1.0104  0.324436    
Power        1 1341.02 1341.02 50.4051 1.928e-07 ***
Temperature  1  252.20  252.20  9.4796  0.004992 ** 
Time         1   40.04   40.04  1.5051  0.231324    
Residuals   25  665.12   26.60                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.13"
> 
> ### * xmp13.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.13
> ### Title: data from Example 13.13
> ### Aliases: xmp13.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.13)
> fm1 <- lm(strength ~ x1 * x2, data = xmp13.13)
> summary(fm1)

Call:
lm(formula = strength ~ x1 * x2, data = xmp13.13)

Residuals:
      1       2       3       4       5       6       7       8       9 
 0.3417 -0.1750 -2.4583 -2.9750 -0.7167  2.0833  0.2000  0.7167  2.9833 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)  
(Intercept)    6.217     30.304   0.205   0.8455  
x1             5.779      2.079   2.779   0.0389 *
x2            51.333     50.434   1.018   0.3555  
x1:x2         -9.357      3.461  -2.704   0.0426 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 2.423 on 5 degrees of freedom
Multiple R-Squared: 0.8946,	Adjusted R-squared: 0.8314 
F-statistic: 14.15 on 3 and 5 DF,  p-value: 0.00706 

> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.15"
> 
> ### * xmp13.15
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.15
> ### Title: data from Example 13.15
> ### Aliases: xmp13.15
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.15)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.16"
> 
> ### * xmp13.16
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.16
> ### Title: data from Example 13.16
> ### Aliases: xmp13.16
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.16)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.18"
> 
> ### * xmp13.18
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.18
> ### Title: data from Example 13.18
> ### Aliases: xmp13.18
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.18)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.19"
> 
> ### * xmp13.19
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.19
> ### Title: data from Example 13.19
> ### Aliases: xmp13.19
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.19)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp13.22"
> 
> ### * xmp13.22
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp13.22
> ### Title: data from Example 13.22
> ### Aliases: xmp13.22
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp13.22)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp14.03"
> 
> ### * xmp14.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp14.03
> ### Title: data from Example 14.3
> ### Aliases: xmp14.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp14.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp14.10"
> 
> ### * xmp14.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp14.10
> ### Title: data from Example 14.10
> ### Aliases: xmp14.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp14.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp14.13"
> 
> ### * xmp14.13
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp14.13
> ### Title: data from Example 14.13
> ### Aliases: xmp14.13
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp14.13)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp14.14"
> 
> ### * xmp14.14
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp14.14
> ### Title: data from Example 14.14
> ### Aliases: xmp14.14
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp14.14)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.01"
> 
> ### * xmp15.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.01
> ### Title: data from Example 15.1
> ### Aliases: xmp15.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.01)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.02"
> 
> ### * xmp15.02
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.02
> ### Title: data from Example 15.2
> ### Aliases: xmp15.02
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.02)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.03"
> 
> ### * xmp15.03
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.03
> ### Title: data from Example 15.3
> ### Aliases: xmp15.03
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.03)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.04"
> 
> ### * xmp15.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.04
> ### Title: data from Example 15.4
> ### Aliases: xmp15.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.06"
> 
> ### * xmp15.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.06
> ### Title: data from Example 15.6
> ### Aliases: xmp15.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.08"
> 
> ### * xmp15.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.08
> ### Title: data from Example 15.8
> ### Aliases: xmp15.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.09"
> 
> ### * xmp15.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.09
> ### Title: data from Example 15.9
> ### Aliases: xmp15.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.09)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp15.10"
> 
> ### * xmp15.10
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp15.10
> ### Title: data from Example 15.10
> ### Aliases: xmp15.10
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp15.10)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.01"
> 
> ### * xmp16.01
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.01
> ### Title: data from Example 16.1
> ### Aliases: xmp16.01
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.01)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.04"
> 
> ### * xmp16.04
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.04
> ### Title: data from Example 16.4
> ### Aliases: xmp16.04
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.04)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.06"
> 
> ### * xmp16.06
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.06
> ### Title: data from Example 16.6
> ### Aliases: xmp16.06
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.06)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.07"
> 
> ### * xmp16.07
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.07
> ### Title: data from Example 16.7
> ### Aliases: xmp16.07
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.07)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.08"
> 
> ### * xmp16.08
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.08
> ### Title: data from Example 16.8
> ### Aliases: xmp16.08
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.08)
> 
> 
> 
> cleanEx(); ..nameEx <- "xmp16.09"
> 
> ### * xmp16.09
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: xmp16.09
> ### Title: data from Example 16.9
> ### Aliases: xmp16.09
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(xmp16.09)
> 
> 
> 
> ### * <FOOTER>
> ###
> cat("Time elapsed: ", proc.time() - get("ptime", env = .CheckExEnv),"\n")
Time elapsed:  10.45 0.65 11.5 0 0 
> grDevices::dev.off()
null device 
          1 
> ###
> ### Local variables: ***
> ### mode: outline-minor ***
> ### outline-regexp: "\\(> \\)?### [*]+" ***
> ### End: ***
> quit('no')