phase_plot | R Documentation |
The function creates a plot showing the three distinct phases of contraction, expansion and stem-radius increment (Deslauriers et al. 2011) for dendrometer series from a data.frame
as produced by phase_def
.
phase_plot(dm.gpf, dm.phase, sensor = NULL, period = NULL, colPhases = NULL, ...)
dm.gpf |
a |
dm.phase |
a |
sensor |
a |
period |
a |
colPhases |
a |
... |
additional graphical parameters (see |
The function plots phases of contraction, expansion and stem-radius increment along (one or more) dendrometer series. If more series are plotted (default), colors for the different lines can be defined using the col
argument for graphical devices (see par
). Note: if there are not enough custom colors, the function will repeat the last one used. If no colors are defined, the current palette
will be used.
The time axis will be automatically labeled depending upon the length of the dendrometer series. If period
is specified using a numeric
, DOY values are displayed on the x-axis. In case a character
of two time stamps is provided, axis labeling will be as follows: if series are longer than 120 days, years and months will be shown. If the length is between 30 and 120 days, months and days, and below 30 days, months, days and hours are displayed.
Plot showing stem-cyclic phases on dendrometer series.
Marko Smiljanic
Deslauriers, A., Rossi, S., Turcotte, A., Morin, H. and Krause, C. (2011) A three-step procedure in SAS to analyze the time series from automatic dendrometers. Dendrochronologia 29: 151-161.
data(dmCD)
dm.phase <- phase_def(dmCD)
phase_plot(dmCD, dm.phase)
# zoom in on the dendrometer series
phase_plot(dmCD, dm.phase, period = c(133, 142))
# customization options
phase_plot(dmCD, dm.phase, period = c("2008-05-12", "2008-05-22"),
colPhases = c("green", "cyan", "orange"),
pch = 4, main = "Dendrometer", ylab = "Values")
# specific sensors may be selected as follows:
data(dmED)
dm.gpf <- fill_gaps(dmED)
dm.phase <- phase_def(dm.gpf)
phase_plot(dm.gpf, dm.phase, sensor = 1)
phase_plot(dm.gpf, dm.phase, sensor = c(2,1))
phase_plot(dm.gpf, dm.phase, sensor = "Beech03")
phase_plot(dm.gpf, dm.phase, sensor = c("Beech03", "Beech04"))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.