Description Usage Arguments Details Value Author(s) Examples

View source: R/combineUncertainty.R

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`

.

1 | ```
combineUncertainty(ts, which='all', nrep=50, ncores='all')
``` |

`ts` |
A ts or zoo object with gcc data. index(ts) must be numeric days of year (doys) or a POSIXct vector |

`which` |
It can be 'all' (default) and all 4 double logistic fits will be calculated (beck, elmore, elosterman, gu), or a vector of subsets of the four fits |

`nrep` |
Number of relications to estimate uncertainty for each single fitting, defaults to 50. |

`ncores` |
Number of processors to be used in parallel computation, defaults to 'all' which will accidentally slow down any other activity on your computer. Otherwise set the number of processors you want to use in parallelization. |

This function 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 modellistic 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.

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.

Gianluca Filippa <gian.filippa@gmail.com>

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ```
## Not run:
require(zoo)
data(bartlett2009.filtered)
combined.fit <- combineUncertainty(na.approx(bartlett2009.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)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.