[Dplr-commits] r1084 - pkg/dplR/vignettes

[noreply at r-forge.r-project.org]

Author: andybunn
Date: 2018-02-01 20:22:45 +0100 (Thu, 01 Feb 2018)
New Revision: 1084

Modified:
   pkg/dplR/vignettes/timeseries-dplR.Rnw
Log:
typos

Modified: pkg/dplR/vignettes/timeseries-dplR.Rnw
===================================================================
--- pkg/dplR/vignettes/timeseries-dplR.Rnw	2018-01-31 09:26:45 UTC (rev 1083)
+++ pkg/dplR/vignettes/timeseries-dplR.Rnw	2018-02-01 19:22:45 UTC (rev 1084)
@@ -175,14 +175,14 @@
 @
 Instead of an AR(4) model, \code{auto.arima} went for an 
 ARMA(1,1) model -- or an ARIMA(1,0,1). The parsimony principle certainly likes 
-a nice simple ARMA(1,1) model! Note that we could look at the residuals
+a nice simple ARMA(1,1) model. Note that we could look at the residuals
 (just the first few), model coefficients, etc. quite easily. And indeed the 
 residuals are quite clean as we would expect.
 \else% If "forecast" is not available
 An example was dropped because \code{"forecast"} is not available.
 \fi% End of conditional
 \section{Frequency Domain}
-There is, at times, and almost manic desire to better characterize the 
+There is, at times, an almost manic desire to better characterize the 
 spectral aspects of a tree-ring series. In dplR, we've implemented two 
 of the most common ways that dendrochronologists go about this and there 
 are a host of other approaches in R that we won't get to in this vignette.
@@ -200,7 +200,7 @@
      type = "n", ylab = "Spectrum (dB)", xlab = "Frequency (1/yr)",
      axes = FALSE)
 grid()
-lines(redf.dat[["freq"]], redf.dat[["gxxc"]], col = "#1B9E77")
+lines(redf.dat[["freq"]], redf.dat[["gxxc"]], col = "#1B9E77",lwd=2)
 lines(redf.dat[["freq"]], redf.dat[["ci99"]], col = "#D95F02")
 lines(redf.dat[["freq"]], redf.dat[["ci95"]], col = "#7570B3")
 lines(redf.dat[["freq"]], redf.dat[["ci90"]], col = "#E7298A")
@@ -239,7 +239,7 @@
 yrs <- as.numeric(rownames(co021.crn))
 out.wave <- morlet(y1 = dat, x1 = yrs, p2 = 8, dj = 0.1,
                    siglvl = 0.99)
-wavelet.plot(out.wave, useRaster=NA)
+wavelet.plot(out.wave, useRaster=NA, reverse.y = TRUE)
 @
 \begin{figure}[ht]
   \centering
@@ -281,8 +281,10 @@
   axis(side=1)
   mtext("Years",side=1,line = 1.25)
   Offset <- 0
+  dat.mra2 <- scale(as.data.frame(dat.mra))
   for(i in nPwrs2:1){
-    x <- scale(dat.mra[[i]]) + Offset
+#    x <- scale(dat.mra[[i]]) + Offset
+    x <- dat.mra2[,i] + Offset
     lines(yrs,x)
     abline(h=Offset,lty="dashed")
     mtext(names(dat.mra)[[i]],side=2,at=Offset,line = 0)