plotSum {phenopix} | R Documentation |
A plotting facility to show the combination of uncertainty
Description
The combineUncertainty
uses greenProcess
to fit all available double logistic equations in the phenopix package and extracts thresholds with all available methods. Then uncertainties can be combined and returned by using summarizePhases
and plotted with plotSum
.
See greenProcess
.
Usage
plotSum(ts, sum, which, v=NULL, quantile=TRUE, ...)
Arguments
ts |
A ts or zoo object with gcc data. index(ts) must be numeric days of year (doys) or a POSIXct vector |
sum |
An object in output from |
which |
One between trs, derivatives, klosterman, gu |
v |
An optional vector of vertical coordinates (in y-axis unit for plot annotation of phase names) |
quantile |
If TRUE, the plotted uncertainty envelope is based on the quantiles, and not min-max, otherwise min-max envelope is plotted |
... |
For the plotting function, a number of parameters from generic |
Details
This function is the last step of a processing chain. It uses greenProcess
to fit all available double logistic equations in the phenopix package and extracts thresholds with all available methods. Then uncertainties can be combined and returned by using summarizePhases
and plotted with plotSum
.
See greenProcess
, summarizePhases
, plotSum
. This function uses a model approach to combine all uncertainties from all available phenopix fittings, as to get an ensemble of phases with different methods, without necessarily choosing any of them.
Value
A named list with dataframes for each phenophase method with all replication for each of the included fitting methods. These data can then be combined with the companion functions summarizePhases
and plotSum
. See examples for details.
Author(s)
Gianluca Filippa <gian.filippa@gmail.com>
Examples
## Not run:
require(zoo)
data(bartlett2009.filtered)
combined.fit <- combineUncertainty(na.approx(filtered.tmp$max.filtered), nrep=100)
# 100 replications for each fitting
names(combined.fit) # a dataframe for each phenoMethod + a list with all fittings
fit.summary <- summarizePhases(combined.fit, across.methods=TRUE)
## again a list with one element for each fitting method + two additional items
## if across.methods is TRUE, which combines gu + klosterman phenophase methods
## in a single method, and the same happens for trs and derivatives
plotSum(bartlett2009.filtered, fit.summary, which='klosterman')
## a plot with original timeseries + phenophases and their uncertainty
## End(Not run)