Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
out.width = "47%",
fig.show = "hold",
fig.ext = "svg",
dev = "svg"
)
# Load problematic correlogram for example
pcorr <- structure(c(0.500528302565943, -0.754836705515775, -1.19039679450543,
-0.178590255566502, -0.60652246495491, 1.82273156499271, -0.172010338128128,
-2.0163653383661, -0.3745369917274, 0.456123617810818, 1.67415100748476,
0.937900901537193, -1.29338484030548, -1.28069618875814, -1.82633767163309,
-0.614386402460464, 0.584995100586833, -0.504037627819764, 0.0226448302823769,
-1.20080447279195, -0.0617993438737834, 0.528551859910452, 1.96448337848478,
0.072350278901269, 0.225152117571247, -1.25610158595963, -0.845732815433606,
0.471680065854876, 2.4943925938456, 1.80772326821158, 0.204809764053051,
-0.150961614287815, -0.0197795008143637, -1.05994768012253, -1.29810950765603,
-1.17853411156364, -0.643232293492797, -2.20411574357046, 0.216829082430902,
-0.941814658830531, -0.781300200015838, 0.273580690612867, 1.09702483374736,
-1.36990396610849, -2.07605582194368, 0.451160815513836, 0.86688665896418,
0.51840092005946, 1.28921143353651, 0.875645430906388, 0.591064637882466,
-0.55966314736262, -0.367191668064688, -2.76744332433181, -1.2443630247938,
-0.154933017518784, 0.0667021828181572, 0.197007379842123, 0.490808179023661,
2.2336743342413, 1.89448119625069, 1.47138565779839, -0.666213196115124,
0.843870629377204, -0.203203819755488, -2.20410902950622, 0.827755498838514,
1.78219490528679, -0.111762758688384, -2.22360645019641, -1.0056276510069,
-0.329364603192913, 1.33056799025663, 0.149121700442515, 0.267676241642407,
-0.679653050876458, 0.0901160031947412, -0.563650739275465, 0.666301686016796,
1.70298732225847, 0.673722209645992, 1.32892165992321, -1.29016755347508,
-2.11926333043274, -0.434170672177059, -1.11952602901843, -0.282918307090735,
-0.296513315892018, -0.215455257897354, -0.229591553459229, 1.90335315883752,
-0.640041683494477, 0.555037751515986, 0.498140844212503, -0.295915459399637,
0.587207819554741, -0.52622399410093, 1.32416109800959, -0.825699074434715,
-2.51977385384327, -0.681811843738689, 1.0515603999965, 1.75547334518565,
1.46352961234024, 0.0972088740341473, -1.41469490750948, -2.39336234402329,
0.596660462796983, 2.29853427605131, 1.65898229920668, -0.0276102807086275,
-1.64374943302099, 0.602395878758738, 0.229442193055737, 2.1687970121186,
0.643910340438762, -0.290718077669911, -1.02619172899549, -0.433490300518998,
-0.978870909554381, -0.0098964640725574, -1.96478096657371, -0.0969197029820769,
-2.08344229475593, -0.293555313118738, -1.29271002623495, 0.768226868279432,
0.714533099091565, -0.739381839162954, 0.486811684636613, -0.0462968073811469,
0.119532030148799, 0.699661964254544, -0.121299551979732, -0.952343329737353,
0.47622151198097, -0.235031006098588, -0.274831514573784, -0.0760339102360359,
1.25475007283353, 1.46495001193988, 0.167876131248247, 0.191014261173205,
1.20954335642528, 1.84131238261517, 1.11873900227391, -0.948077984555249,
-1.06745345491371, -1.4245470888048, -0.538478752695991, -0.277037416300269,
-0.887735440660801, -0.960970972045782, -1.17008385676073, 0.291022149761318,
1.57923122715123, 2.55155997330415, 0.225931504246328, 0.00220891571648253,
-0.416842778537851, -0.607033513606461, -0.40792204210633, 0.0685179557733119,
1.14003218398751, 1.82522025910671, 1.35481822913805, 1.52147329773245,
0.328519767736851, -0.0225431382950235, -0.0632294307946625,
-0.34692001495444, 0.609522181403036, 1.18523376254118, -0.2922169564903,
0.486710723126883, 0.0171559573700382, -0.841704196640039, -0.170564928787282,
1.09088790395037, 1.01711383301683, -0.57576763235922, -0.813788055144079,
0.845138966533092, -0.506617097735509, -0.180391805588545, -0.0731710033636772,
0.650969458632428, -0.544567330307273, -0.697447901212486, -0.843571490826714,
-1.04624243603828, 0.862168240020156, 0.393550845622028, -1.48053698592407,
-1.750632331433, 1.47944447592151, 1.1384887083871, 1.19239984208863,
1.05781365301425, 0.231828681437716), .Tsp = c(1, 200, 1), class = "ts")
## ----setup--------------------------------------------------------------------
library(blocklength)
set.seed(32)
## ----Inspect Parameters-------------------------------------------------------
# Generate an AR(1) simulation
sim <- stats::arima.sim(list(order = c(1, 0, 0), ar = 0.5),
n = 500, innov = rnorm(500))
# Run the PWSD method and suppress output of the correlogram
b <- pwsd(sim, correlogram = FALSE)
# Inspect interim parameters
b$parameters
## ----Inpect b_max-------------------------------------------------------------
# Print block-length estimates
b$BlockLength
# Test if optimal block-length is censored by b_max
b$parameters[, "b_max"] == b$BlockLength
# Print the upper-bound
b$parameters[, "b_max"]
## ----Correlogram--------------------------------------------------------------
# Print rho_k_critical
b$parameters[, "rho_k_critical"]
# Plot correlogram from previous selection
plot(b, main = "c = qnorm(0.975)", ylim = c(-0.2, 1))
# Expand confidence level and plot again
b1 <- pwsd(sim, c = 4, correlogram = FALSE)
plot(b1, main = "c = 4", ylim = c(-0.2, 1))
## ----Convex MSE---------------------------------------------------------------
# Select a block-length using HHJ method
b2 <- hhj(sim, sub_sample = 10, k = "bias/variance", plots = FALSE)
plot(b2)
## ----Concave MSE--------------------------------------------------------------
# Generate a AR(2) simulation
sim2 <- stats::arima.sim(list(order = c(2, 0, 0), ar = c(0.5, 0.3)),
n = 500, innov = rnorm(500))
# Select a block-length using HHJ method
b3 <- hhj(sim2, plots = FALSE)
plot(b3)
## ----Problematic Correlogram, fig.align = "center", out.width = "80%"---------
# Select a block-length using PWSD method
b4 <- pwsd(pcorr, K_N = 5, correlogram = FALSE)
plot(b4, main = "Problematic Correlogram")
# Output m_hat to console
b4$parameters[, "m_hat"]
## ----Problematic Correlogram 2, fig.align = "center", out.width = "80%"-------
# Select block-length using PWSD method
b5 <- pwsd(pcorr, K_N = 6, correlogram = FALSE)
plot(b5, main = "Problematic Correlogram 2")
# Output m_hat to console
b5$parameters[, "m_hat"]
## ----Problematic Correlogram 3, fig.align = "center", out.width = "80%"-------
# Select a block-length using PWSD method
b6 <- pwsd(pcorr, c = qnorm(0.995), correlogram = FALSE)
plot(b6, main = "Problematic Correlogram 3")
# Output m_hat to console
b6$parameters[, "m_hat"]
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.