| daily.data | R Documentation |
Computes daily statistics from high-frequency dendrometer time series. For each day, it extracts the minimum and maximum values, their times of occurrence, daily mean and median, daily amplitude, the signed lag between the time of maximum and minimum, the day-to-day change in the daily maximum, and a daily status indicating whether the day is growing, shrinking, or stable relative to the previous day.
daily.data(df, TreeNum)
df |
A data frame with the first column containing date-time stamps
convertible to |
TreeNum |
Integer. The index of the tree (column) to analyze.
|
The function requires a data frame with a time column in the first column
and one or more dendrometer series in the following columns. The user
selects the series using TreeNum.
The returned object has class "daily_output", so it can be plotted
directly with plot().
The column Max_diff is computed as:
Max\_diff_t = Max_t - Max_{t-1}
The column Day_status is derived from Max_diff:
"growing" if Max_diff > 0
"shrinking" if Max_diff < 0
"stable" if Max_diff = 0
The first day has NA for Max_diff and Day_status.
A tibble of class "daily_output" containing:
Calendar date.
Daily minimum value.
Time of day of minimum value.
Daily maximum value.
Time of day of maximum value.
Daily mean value.
Daily median value.
Daily amplitude = Max - Min.
Time of minimum expressed in decimal hours.
Time of maximum expressed in decimal hours.
Signed difference in hours: Time_max_h - Time_min_h.
"*" if Time_max > Time_min, otherwise "".
Difference between today's maximum and the previous day's maximum.
"growing", "shrinking", "stable", or NA.
The object returned by daily.data() can be plotted using
plot() because it is assigned class "daily_output".
King G, Fonti P, Nievergelt D, Büntgen U, Frank D (2013) Climatic drivers of hourly to yearly tree radius variations along a 6°C natural warming gradient. Agricultural and Forest Meteorology 168:36–46. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/j.agrformet.2012.08.002")}
data(nepa17)
daily_stats <- daily.data(df = nepa17[1:1000, ], TreeNum = 1)
head(daily_stats, 10)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.