RDX2 X   g9.1 source function(device=""){ - if(device!="")hardcopy(width=4, height=2, ) device=device) < oldpar <- par(mar=par()$mar-c(.5,0,1.5,0), mfrow=c(1,1))  on.exit(par(oldpar))  data(LakeHuron) D plot(LakeHuron,ylab="Depth (in feet)", xlab = "Time (in years)")  if(device!="")dev.off()  }þý device þ srcref  srcfileò encoding native.enc timestampAÒ³J)€ class POSIXt POSIXctþ filename ~/r-book/ed2/figures/figs9.R wd /Users/johnm/r/ch9þþ ÿ srcfileþ ÿ srcrefþ ))ÿÿ ÿ srcrefþ <<ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ DDÿÿ ÿ srcrefþ   ÿÿ ÿ srcrefþÿÿþ { if !=ÿ þ hardcopy width@ height@ÿÿþþ <- oldpar par mar - $ÿþÿþ c?à?øþþ mfrowÿ?ð?ðþþþ on.exitÿÿþþ data LakeHuronþ plotÿ ylab Depth (in feet) xlab Time (in years)þ ÿÿÿ þ dev.offþþþ g9.2ÿ function(device=""){ A if(device!="")hardcopy(width=5.0, height=1.65, device=device)  oldpar <- par(pty="s")  on.exit(par(oldpar))  data(LakeHuron) F lag.plot(LakeHuron, set.lags=1:4,do.lines=F, oma=c(0,1.5,2.6,1.5), ) layout=c(1,4), cex.lab=1.15) 0 mtext(side=3, line=3, "A: Lag plots", adj=0)  if(device!="")dev.off()  }þýÿ þÿ   ÿÿ ÿ srcrefþ AAÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ))ÿÿ ÿ srcrefþ 00ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ÿ?úffffffÿÿþþÿÿÿ pty sþþÿÿÿþþÿÿþ lag.plotÿ set.lags :?ð@þ do.lines F omaÿ?ø@ÌÌÌÌÌÍ?øþ layoutÿ?ð@þ cex.lab?òffffffþ mtext side@ line@ A: Lag plots adjþ ÿÿÿ þ ÿþþþ g9.2bÿ function(device=""){ ? if(device!="")hardcopy(width=3.25, height=2, device=device) 7 oldpar <- par(mar=c(4.1,4.1,3.6,1.1), mfrow=c(1,1), ' mgp=c(2.25, 0.75, 0))  on.exit(par(oldpar))  data(LakeHuron)  acf(LakeHuron, main="") D mtext(side=3, line=0.75, "B: Estimated autocorrelations", adj=0)  if(device!="")dev.off()  }þýÿ þÿ ÿÿ ÿ srcrefþ ??ÿÿ ÿ srcrefþ ''ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ   DDÿÿ ÿ srcrefþ !!ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ ÿ@ÿÿþþÿÿÿÿÿ@ffffff@ffffff@ ÌÌÌÌÌÍ?ñ™™™™™šþÿÿ?ð?ðþ mgpÿ@?èþþþÿÿÿþþÿÿþ acfÿ main þ+ÿ,ÿ@-ÿ?è B: Estimated autocorrelations.ÿþ ÿÿÿ þ ÿþþþ g9.3ÿ function( device=""){ 2 if(device!="")hardcopy(width=3.25, height=1.8, ) device=device) M oldpar <- par(mar=c(4.1,4.1,1.1,1.1), mgp=c(2.25, 0.75, 0), mfrow=c(1,1))  on.exit(par(oldpar))  data(LakeHuron) + acf(LakeHuron, type="partial", main="")  if(device!="")dev.off()  }þýÿ þÿ %%ÿÿ ÿ srcrefþ &'))ÿÿ ÿ srcrefþ ((MMÿÿ ÿ srcrefþ ))ÿÿ ÿ srcrefþ **ÿÿ ÿ srcrefþ ++++ÿÿ ÿ srcrefþ ,,ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ ÿ?üÌÌÌÌÌÍÿÿþþÿÿÿÿÿ@ffffff@ffffff?ñ™™™™™š?ñ™™™™™šþ0ÿÿ@?èþÿÿ?ð?ðþþþÿÿÿþþÿÿþ1ÿÿ type partial2ÿ þ ÿÿÿ þ ÿþþþ g9.4ÿ# function(device="", lag.max=11){ 3 if(device!="")hardcopy(width=3.75, height=2.75, F device=device, pointsize=c(7,7), trellis=T) ? trellis.par.set(superpose.line=list(lwd=c(2,1), lty=c(1,1), - col=c("black","gray")))  library(lattice) 1 four <- c(ARMAacf(ma=c(0.5),lag.max=lag.max), 6 ARMAacf(ma=c(0,0,0.5), lag.max=lag.max), ; ARMAacf(ma=c(0,0,0.5,0,0.5),lag.max=lag.max), 5 ARMAacf(ma=rep(0.5,6),lag.max=lag.max)) : pfour <- c(ARMAacf(ma=c(0.5),lag.max=lag.max, pacf=T), ? ARMAacf(ma=c(0,0,0.5), lag.max=lag.max, pacf=T), D ARMAacf(ma=c(0,0,0.5,0,0.5),lag.max=lag.max, pacf=T), > ARMAacf(ma=rep(0.5,6),lag.max=lag.max, pacf=T)) J xy <- data.frame(acf = c(four,pfour), lag = c(rep(0:lag.max, 4)-0.095, D rep(1:lag.max,4)+0.095), K gp=rep(c("acf", "pacf"), c((lag.max+1)*4, lag.max*4)),  series = ( c(rep(c("ma = 0.5", 1 "ma = c(0, 0, 0.5)", 9 "ma = c(0, 0, 0.5, 0, 0.5)", E "ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)"), - rep(lag.max+1,4)), ( rep(c("ma = 0.5", 1 "ma = c(0, 0, 0.5)", 9 "ma = c(0, 0, 0.5, 0, 0.5)", E "ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)"), , rep(lag.max,4))))  > print(xyplot(acf ~ lag|series, data = xy, layout = c(2,2), 5 scale = list(cex = 0.75), groups=gp, 3 par.strip.text = list(cex = 0.65), 7 type = "h", ylab = "Autocorrelation"))  if(device!="") dev.off()  }þýÿ  lag.max@&þÿ  0 0 ÿÿ ÿ srcrefþ 12FFÿÿ ÿ srcrefþ 34--ÿÿ ÿ srcrefþ 55ÿÿ ÿ srcrefþ 6955ÿÿ ÿ srcrefþ :=>>ÿÿ ÿ srcrefþ >K,,ÿÿ ÿ srcrefþ MP77ÿÿ ÿ srcrefþ QQÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ÿ@ÿÿ pointsizeÿ@@þ trellis Tþþ trellis.par.set superpose.line list lwdÿ@?ðþ ltyÿ?ð?ðþ colÿ black grayþþþ library latticeþÿ fourÿ ARMAacf maÿ?àþ6ÿ6ÿþCÿDÿÿ?àþ6ÿ6ÿþCÿDÿÿ?à?àþ6ÿ6ÿþCÿDÿ rep?à@þ6ÿ6ÿþþþÿ pfourÿCÿDÿÿ?àþ6ÿ6ÿ pacf9ÿþCÿDÿÿ?àþ6ÿ6ÿGÿ9ÿþCÿDÿÿ?à?àþ6ÿ6ÿGÿ9ÿþCÿDÿEÿ?à@þ6ÿ6ÿGÿ9ÿþþþÿ xy data.frame1ÿÿBÿFÿþ lagÿÿEÿ%ÿ6ÿþ@þ?¸Që…¸Rþ +Eÿ%ÿ?ð6ÿþ@þ?¸Që…¸Rþþ gpEÿÿ acf pacfþÿ * (Kÿ6ÿ?ðþþ@þMÿ6ÿ@þþþ seriesÿEÿÿ ma = 0.5 ma = c(0, 0, 0.5) ma = c(0, 0, 0.5, 0, 0.5) $ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)þEÿKÿ6ÿ?ðþ@þþEÿÿ ma = 0.5 ma = c(0, 0, 0.5) ma = c(0, 0, 0.5, 0, 0.5) $ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)þEÿ6ÿ@þþþþþ print xyplot ~1ÿ |JÿOÿþþÿHÿ)ÿÿ@@þ scale<ÿ cex?èþ groupsLÿ par.strip.text<ÿUÿ?äÌÌÌÌÌÍþ4ÿ hÿ Autocorrelationþþ ÿÿÿ þ ÿþþþ g9.5ÿ function(device=""){ B if(device!="")hardcopy(width=3.25, height=3.25, device=device) ? oldpar <- par(mar=c(4.1,4.1,1.1,1.1), mgp=c(2.75, 0.75,0))  on.exit(par(oldpar))  if(!exists("xbomsoi")){ J xbomsoi <- with(bomsoi, data.frame(SOI=SOI, cuberootRain=avrain^0.33)) - xbomsoi$trendSOI <- lowess(xbomsoi$SOI)$y 7 xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y E with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) " xbomsoi$detrendRain <- I with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) ! xbomsoi$detrendSOI <- : with(xbomsoi, SOI - trendSOI + mean(trendSOI)) C xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y  } ? ## Plot time series avrain and SOI: ts object bomsoi (DAAG) = plot(ts(xbomsoi[, c("cuberootRain","SOI")], start=1900), = panel=function(y,...)panel.smooth(1900:2004, y,...), B xlab = "Year", main="", ylim=list(c(250, 800),c(-20,25)))  if(device!="") dev.off()  invisible(xbomsoi)  }þýÿ þÿ UUÿÿ ÿ srcrefþ VVBBÿÿ ÿ srcrefþ WW??ÿÿ ÿ srcrefþ XXÿÿ ÿ srcrefþ Ycÿÿ ÿ srcrefþ egBBÿÿ ÿ srcrefþ hhÿÿ ÿ srcrefþ iiÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ ÿ@ ÿÿþþÿÿÿÿÿ@ffffff@ffffff?ñ™™™™™š?ñ™™™™™šþ0ÿÿ@?èþþþÿÿÿþþ ÿ ! exists xbomsoiþþÿ YYÿÿ ÿ srcrefþ ZZJJÿÿ ÿ srcrefþ [[--ÿÿ ÿ srcrefþ \\77ÿÿ ÿ srcrefþ ] ]E Eÿÿ ÿ srcrefþ ^ _E Eÿÿ ÿ srcrefþ ` a: :ÿÿ ÿ srcrefþ b b? ?ÿÿ ÿ srcrefþÿÿþ ÿÿ xbomsoi with bomsoiIÿ SOI^ÿ cuberootRain ^ avrain?Õ¸Që…þþþþÿÿ[ÿ trendSOIþÿ lowessÿ[ÿ^ÿþþ yþþÿÿ[ÿ trendRainþÿcÿÿ[ÿ_ÿþþdÿþþ\ÿ[ÿKÿÿ_ÿeÿþ meaneÿþþþÿÿ[ÿ detrendRainþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþþÿÿ[ÿ detrendSOIþ\ÿ[ÿKÿÿ^ÿbÿþfÿbÿþþþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþþþÿ ts [[ÿûÿ cuberootRain SOIþþ start@°þ panel functiondÿû ...ûþ panel.smooth%ÿ@°@ŸPþdÿnÿþ -function(y,...)panel.smooth(1900:2004, y,...)þÿ Year2ÿ  ylim<ÿÿ@o@@‰þÿÿ@4þ@9þþþ ÿÿÿ þ ÿþþ invisible[ÿþþ g9.6ÿ function(device=""){ @ if(device!="")hardcopy(width=2.4, height=2.4, device=device) > oldpar <- par(mar=c(4.1,4.1,1.1,1.1), mgp=c(2.75,0.75,0))  on.exit(par(oldpar))  if(!exists("xbomsoi")){ J xbomsoi <- with(bomsoi, data.frame(SOI=SOI, cuberootRain=avrain^0.33)) - xbomsoi$trendSOI <- lowess(xbomsoi$SOI)$y 7 xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y E with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) " xbomsoi$detrendRain <- I with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) ! xbomsoi$detrendSOI <- : with(xbomsoi, SOI - trendSOI + mean(trendSOI)) C xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y  } ' rainpos <- pretty(bomsoi$avrain, 5)  with(xbomsoi, 0 {plot(cuberootRain ~ SOI, xlab = "SOI", = ylab = "Rainfall (cube root scale)", yaxt="n") ; axis(2, at = rainpos^0.33, labels=paste(rainpos)) 5 ## Relative changes in the two trend curves + lines(lowess(cuberootRain ~ SOI)) E lines(lowess(trendRain ~ trendSOI), lwd=1.75, col="gray40")  })  if(device!="") dev.off()  invisible(xbomsoi)  }þýÿ þÿ  mmÿÿ ÿ srcrefþ nn@@ÿÿ ÿ srcrefþ oo>>ÿÿ ÿ srcrefþ ppÿÿ ÿ srcrefþ q{ÿÿ ÿ srcrefþ ||''ÿÿ ÿ srcrefþ }„  ÿÿ ÿ srcrefþ ……ÿÿ ÿ srcrefþ ††ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@333333ÿ@333333ÿÿþþÿÿÿÿÿ@ffffff@ffffff?ñ™™™™™š?ñ™™™™™šþ0ÿÿ@?èþþþÿÿÿþþ ÿYÿZÿ xbomsoiþþÿ qqÿÿ ÿ srcrefþ rrJJÿÿ ÿ srcrefþ ss--ÿÿ ÿ srcrefþ tt77ÿÿ ÿ srcrefþ u uE Eÿÿ ÿ srcrefþ v wE Eÿÿ ÿ srcrefþ x y: :ÿÿ ÿ srcrefþ z z? ?ÿÿ ÿ srcrefþÿÿþ ÿÿ[ÿ\ÿ]ÿIÿ^ÿ^ÿ_ÿ`ÿaÿ?Õ¸Që…þþþþÿÿ[ÿbÿþÿcÿÿ[ÿ^ÿþþdÿþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþÿÿ[ÿgÿþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþþÿÿ[ÿhÿþ\ÿ[ÿKÿÿ^ÿbÿþfÿbÿþþþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþþþÿ rainpos prettyÿ]ÿaÿþ@þþ\ÿ[ÿÿ ~ ~ ÿÿ ÿ srcrefþ ~ = =ÿÿ ÿ srcrefþ € €; ;ÿÿ ÿ srcrefþ ‚ ‚+ +ÿÿ ÿ srcrefþ ƒ ƒE Eÿÿ ÿ srcrefþÿÿþ ÿÿRÿ_ÿ^ÿþÿ SOIÿ Rainfall (cube root scale) yaxt nþ axis@ at`ÿsÿ?Õ¸Që…þ labels pastesÿþþ linescÿRÿ_ÿ^ÿþþþzÿcÿRÿeÿbÿþþ=ÿ?ü?ÿ gray40þþþ ÿÿÿ þ ÿþþqÿ[ÿþþ g9.7ÿ function(device=""){ ? if(device!="")hardcopy(width=4.25, height=2, device=device) M oldpar <- par(mar=c(4.1,4.6,2.1,1.1), mfrow = c(1,2), mgp=c(2.5,0.75,0), ' oma = c(0, 0, 0, 0))  on.exit(par(oldpar))  if(!exists("xbomsoi")){  xbomsoi <- G with(bomsoi, data.frame(SOI=SOI, cuberootRain=avrain^0.33)) 5 xbomsoi$trendSOI <- lowess(xbomsoi$SOI)$y " xbomsoi$detrendRain <- E with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) ! xbomsoi$detrendSOI <- : with(xbomsoi, SOI - trendSOI + mean(trendSOI)) ? xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y } 0 rainpos <- pretty(xbomsoi$cuberootRain^3, 6) - plot(cuberootRain ~ SOI, data = xbomsoi, 7 ylab = "Rainfall (cube root scale)", yaxt="n") 5 axis(2, at = rainpos^0.33, labels=paste(rainpos)) . mtext(side = 3, line = 0.25, "A", adj = 0) 4 with(xbomsoi, lines(lowess(cuberootRain ~ SOI))) 2 plot(detrendRain ~ detrendSOI, data = xbomsoi, E xlab="Detrended SOI", ylab = "Detrended rainfall", yaxt="n") 6 axis(2, at = rainpos^0.33, labels=paste(rainpos)) : with(xbomsoi, lines(lowess(detrendRain ~ detrendSOI))) . mtext(side = 3, line = 0.25, "B", adj = 0)  if(device!="") dev.off()  }þýÿ þÿ ŠŠÿÿ ÿ srcrefþ ‹‹??ÿÿ ÿ srcrefþ Œ''ÿÿ ÿ srcrefþ ŽŽÿÿ ÿ srcrefþ ˜  ÿÿ ÿ srcrefþ ™™00ÿÿ ÿ srcrefþ š›77ÿÿ ÿ srcrefþ œœ55ÿÿ ÿ srcrefþ ..ÿÿ ÿ srcrefþ žž44ÿÿ ÿ srcrefþ Ÿ EEÿÿ ÿ srcrefþ ¡¡55ÿÿ ÿ srcrefþ ¢¢::ÿÿ ÿ srcrefþ ££..ÿÿ ÿ srcrefþ ¤¤ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ÿ@ÿÿþþÿÿÿÿÿ@ffffff@ffffff@ÌÌÌÌÌÍ?ñ™™™™™šþÿÿ?ð@þ0ÿÿ@?èþ(ÿÿþþþÿÿÿþþ ÿYÿZÿ xbomsoiþþÿ ÿÿ ÿ srcrefþ  ‘G Gÿÿ ÿ srcrefþ ’ ’5 5ÿÿ ÿ srcrefþ “ ”E Eÿÿ ÿ srcrefþ • –: :ÿÿ ÿ srcrefþ — —? ?ÿÿ ÿ srcrefþÿÿþ ÿÿ[ÿ\ÿ]ÿIÿ^ÿ^ÿ_ÿ`ÿaÿ?Õ¸Që…þþþþÿÿ[ÿbÿþÿcÿÿ[ÿ^ÿþþdÿþþÿÿ[ÿgÿþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþþÿÿ[ÿhÿþ\ÿ[ÿKÿÿ^ÿbÿþfÿbÿþþþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþþþÿsÿtÿ`ÿÿ[ÿ_ÿþ@þ@þþÿRÿ_ÿ^ÿþÿ[ÿÿ Rainfall (cube root scale)uÿ nþvÿ@wÿ`ÿsÿ?Õ¸Që…þxÿyÿsÿþþ+ÿ,ÿ@-ÿ?Ð A.ÿþ\ÿ[ÿzÿcÿRÿ_ÿ^ÿþþþþÿRÿgÿhÿþÿ[ÿÿ Detrended SOIÿ Detrended rainfalluÿ nþvÿ@wÿ`ÿsÿ?Õ¸Që…þxÿyÿsÿþþ\ÿ[ÿzÿcÿRÿgÿhÿþþþþ+ÿ,ÿ@-ÿ?Ð B.ÿþ ÿÿÿ þ ÿþþþ g9.8ÿ function(device=""){ A if(device!="")hardcopy(width=3.25, height=1.8, device=device) M oldpar <- par(mar=c(4.1,4.1,1.1,1.1), mfrow=c(1,1), mgp=c(2.25, 0.75, 0))  on.exit(par(oldpar))  if(!exists("xbomsoi")){ J xbomsoi <- with(bomsoi, data.frame(SOI=SOI, cuberootRain=avrain^0.33)) - xbomsoi$trendSOI <- lowess(xbomsoi$SOI)$y 7 xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y E with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) " xbomsoi$detrendRain <- I with(xbomsoi, cuberootRain - trendRain + mean(trendRain)) ! xbomsoi$detrendSOI <- : with(xbomsoi, SOI - trendSOI + mean(trendSOI)) C xbomsoi$trendRain <- lowess(xbomsoi$cuberootRain)$y  } E acf(resid(lm(detrendRain ~ detrendSOI, data = xbomsoi)), main="")  if(device!="")dev.off()  }þýÿ þÿ ¨¨ÿÿ ÿ srcrefþ ©©AAÿÿ ÿ srcrefþ ªªMMÿÿ ÿ srcrefþ ««ÿÿ ÿ srcrefþ ¬¶ÿÿ ÿ srcrefþ ··EEÿÿ ÿ srcrefþ ¸¸ÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ ÿ?üÌÌÌÌÌÍÿÿþþÿÿÿÿÿ@ffffff@ffffff?ñ™™™™™š?ñ™™™™™šþÿÿ?ð?ðþ0ÿÿ@?èþþþÿÿÿþþ ÿYÿZÿ xbomsoiþþÿ ¬¬ÿÿ ÿ srcrefþ ­­JJÿÿ ÿ srcrefþ ®®--ÿÿ ÿ srcrefþ ¯¯77ÿÿ ÿ srcrefþ ° °E Eÿÿ ÿ srcrefþ ± ²E Eÿÿ ÿ srcrefþ ³ ´: :ÿÿ ÿ srcrefþ µ µ? ?ÿÿ ÿ srcrefþÿÿþ ÿÿ[ÿ\ÿ]ÿIÿ^ÿ^ÿ_ÿ`ÿaÿ?Õ¸Që…þþþþÿÿ[ÿbÿþÿcÿÿ[ÿ^ÿþþdÿþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþÿÿ[ÿgÿþ\ÿ[ÿKÿÿ_ÿeÿþfÿeÿþþþþÿÿ[ÿhÿþ\ÿ[ÿKÿÿ^ÿbÿþfÿbÿþþþþÿÿ[ÿeÿþÿcÿÿ[ÿ_ÿþþdÿþþþþ1ÿ resid lmRÿgÿhÿþÿ[ÿþþ2ÿ þ ÿÿÿ þ ÿþþþ g9.9ÿ$ function(device="", lag.max=11){ 3 if(device!="")hardcopy(width=3.75, height=2.75, F device=device, pointsize=c(7,7), trellis=T) ? trellis.par.set(superpose.line=list(lwd=c(2,1), lty=c(1,1), - col=c("black","gray")))  on.exit(par(oldpar))  library(lattice) 1 four <- c(ARMAacf(ma=c(0.5),lag.max=lag.max), 6 ARMAacf(ma=c(0,0,0.5), lag.max=lag.max), ; ARMAacf(ma=c(0,0,0.5,0,0.5),lag.max=lag.max), 5 ARMAacf(ma=rep(0.5,6),lag.max=lag.max)) : pfour <- c(ARMAacf(ma=c(0.5),lag.max=lag.max, pacf=T), ? ARMAacf(ma=c(0,0,0.5), lag.max=lag.max, pacf=T), D ARMAacf(ma=c(0,0,0.5,0,0.5),lag.max=lag.max, pacf=T), > ARMAacf(ma=rep(0.5,6),lag.max=lag.max, pacf=T)) J xy <- data.frame(acf = c(four,pfour), lag = c(rep(0:lag.max, 4)-0.095, D rep(1:lag.max,4)+0.095), K gp=rep(c("acf", "pacf"), c((lag.max+1)*4, lag.max*4)),  series = ( c(rep(c("ma = 0.5", 1 "ma = c(0, 0, 0.5)", 9 "ma = c(0, 0, 0.5, 0, 0.5)", E "ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)"), - rep(lag.max+1,4)), ( rep(c("ma = 0.5", 1 "ma = c(0, 0, 0.5)", 9 "ma = c(0, 0, 0.5, 0, 0.5)", E "ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)"), , rep(lag.max,4))))  > print(xyplot(acf ~ lag|series, data = xy, layout = c(2,2), 5 scale = list(cex = 0.75), groups=gp, 3 par.strip.text = list(cex = 0.65), 7 type = "h", ylab = "Autocorrelation"))  if(device!="") dev.off()  }þýÿ 6ÿ@&þÿ  ¼ ¼ ÿÿ ÿ srcrefþ ½¾FFÿÿ ÿ srcrefþ ¿À--ÿÿ ÿ srcrefþ ÁÁÿÿ ÿ srcrefþ ÂÂÿÿ ÿ srcrefþ ÃÆ55ÿÿ ÿ srcrefþ ÇÊ>>ÿÿ ÿ srcrefþ ËØ,,ÿÿ ÿ srcrefþ ÚÝ77ÿÿ ÿ srcrefþ ÞÞÿÿ ÿ srcrefþÿÿþ ÿ ÿÿÿ þÿÿ@ÿ@ÿÿ7ÿÿ@@þ8ÿ9ÿþþ:ÿ;ÿ<ÿ=ÿÿ@?ðþ>ÿÿ?ð?ðþ?ÿÿ black grayþþþÿÿÿþþ@ÿAÿþÿBÿÿCÿDÿÿ?àþ6ÿ6ÿþCÿDÿÿ?àþ6ÿ6ÿþCÿDÿÿ?à?àþ6ÿ6ÿþCÿDÿEÿ?à@þ6ÿ6ÿþþþÿFÿÿCÿDÿÿ?àþ6ÿ6ÿGÿ9ÿþCÿDÿÿ?àþ6ÿ6ÿGÿ9ÿþCÿDÿÿ?à?àþ6ÿ6ÿGÿ9ÿþCÿDÿEÿ?à@þ6ÿ6ÿGÿ9ÿþþþÿHÿIÿ1ÿÿBÿFÿþJÿÿÿEÿ%ÿ6ÿþ@þ?¸Që…¸RþKÿEÿ%ÿ?ð6ÿþ@þ?¸Që…¸RþþLÿEÿÿ acf pacfþÿMÿNÿKÿ6ÿ?ðþþ@þMÿ6ÿ@þþþOÿÿEÿÿ ma = 0.5 ma = c(0, 0, 0.5) ma = c(0, 0, 0.5, 0, 0.5) $ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)þEÿKÿ6ÿ?ðþ@þþEÿÿ ma = 0.5 ma = c(0, 0, 0.5) ma = c(0, 0, 0.5, 0, 0.5) $ma = c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5)þEÿ6ÿ@þþþþþPÿQÿRÿ1ÿSÿJÿOÿþþÿHÿ)ÿÿ@@þTÿ<ÿUÿ?èþVÿLÿWÿ<ÿUÿ?äÌÌÌÌÌÍþ4ÿ hÿ Autocorrelationþþ ÿÿÿ þ ÿþþþ gdumpÿ 7function(fnam=NULL, prefix="~/r-book/ed2/figures/figs",  splitchar="/ch"){  if(is.null(fnam)){  path <- getwd() 7 pathtag <- strsplit(path, "/ch", fixed=TRUE)[[1]] 5 fnam <- paste(prefix, pathtag[length(pathtag)], ! ".R", sep="")  } - else fnam <- paste(prefix, fnam, sep="/") H objnames <- c(objects(pattern="^g", envir=sys.frame(0)), "hardcopy") & cat("\nDump to file:", fnam, "\n")  print(objnames)  dump(objnames, fnam)  }þý fnamþ prefix ~/r-book/ed2/figures/figs splitchar /chþÿ ããÿÿ ÿ srcrefþ äê--ÿÿ ÿ srcrefþ ëëHHÿÿ ÿ srcrefþ ìì&&ÿÿ ÿ srcrefþ ííÿÿ ÿ srcrefþ îîÿÿ ÿ srcrefþÿÿþ ÿ ÿ is.nullÿþÿ ääÿÿ ÿ srcrefþ ååÿÿ ÿ srcrefþ ææ77ÿÿ ÿ srcrefþ çè!!ÿÿ ÿ srcrefþÿÿþ ÿÿ path getwdþþÿ pathtag [[ strsplit…ÿ /ch fixed þ?ðþþÿÿyÿ‚ÿjÿ‡ÿ length‡ÿþþ .R sep þþþÿÿyÿ‚ÿÿŒÿ /þþþÿ objnamesÿ objects pattern ^g envir sys.frameþþ hardcopyþþ cat  Dump to file:ÿ  þPÿÿþ dumpÿÿþþ gsaveÿ 7function(fnam=NULL, prefix="~/r-book/ed2/figures/figs", L splitchar="/ch", xtras=c("hardcopy", "renum.fun","renum.files")){  if(is.null(fnam)){  path <- getwd() 7 pathtag <- strsplit(path, "/ch", fixed=TRUE)[[1]] 5 fnam <- paste(prefix, pathtag[length(pathtag)], % ".RData", sep="")  } - else fnam <- paste(prefix, fnam, sep="/") C objnames <- c(objects(pattern="^g", envir=sys.frame(0)), xtras) & cat("\nDump to file:", fnam, "\n")  print(objnames) " save(list=objnames, file=fnam)  }þýÿþ‚ÿ ~/r-book/ed2/figures/figsƒÿ /ch xtrasÿ hardcopy renum.fun renum.filesþþÿ óLóLLLÿÿ ÿ srcrefþ ôú--ÿÿ ÿ srcrefþ ûûCCÿÿ ÿ srcrefþ üü&&ÿÿ ÿ srcrefþ ýýÿÿ ÿ srcrefþ þþ""ÿÿ ÿ srcrefþÿÿþ ÿ ÿ„ÿÿþÿ ôôÿÿ ÿ srcrefþ õõÿÿ ÿ srcrefþ öö77ÿÿ ÿ srcrefþ ÷ø%%ÿÿ ÿ srcrefþÿÿþ ÿÿ…ÿ†ÿþþÿ‡ÿˆÿ‰ÿ…ÿ /chŠÿ þ?ðþþÿÿyÿ‚ÿjÿ‡ÿ‹ÿ‡ÿþþ .RDataŒÿ þþþÿÿyÿ‚ÿÿŒÿ /þþþÿÿÿŽÿÿ ^gÿ‘ÿþþ•ÿþþ’ÿ  Dump to file:ÿ  þPÿÿþ save<ÿÿ fileÿþþÿÿ4 =function(width=3.75, height=3.75, color=FALSE, trellis=FALSE, H device=c("","pdf","ps"), path="~/r-book/ed2/Art/", file=NULL, D format=c("nn-nn", "name"), split="\\.", pointsize=c(8,4),  fonts=NULL,  horiz=FALSE, ...){ ) if(!trellis)pointsize <- pointsize[1]  funtxt <- sys.call(1) B nam <- strsplit(as.character(funtxt), "(", fixed=TRUE)[[1]][1] 3 suffix <- switch(device, ps=".eps", pdf=".pdf") N if(is.character(path) & nchar(path)>1 & substring(path, nchar(path))!="/") & path <- paste(path, "/", sep="") - if(is.null(file)) if(format[1]=="nn-nn"){ C if(!is.null(split))dotsplit <- strsplit(nam, split)[[1]] else  dotsplit <- nam 8 if(length(dotsplit)==1)dotsplit <- c("", dotsplit) E nn2 <- paste(if(nchar(dotsplit[2])==1)"0" else "", dotsplit[2],  sep="")  if(nchar(dotsplit[1])>0){ O numstart <- which(unlist(strsplit(dotsplit[1], "")) %in% paste(0:9))[1] / nn1 <- substring(dotsplit[1], numstart) G nn1 <- paste(if(nchar(nn1) == 1) "0" else "", nn1, "-", sep="")  } else nn1 <- "" % file <- paste(nn1, nn2, sep="")  } else file <- nam L if(nchar(file)>4 & substring(file, nchar(file)-nchar(suffix)+1)==suffix)  suffix <- "" - file <- paste(path, file, suffix, sep="") : print(paste("Output will be directed to file:", file))  dev.out <- device[1] 6 dev.fun <- switch(dev.out, pdf=pdf, ps=postscript)  if(trellis){  library(lattice)  if(device=="ps") 1 trellis.device(file=file, device=dev.fun, ? color = color, horiz=horiz, fonts=fonts, < width=width, height=height, ...) else < trellis.device(file=file, device=dev.fun, fonts=fonts, D color = color, width=width, height=height, ...) R trellis.par.set(list(fontsize=list(text=pointsize[1], points=pointsize[2]))) } else  if (dev.out!=""){  print(c(width, height))  if(device=="ps") F dev.fun(file=file, paper="special", horiz=horiz, fonts=fonts, M width=width, height=height, pointsize=pointsize[1], ...) else 6 dev.fun(file=file, paper="special", fonts=fonts, L width=width, height=height, pointsize=pointsize[1], ...)  } D if(trellis)trellis.par.set(list(fontsize=list(text=pointsize[1], < points=pointsize[2])))  }þýÿ@ÿ@ color 8ÿ ÿÿ  pdf psþ…ÿ ~/r-book/ed2/Art/—ÿþ formatÿ nn-nn nameþ split \.7ÿÿ@ @þ fontsþ horiz nÿûþÿ ÿÿ ÿ srcrefþ ))ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ   BBÿÿ ÿ srcrefþ   33ÿÿ ÿ srcrefþ   &&ÿÿ ÿ srcrefþ  ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ --ÿÿ ÿ srcrefþ ::ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ 66ÿÿ ÿ srcrefþ  2ÿÿ ÿ srcrefþ 34<<ÿÿ ÿ srcrefþÿÿþ ÿ ÿYÿ8ÿþÿ7ÿjÿ7ÿ?ðþþþÿ funtxt sys.call?ðþþÿ namjÿˆÿ‰ÿ as.characterÿþ (Šÿ þ?ðþ?ðþþÿ suffix switchÿ ps .eps pdf .pdfþþ ÿ &¥ÿ is.character…ÿþ > nchar…ÿþ?ðþþÿ substring…ÿ¨ÿ…ÿþþ /þþÿ…ÿyÿ…ÿ /Œÿ þþþ ÿ„ÿ—ÿþ ÿ ==jÿ™ÿ?ðþ nn-nnþÿ  - ---ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ 88ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ ÿÿ ÿ srcrefþ %%ÿÿ ÿ srcrefþÿÿþ ÿ ÿYÿ„ÿšÿþþÿ dotsplitˆÿ‰ÿŸÿšÿþ?ðþþÿ«ÿŸÿþþ ÿªÿ‹ÿ«ÿþ?ðþÿ«ÿÿ «ÿþþþÿ nn2yÿ ÿªÿ¨ÿjÿ«ÿ@þþ?ðþ 0 þjÿ«ÿ@þŒÿ þþ ÿ§ÿ¨ÿjÿ«ÿ?ðþþþÿ ÿÿ ÿ srcrefþ  O Oÿÿ ÿ srcrefþ  / /ÿÿ ÿ srcrefþ  G Gÿÿ ÿ srcrefþÿÿþ ÿÿ numstartjÿ which %in% unlist‰ÿjÿ«ÿ?ðþ þþyÿ%ÿ@"þþþþ?ðþþÿ nn1©ÿjÿ«ÿ?ðþ­ÿþþÿ±ÿyÿ ÿªÿ¨ÿ±ÿþ?ðþ 0 þ±ÿ -Œÿ þþþÿ±ÿ þþÿ—ÿyÿ±ÿ¬ÿŒÿ þþþÿ—ÿŸÿþþþ ÿ¥ÿ§ÿ¨ÿ—ÿþ@þªÿ©ÿ—ÿKÿÿ¨ÿ—ÿþ¨ÿ¡ÿþþ?ðþþ¡ÿþþÿ¡ÿ þþÿ—ÿyÿ…ÿ—ÿ¡ÿŒÿ þþPÿyÿ Output will be directed to file:—ÿþþÿ dev.outjÿÿ?ðþþÿ dev.fun¢ÿ²ÿ¤ÿ¤ÿ£ÿ postscriptþþ ÿ8ÿÿ   ÿÿ ÿ srcrefþ !!ÿÿ ÿ srcrefþ "'DDÿÿ ÿ srcrefþ ))RRÿÿ ÿ srcrefþÿÿþ ÿ@ÿAÿþ ÿªÿÿ psþ trellis.device—ÿ—ÿÿ³ÿ˜ÿ˜ÿœÿœÿ›ÿ›ÿÿÿÿÿnÿþµÿ—ÿ—ÿÿ³ÿ›ÿ›ÿ˜ÿ˜ÿÿÿÿÿnÿþþ:ÿ<ÿ fontsize<ÿ textjÿ7ÿ?ðþ pointsjÿ7ÿ@þþþþþ ÿÿ²ÿ þÿ ++ÿÿ ÿ srcrefþ ,,ÿÿ ÿ srcrefþ -1FFÿÿ ÿ srcrefþÿÿþ ÿPÿÿÿÿþþ ÿªÿÿ psþ³ÿ—ÿ—ÿ paper specialœÿœÿ›ÿ›ÿÿÿÿÿ7ÿjÿ7ÿ?ðþnÿþ³ÿ—ÿ—ÿ¹ÿ special›ÿ›ÿÿÿÿÿ7ÿjÿ7ÿ?ðþnÿþþþþþ ÿ8ÿ:ÿ<ÿ¶ÿ<ÿ·ÿjÿ7ÿ?ðþ¸ÿjÿ7ÿ@þþþþþþ renum.funÿ Cfunction(from.prefix=NULL, to.prefix="g",from=4:1, to=4:1, doit=F){  path <- getwd() 5 pathtag <- strsplit(path, "/ch", fixed=TRUE)[[1]] & endbit <- pathtag[length(pathtag)] 1 to.prefix <- paste(to.prefix, endbit, sep="") 4 if(is.null(from.prefix))from.prefix <- to.prefix  for(i in 1:length(to)) M {txt<-paste(to.prefix,".",to[i]," <- ", from.prefix,".",from[i],sep="") 7 if(doit)eval(parse(text=txt),envir=sys.frame(0))  print(txt) " if(from.prefix!=to.prefix){ B rm.txt <- paste("rm(",from.prefix,".",from[i],")",sep="") < if(doit)eval(parse(text=rm.txt),envir=sys.frame(0))  print(rm.txt)  }  }  }þý from.prefixþ to.prefix g$ from%ÿ@?ðþ$ to%ÿ@?ðþ$ doit'ÿþ ÿÿ…ÿ†ÿþþÿ‡ÿˆÿ‰ÿ…ÿ /chŠÿ þ?ðþþÿ endbitjÿ‡ÿ‹ÿ‡ÿþþþÿ¼ÿyÿ¼ÿÀÿŒÿ þþ ÿ„ÿ»ÿþÿ»ÿ¼ÿþþ for i%ÿ?ð‹ÿ¾ÿþþ ÿÿ txtyÿ¼ÿ .jÿ¾ÿÂÿþ  <- »ÿ .jÿ½ÿÂÿþŒÿ þþ ÿ¿ÿ eval parse·ÿÃÿþÿ‘ÿþþþPÿÃÿþ ÿÿ»ÿ¼ÿþ ÿÿ rm.txtyÿ rm(»ÿ .jÿ½ÿÂÿþ )Œÿ þþ ÿ¿ÿÄÿÅÿ·ÿÆÿþÿ‘ÿþþþPÿÆÿþþþþþþ renum.filesÿ Nfunction(from.prefix="~/r-book/second/Art/", to.prefix="~/r-book/second/Art/", ' from=20:7, to=21:8, doit=F){  path <- getwd() 5 pathtag <- strsplit(path, "/ch", fixed=TRUE)[[1]] & endbit <- pathtag[length(pathtag)] 2 if(nchar(endbit)==2)chap <- paste(endbit) else ' chap <- paste("0",endbit,sep="") 8 from.prefix <- paste(from.prefix, chap, "-", sep="") 4 to.prefix <- paste(to.prefix, chap, "-", sep="")  for(i in 1:length(from)){ 8 if (from[i]<=9) ltext <- paste("0",from[i],sep="") " else ltext <- paste(from[i]) 4 if (to[i]<=9) rtext <- paste("0",to[i],sep="") , else rtext <- paste(to[i]) 8 txt<-paste("mv ", from.prefix, ltext, ".eps", " ", 2 to.prefix, rtext, ".eps", sep="") ; backup<-paste("cp ", from.prefix, ltext, ".eps", " ", # "archive", sep="")  if(doit)system(backup)  if(doit)system(txt)  print(backup)  print(txt)  }  }þý»ÿ ~/r-book/second/Art/¼ÿ ~/r-book/second/Art/½ÿ%ÿ@4@þ¾ÿ%ÿ@5@ þ¿ÿ'ÿþ ÿÿ…ÿ†ÿþþÿ‡ÿˆÿ‰ÿ…ÿ /chŠÿ þ?ðþþÿÀÿjÿ‡ÿ‹ÿ‡ÿþþþ ÿªÿ¨ÿÀÿþ@þÿ chapyÿÀÿþþÿÈÿyÿ 0ÀÿŒÿ þþþÿ»ÿyÿ»ÿÈÿ -Œÿ þþÿ¼ÿyÿ¼ÿÈÿ -Œÿ þþÁÿÂÿ%ÿ?ð‹ÿ½ÿþþ ÿ ÿ <=jÿ½ÿÂÿþ@"þÿ ltextyÿ 0jÿ½ÿÂÿþŒÿ þþÿÊÿyÿjÿ½ÿÂÿþþþþ ÿÉÿjÿ¾ÿÂÿþ@"þÿ rtextyÿ 0jÿ¾ÿÂÿþŒÿ þþÿËÿyÿjÿ¾ÿÂÿþþþþÿÃÿyÿ mv »ÿÊÿ .eps  ¼ÿËÿ .epsŒÿ þþÿ backupyÿ cp »ÿÊÿ .eps   archiveŒÿ þþ ÿ¿ÿ systemÌÿþþ ÿ¿ÿÍÿÃÿþþPÿÌÿþPÿÃÿþþþþþ