|
|
|
||||
|
|
|||||
Figure 2.1
qqmath(~ height | voice.part,
distribution=qunif,
data=singer,
panel = function(x, y) {
panel.grid()
panel.xyplot(x, y)
},
layout=c(2,4),
aspect=1,
sub = list("Figure 2.1",cex=.8),
xlab = "f-value",
ylab="Height (inches)")
Figure 2.2
qqmath(~ sort(singer$height[singer$voice.part=="Tenor 1"]),
distribution = qunif,
panel = function(x, y) {
panel.qqmath(c(0, x, 1), c(min(y), y, max(y)), type = "l")
panel.qqmath(x, y, col = 0, pch = 16)
panel.qqmath(x, y)
},
aspect = 1,
sub = list("Figure 2.2",cex=.8),
xlab = "f-value",
ylab = "Tenor 1 Height (inches)")
Figure 2.3
voice.part <- ordered(singer$voice.part,
c("Soprano 1", "Soprano 2", "Alto 1", "Alto 2",
"Tenor 1", "Tenor 2", "Bass 1", "Bass 2"))
qq(voice.part ~ singer$height,
subset=voice.part=="Bass 2" | voice.part=="Tenor 1",
aspect=1,
sub = list("Figure 2.3",cex=.8),
xlab = "Tenor 1 Height (inches)",
ylab = "Base 2 Height (inches)")
Figure 2.4
voice.part <- ordered(singer$voice.part,
c("Soprano 1", "Soprano 2", "Alto 1", "Alto 2",
"Tenor 1", "Tenor 2", "Bass 1", "Bass 2"))
bass.tenor.qq <- qq(voice.part ~ singer$height,
subset=voice.part=="Bass 2" | voice.part=="Tenor 1")
tmd(bass.tenor.qq,
aspect=1,
ylab = "Difference (inches)",
sub = list("Figure 2.4",cex=.8),
xlab = "Mean (inches)")
Figure 2.5 missing
Figure 2.6
oldpty <- par("pty")
par(pty = "s")
data <-
c(0.9, 1.6, 2.26305, 2.55052, 2.61059, 2.69284, 2.78511, 2.80955,
2.94647, 2.96043, 3.05728, 3.15748, 3.18033, 3.20021,
3.20156, 3.24435, 3.33231, 3.34176, 3.3762, 3.39578, 3.4925,
3.55195, 3.56207, 3.65149, 3.72746, 3.73338, 3.73869,
3.80469, 3.85224, 3.91386, 3.93034, 4.02351, 4.03947,
4.05481, 4.10111, 4.26249, 4.28782, 4.37586, 4.48811,
4.6001, 4.65677, 4.66167, 4.73211, 4.80803, 4.9812, 5.17246,
5.3156, 5.35086, 5.36848, 5.48167, 5.68, 5.98848, 6.2, 7.1,
7.4)
boxplot(data, rep(NA, length(data)), ylab = "Data")
usr <- par("usr")
x <- usr[1] + (usr[2] - usr[1]) * 0.5
at <- c(0.9, 1.6, 3.2, 3.8, 4.65, 6.2, 7.2)
arrows(rep(x * 1.15, 7), at, rep(x, 7), at)
mtext("Figure 2.6",1,1,cex=.8)
text(rep(x * 1.2, 7), at, adj = 0,
labels = c("outside value", "lower adjacent value",
"lower quartile", "median", "upper quartile",
"upper adjacent value", "outside values"))
par(pty = oldpty)
invisible()
Figure 2.7
data <- round(c(0.9, 1.6, 2.263047,
2.550518, 2.610592, 2.69284, 2.785113,
2.809547, 2.946467, 2.96044, 3.057283,
3.15748, 3.180327, 3.200206,
3.20156, 3.244347, 3.332312,
3.341763, 3.3762, 3.395778, 3.492497,
3.551945, 3.562066, 3.65149,
3.7274632, 3.73338, 3.738686, 3.80469,
3.85224, 3.91386, 3.93034,
4.02351, 4.039466, 4.05481, 4.101108, 4.262486,
4.28782, 4.375864, 4.48811, 4.6001,
4.656775, 4.661673, 4.73211,
4.80803, 4.9812, 5.172464,
5.3156, 5.35086, 5.36848,
5.48167, 5.68, 5.98848, 6.2,
7.1, 7.4),5)
uq <- quantile(data,.75)
lq <- quantile(data,.25)
r <- 1.5*(uq-lq)
h <- c(lq-r,1.6,lq,uq,6.2,uq+r)
writing <- c("lower quartile - 1.5 r",
"lower adjacent value",
"lower quartile",
"upper quartile",
"upper adjacent value",
"upper quartile + 1.5 r")
qqmath(~ data,
distribution = qunif,
panel = substitute(function(x, y) {
reference.line <- trellis.par.get("reference.line")
panel.abline(h = h, lwd = reference.line$lwd, lty = reference.line$lty, col = reference.line$col)
panel.qqmath(x, y, col = 0, pch = 16)
panel.qqmath(x, y)
text(rep(0,3), h[4:6], writing[4:6], adj=0)
text(rep(1,3), h[1:3], writing[1:3], adj=1)
}),
aspect = 1,
sub = list("Figure 2.7",cex=.8),
xlab = "f-value",
ylab = "Data")
Figure 2.8
bwplot(voice.part ~ height,
data=singer,
aspect=1,
sub = list("Figure 2.8",cex=.8),
xlab="Height (inches)")
Figure 2.9
data <- sort(singer$height[singer$voice.part=="Alto 1"])
qqmath(~ data,
distribution = qunif,
panel = function(x, y) {
panel.grid()
panel.qqmath(c(0, x, 1), c(min(y), y, max(y)), type = "l")
panel.qqmath(x, y, col = 0, pch = 16)
panel.qqmath(x, y)
},
aspect = 1,
ylim = range(data, qnorm(ppoints(data), mean(data), sqrt(var(data)))),
sub = list("Figure 2.9",cex=.8),
xlab = "f-value",
ylab = "Alto 1 Height (inches)")
Figure 2.10
data <- sort(singer$height[singer$voice.part=="Alto 1"])
x <- ppoints(data)
y <- qnorm(x, mean(data), sqrt(var(data)))
xyplot(y ~ x,
panel = function(x, y){
panel.grid()
panel.xyplot(x, y, type = "l")
},
ylim = range(data, y),
aspect = 1,
sub = list("Figure 2.10",cex=.8),
xlab = "f-value",
ylab = "Normal Quantile Function")
Figure 2.11
qqmath(~ height | voice.part,
data=singer,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
layout=c(2,4),
aspect=1,
sub = list("Figure 2.11",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Height (inches)")
Figure 2.12
dotplot(tapply(singer$height,singer$voice.part,mean),
aspect=1,
sub = list("Figure 2.12",cex=.8),
xlab="Mean Height (inches)")
Figure 2.13
bwplot(voice.part ~ oneway(height~voice.part, spread = 1)$residuals,
data = singer,
aspect=0.75,
panel = function(x,y){
panel.bwplot(x,y)
panel.abline(v=0)
},
sub = list("Figure 2.13",cex=.8),
xlab = "Residual Height (inches)")
Figure 2.14
res.height <- oneway(height ~ voice.part, data = singer, spread = 1)$residuals
qqmath(~ res.height | singer$voice.part,
distribution = substitute(function(p) quantile(res.height, p)),
panel=function(x,y){
panel.grid()
panel.abline(0, 1)
panel.qqmath(x, y)
},
aspect=1,
layout=c(2,4),
sub = list("Figure 2.14",cex=.8),
xlab = "Pooled Residual Height (inches)",
ylab = "Residual Height (inches)")
Figure 2.15
qqmath(~ oneway(height ~ voice.part, spread = 1)$residuals,
data = singer,
distribution = qunif,
aspect = 1, sub = list("Figure 2.15",cex=.8),
xlab = "f-value",
ylab = "Residual Height (inches)")
Figure 2.16
qqmath(~ oneway(height~voice.part, spread = 1)$residuals,
data = singer,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect=1,
sub = list("Figure 2.16",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Residual Height (inches)")
Figure 2.17
rfs(oneway(height~voice.part, data = singer, spread = 1),
aspect=1,
sub = list("Figure 2.17",cex=.8),
ylab = "Height (inches)")
Figure 2.18 missing
Figure 2.19
qqmath(~ time | nv.vv,
data=fusion.time,
distribution = qunif,
panel = function(x, y) {
panel.grid()
panel.qqmath(x, y)
},
aspect=1,
layout=c(2,1),
sub = list("Figure 2.19",cex=.8),
xlab = "f-value",
ylab="Time (seconds)")
Figure 2.20
data <- 5+qnorm(ppoints(25))
qqmath(~data,
distribution = qunif,
panel = function(x, y) {
reference.line <- trellis.par.get("reference.line")
m <- median(y)
segments(c(.1, .9), m, c(.1, .9), quantile(y, c(.1, .9)),
lwd = reference.line$lwd, lty = reference.line$lty, col = reference.line$col)
panel.qqmath(c(0, x, 1), c(min(y), y, max(y)), type = "l")
panel.qqmath(x, y, col = 0, pch = 16)
panel.qqmath(x, y)
panel.abline(h = m)
text(.05, 4.25, "d(0.1)", srt = 90, adj = 0)
text(.85, 5.25, "d(0.9)", srt = 90, adj = 0)
},
aspect = 1,
sub = list("Figure 2.20",cex=.8),
xlab = "f-value",
ylab = "Data")
Figure 2.21
data <- 2 ^ (5 + qnorm(ppoints(25)))
qqmath(~ data,
distribution = qunif,
panel = function(x, y) {
reference.line <- trellis.par.get("reference.line")
m <- median(y)
segments(c(.1, .9), m, c(.1, .9), quantile(y, c(.1, .9)),
lwd = reference.line$lwd, lty = reference.line$lty, col = reference.line$col)
panel.qqmath(c(0, x, 1), c(min(y), y, max(y)), type = "l")
panel.qqmath(x, y, col = 0, pch = 16)
panel.qqmath(x, y)
panel.abline(h = m)
text(.05, 15, "d(0.1)", srt = 90, adj = 0)
text(.85, 40, "d(0.9)", srt = 90, adj = 0)
},
aspect = 1,
sub = list("Figure 2.21",cex=.8),
xlab = "f-value",
ylab = "Data")
Figure 2.22
qqmath(~ time | nv.vv,
data=fusion.time,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect=1,
layout=c(2,1),
sub = list("Figure 2.22",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Time (seconds)")
Figure 2.23
function()
bwplot(nv.vv ~ time,
data=fusion.time,
aspect = .5,
sub = list("Figure 2.23",cex=.8),
xlab="Time (seconds)")
Figure 2.24
qqmath(~ logb(time,2) | nv.vv,
data=fusion.time,
prepanel=prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect=1,
layout=c(2,1),
sub = list("Figure 2.24",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Log Time (log 2 seconds)")
Figure 2.25
fusion.time.m <- oneway(time ~ nv.vv, data=fusion.time, location=median, spread=1)
xyplot(sqrt(abs(residuals(fusion.time.m)))~jitter(fitted.values(fusion.time.m),factor=3),
aspect=1,
panel=substitute(function(x,y){
panel.xyplot(x,y)
srmads <- sqrt(tapply(abs(residuals(fusion.time.m)),
fusion.time$nv.vv, median))
lines(fusion.time.m$location,srmads)
}),
sub = list("Figure 2.25",cex=.8),
xlab="Jittered Median Time (sec)",
ylab="Square Root Absolute Residual Time (square root sec)")
Figure 2.26
fusion.time.m <- oneway(logb(time,2) ~ nv.vv,data=fusion.time, location = median, spread=1)
xyplot(sqrt(abs(residuals(fusion.time.m))) ~ jitter(fitted.values(fusion.time.m),factor=3),
aspect=1,
panel=substitute(function(x,y){
panel.xyplot(x,y)
srmads <- tapply(abs(residuals(fusion.time.m)),
fusion.time$nv.vv,median)
lines(fusion.time.m$location,srmads)
}),
sub = list("Figure 2.26",cex=.8),
xlab="Jittered Median Log Time (log 2 sec)",
ylab="Square Root Absolute Residual Log Time (square root absolute log 2 sec)")
Figure 2.27
qq(nv.vv ~ time,
data = fusion.time,
aspect = 1,
sub = list("Figure 2.27",cex=.8),
xlab="NV Time (seconds)",
ylab="VV Time (seconds)")
Figure 2.28
qq(nv.vv ~ logb(time, 2),
data = fusion.time,
aspect = 1,
sub = list("Figure 2.28",cex=.8),
xlab = "Log NV Time (log 2 seconds)",
ylab = "Log VV Time (log 2 seconds)")
Figure 2.29
tmd(qq(nv.vv ~ logb(time, 2), data = fusion.time),
aspect = 1,
sub = list("Figure 2.29",cex=.8),
xlab = "Mean (log 2 seconds)",
ylab = "Difference (log 2 seconds)")
Figure 2.30
res <- oneway(logb(time,2)~nv.vv, data = fusion.time, spread = 1)$residuals
qqmath(~ res | fusion.time$nv.vv,
distribution = substitute(function(p) quantile(res, p)),
panel=function(x,y){
panel.grid()
panel.abline(0, 1)
panel.qqmath(x, y)
},
aspect=1,
layout=c(2,1),
sub = list("Figure 2.30",cex=.8),
xlab = "Pooled Residual Log Time (log 2 seconds)",
ylab = "Residual Log Time (log 2 seconds)")
Figure 2.31
qqmath(~ oneway(logb(time,2)~nv.vv, data = fusion.time, spread = 1)$residuals,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect = 1,
sub = list("Figure 2.31",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Residual Log Time (log 2 seconds)")
Figure 2.32
rfs(oneway(logb(time, 2)~nv.vv, data = fusion.time, spread = 1),
aspect=1,
sub = list("Figure 2.32",cex=.8),
ylab = "Log Time (log 2 seconds)")
Figure 2.33
attach(fusion.time)
vvtime <- time[nv.vv=="VV"]
transformed <- cbind(outer(vvtime,c(-1,-1/2,-1/4),"^"),log(vvtime),
(outer(vvtime,c(1/4,1/2,1),"^")))
fusion.time.power <- data.frame(transformed=c(transformed),
lambda = factor(rep(c(-1,-1/2,-1/4,0,1/4,1/2,1),rep(length(vvtime),7))))
ans <- qqmath(~ transformed | lambda,
data=fusion.time.power,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid(h = 0)
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect=1,
scale = list(y = "free"),
layout=c(2,4),
sub = list("Figure 2.33",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "VV Time")
detach()
ans
Figure 2.34
qqmath(~ mean.length | dimension,
distribution = qunif,
data=food.web,
panel = function(x, y) {
panel.grid()
panel.qqmath(x, y)
},
layout=c(1,3),
aspect=1,
sub = list("Figure 2.34",cex=.8),
xlab = "f-value",
ylab="Chain Length")
Figure 2.35
foo.m <- oneway(mean.length~dimension, data = food.web, location = median, spread=1)
set.seed(19)
xyplot(sqrt(abs(residuals(foo.m))) ~ jitter(fitted.values(foo.m),factor=3),
aspect=1,
panel = substitute(function(x,y){
panel.xyplot(x,y)
srmads <- tapply(abs(residuals(foo.m)),
food.web$dimension, median)
lines(foo.m$location,srmads)
}),
sub = list("Figure 2.35",cex=.8),
xlab="Jittered Median Chain Length",
ylab="Square Root Absolute Residual Chain Length")
Figure 2.36
qqmath(~ mean.length | dimension,
data=food.web,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
layout=c(1,3),
aspect=1,
sub = list("Figure 2.36",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Chain Length")
Figure 2.37
foo.m <- oneway(logb(mean.length, 2) ~ dimension, data = food.web, location = median, spread = 1)
set.seed(19)
xyplot(sqrt(abs(residuals(foo.m))) ~ jitter(fitted.values(foo.m), factor = 3),
panel = substitute(function(x, y) {
panel.xyplot(x, y)
add.line <- trellis.par.get("add.line")
lines(foo.m$location, tapply(y, food.web$dimension, median),
lty = add.line$lty, lwd = add.line$lwd, col = add.line$col)
}),
aspect = 1,
sub = list("Figure 2.37",cex=.8),
xlab = "Jittered Median Log 2 Chain Length",
ylab = "Square Root Absolute Residual Log 2 Chain Length")
Figure 2.38
qqmath(~ logb(mean.length,2) | dimension,
data=food.web,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
layout=c(1,3),
aspect=1,
sub = list("Figure 2.38",cex=.8),
xlab = "Unit Normal Quantile",
ylab="Log 2 Chain Length")
Figure 2.39
foo.m <- oneway(1/mean.length ~ dimension, data = food.web, location = median, spread = 1)
set.seed(19)
xyplot(sqrt(abs(residuals(foo.m))) ~ jitter(fitted.values(foo.m), factor = 3),
panel = substitute(function(x,y) {
panel.xyplot(x,y)
add.line <- trellis.par.get("add.line")
lines(foo.m$location, tapply(y, food.web$dimension, median),
lty = add.line$lty, lwd = add.line$lwd, col = add.line$col)
}),
aspect = 1,
sub = list("Figure 2.39",cex=.8),
xlab = "Jittered Median Link Fraction",
ylab = "Square Root Absolute Residual Link Fraction")
Figure 2.40
qqmath(~ (1/mean.length) | dimension,
data = food.web,
panel = function(x, y){
panel.grid()
panel.xyplot(x, y)
panel.qqmathline(y, distribution = qnorm)
},
layout = c(1, 3),
aspect = 1,
sub = list("Figure 2.40",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Link Fraction")
Figure 2.41
res <- oneway((1/mean.length)~dimension, data = food.web, spread = 1)$residuals
qqmath(~ res | food.web$dimension,
distribution = substitute(function(p) quantile(res, p)),
panel=function(x,y){
panel.grid()
panel.abline(0, 1)
panel.qqmath(x, y)
},
layout=c(1,3),
aspect=1,
sub = list("Figure 2.41",cex=.8),
xlab = "Pooled Residual Link Fraction",
ylab = "Residual Link Fraction")
Figure 2.42
rfs(oneway((1/mean.length)~dimension, data = food.web, spread = 1),
sub = list("Figure 2.42",cex=.8),
aspect=1,
ylab = "Link Fraction")
Figure 2.43
bwplot(factor(number.runs) ~ logb(empty.space,2),
data=bin.packing,
aspect=1,
sub = list("Figure 2.43",cex=.8),
xlab="Log 2 Empty Space")
Figure 2.44
qqmath(~ logb(empty.space,2) | factor(number.runs),
data = bin.packing,
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
layout = c(3, 4),
sub = list("Figure 2.44",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Log 2 Empty Space")
Figure 2.45
res <- oneway(logb(empty.space,2) ~ number.runs, data = bin.packing,
location = median,
# This next line is the way it should be.
# spread = function(x) median(abs(x-median(x)))
# This next line is the way it is now.
spread = function(x) (quantile(x,.75)-quantile(x,.25))/1.33)$scaled.residuals
qqmath(~ res | factor(bin.packing$number.runs),
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
layout = c(3,4),
sub = list("Figure 2.45",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Spread-Standardized Residual Log 2 Empty Space")
Figure 2.46
attach(bin.packing)
data <- logb(empty.space,2)
bin.packing.m <- oneway(data~number.runs,location = median,
spread=function(x) diff(quantile(x,c(.25,.75)))/1.35)
res <- bin.packing.m$scaled.res[number.runs>1000]
gr <- factor(number.runs[number.runs>1000])
ans <- qqmath(~ res | gr,
distribution = substitute(function(p) quantile(res, p)),
panel = function(x, y){
panel.grid()
panel.abline(0, 1)
panel.qqmath(x, y)
},
aspect = 1,
layout = c(2, 4),
sub = list("Figure 2.46",cex=.8),
xlab = "Pooled Spread-Standardized Residual Log 2 Empty Space",
ylab = "Spread-Standardized Residual Log 2 Empty Space")
detach()
ans
Figure 2.47
bin.packing.m <- oneway(logb(empty.space,2) ~ number.runs, data = bin.packing,
location = median,
spread = function(x) diff(quantile(x,c(.25,.75)))/1.35)
qqmath(~ bin.packing.m$scaled.res[bin.packing$number.runs > 1000],
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y)
},
aspect = 1,
sub = list("Figure 2.47",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Spread-Standardized Residual Log 2 Empty Space")
Figure 2.48
attach(bin.packing)
data <- logb(empty.space,2)
mq <- tapply(data, number.runs, median)
nw <- logb(sort(unique(number.runs)),2)
ans <- xyplot(mq ~ nw,
panel = function(x, y){
panel.xyplot(x, y)
panel.abline(y[11] - x[11]/3, 1/3)
},
aspect = 1,
sub = list("Figure 2.48",cex=.8),
xlab = "Log 2 Number of Weights",
ylab = "Median Log 2 Empty Space")
detach()
ans
Figure 2.49
attach(bin.packing)
bin.packing.m <- oneway(logb(empty.space, 2) ~ number.runs, location = median, spread = 1)
srmads <- tapply(abs(residuals(bin.packing.m)), number.runs, median)
ans <- xyplot(logb(srmads, 2) ~ logb(sort(unique(number.runs)), 2),
aspect = 1,
sub = list("Figure 2.49",cex=.8),
xlab = "Log 2 Number of Weights",
ylab = "Log 2 Mad of Log 2 Empty Space")
detach()
ans
Figure 2.50
attach(bin.packing)
bin.packing.m <- oneway(logb(empty.space,2) ~ number.runs, location = median, spread = 1)
srmads <- tapply(abs(residuals(bin.packing.m)), number.runs, median)
ans <- xyplot(logb(srmads/min(srmads), 2) ~ bin.packing.m$location,
aspect = 1,
sub = list("Figure 2.50",cex=.8),
xlab = "Median Log 2 Empty Space",
ylab = "Log 2 Relative Spread")
detach()
ans
Figure 2.51
attach(bin.packing)
bin.packing.m <- oneway(empty.space~number.runs,location = median, spread=1)
srmads <- tapply(abs(residuals(bin.packing.m)),number.runs,median)
log.srmads <- logb(srmads/min(srmads),2)
ans <- xyplot(log.srmads ~ bin.packing.m$location,
aspect=1,
sub = list("Figure 2.51",cex=.8),
xlab="Median Empty Space",
ylab="Log 2 Relative Spread")
detach()
ans
Figure 2.52
res <- oneway(logb(empty.space,2) ~ number.runs, data = bin.packing,
location = median,
spread = function(x) median(abs(x-median(x))))$scaled.residuals/1.68
qqmath(~ res | factor(bin.packing$number.runs),
prepanel = prepanel.qqmathline,
panel = function(x, y) {
panel.grid()
panel.qqmathline(y, distribution = qnorm)
panel.qqmath(x, y, err=-1) # no warnings for out of bounds
},
ylim = c(-2.75, 2.75),
layout = c(3, 4),
sub = list("Figure 2.52",cex=.8),
xlab = "Unit Normal Quantile",
ylab = "Spread-Standardized Log 2 Empty Space")
UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services