Nothing
## ----sel0, eval=TRUE, echo=FALSE, message=FALSE--------------------------
library(WiSEBoot)
## ----sel1, eval=TRUE, echo=TRUE------------------------------------------
data(SimulatedSmoothSeries)
dim(SimulatedSmoothSeries)
SimulatedSmoothSeries[1:3, ]
## ----sel2, eval=TRUE, echo=FALSE, out.width='12cm', out.height='8cm'-----
plot(seq(1, 2^10), SimulatedSmoothSeries[,4], ty="l",
main=bquote(paste("Smooth Series with Threshold ", J[0],"=3")),
xlab="Time", ylab="Observation")
## ----sel3, eval=TRUE, echo=FALSE-----------------------------------------
data(SimulatedSNR5Series); data(SimulatedSNR9Series); data(SimulatedSNR15Series); data(SimulatedSNR25Series)
par(mfrow=c(2,2))
plot(seq(1, 2^10), SimulatedSNR5Series[,4], ty="l",
main=bquote(paste("SNR=5 Data, Threshold ", J[0],"=3")),
xlab="Time", ylab="Observation", col="gray", ylim=c(-.5, .6))
lines(seq(1, 2^10), SimulatedSmoothSeries[,4], col="red", lwd=2)
plot(seq(1, 2^10), SimulatedSNR9Series[,4], ty="l",
main=bquote(paste("SNR=9 Data, Threshold ", J[0],"=3")),
xlab="Time", ylab="Observation", col="gray", ylim=c(-.5, .6))
lines(seq(1, 2^10), SimulatedSmoothSeries[,4], col="red", lwd=2)
plot(seq(1, 2^10), SimulatedSNR15Series[,4], ty="l",
main=bquote(paste("SNR=15 Data, Threshold ", J[0],"=3")),
xlab="Time", ylab="Observation", col="gray", ylim=c(-.5, .6))
lines(seq(1, 2^10), SimulatedSmoothSeries[,4], col="red", lwd=2)
plot(seq(1, 2^10), SimulatedSNR25Series[,4], ty="l",
main=bquote(paste("SNR=25 Data, Threshold ", J[0],"=3")),
xlab="Time", ylab="Observation", col="gray", ylim=c(-.5, .6))
lines(seq(1, 2^10), SimulatedSmoothSeries[,4], col="red", lwd=2)
## ----sel4, echo=TRUE, eval=TRUE, out.width='1.1\\textwidth', out.height='17cm'----
smoothPlot <- smoothTimeSeries(SimulatedSNR15Series[ ,4], plot="all")
## ----sel5, echo=TRUE, eval=TRUE------------------------------------------
set.seed(1414)
SNR15Boot <- WiSEBoot(SimulatedSNR15Series[ ,4], R=10)
SNR15Boot$MSECriteria
## ----sel6, echo=TRUE, eval=TRUE, out.width='13cm', out.height='7cm'------
par(mfrow=c(1,2))
boxplot(SNR15Boot$BootIntercept,
main=expression(paste("R=10 Bootstrap Estimates of ", gamma[0])),
ylab=expression(hat(gamma)[0][b]))
abline(h=0, col="red")
boxplot(SNR15Boot$BootSlope,
main=expression(paste("R=10 Bootstrap Estimates of ", gamma[1])),
ylab=expression(hat(gamma)[1][b]))
abline(h=0, col="red")
## ----sel7, echo=TRUE, eval=TRUE, out.width='15cm', out.height='7cm'------
par(mfrow=c(1,2))
boxplot(SNR15Boot$BootWavelet[,3],
main=expression(paste("R=10 Boot. Est. of lvl=1, coef=1, ", gamma)),
ylab=expression(hat(gamma)[b]))
boxplot(SNR15Boot$BootWavelet[,4],
main=expression(paste("R=10 Boot. Est. of lvl=1, coef=2, ", gamma)),
ylab=expression(hat(gamma)[b]))
## ----hyp1, echo=TRUE, eval=TRUE------------------------------------------
data(CM20N20S60E)
CM20N20S60E[1:3,]
## ----hyp2, echo=TRUE, eval=TRUE, out.width='15cm', out.height='15cm'-----
par(mfrow=c(2,1))
plot.ts(CM20N20S60E[,1], main="AIRS", ylab="Obs. Climate")
plot.ts(CM20N20S60E[,10], main="MIROC5, run 5", ylab="Model Climate")
## ----hyp3, echo=TRUE, eval=TRUE------------------------------------------
pad60E <- padMatrix(CM20N20S60E)
dim(pad60E$xPad)
## ----hyp4, ehco=TRUE, eval=TRUE------------------------------------------
hypObj <- WiSEHypothesisTest(pad60E$xPad[,1], pad60E$xPad[,10], R=10, J0=5,
XParam=pad60E$linearParam[,1],
YParam=pad60E$linearParam[,10])
hypObj$AsymptoticPValue
hypObj$BootstrapPValue
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.