Print a sea
object.
# S3 method for sea
print(x, ...)
A sea
object.
Additional arguments that are tossed.
sea()
creates a sea
object.
plot_sealags()
basic plot of sea
object lags.
if (FALSE) {
# Read in the Cook and Krusic (2004; The North American Drought Atlas)
# reconstruction of Palmer Drought Severity Index (PDSI) for the Jemez
# Mountains area (gridpoint 133).
target_url <- paste0(
"http://iridl.ldeo.columbia.edu",
"/SOURCES/.LDEO/.TRL/.NADA2004",
"/pdsiatlashtml/pdsiwebdata/1050w_350n_133.txt"
)
pdsi <- read.table(target_url, header = TRUE, row.names = 1)
pdsi <- subset(pdsi, select = "RECON")
# Run SEA on Peggy Mesa (pgm) data
data(pgm)
pgm_comp <- composite(pgm)
pgm_sea <- sea(pdsi, pgm_comp)
# See basic results:
print(pgm_sea)
# Basic plot:
plot(pgm_sea)
}