Description Usage Arguments Details Value Author(s) References See Also Examples
Computes paleo-fire frequency for a set of fire events (or frequency from other events types, see examples) using a gaussian kernel density estimation procedure based on a defined bandwidth (see Mudelsee 2004 for details). Pseudo-replicated values are used to correct for edge bias, equivalent to "minimum slope" correction in Mann (2004).
1 2 3 4 5 6 7 8 9 10 11 12 13 |
fevent |
Numeric vector, set of dates |
up |
Numeric, upper age for fire frequency calculus |
lo |
Numeric, lower age for fire frequency calculus |
interval |
Numeric, interval between two points for fire frequency calculus (default 10 years) |
bandwidth |
Numeric, bandwidth in years, or character for automatic
bandwidth calculation (e.g. "bw.ucv" for unbiased cross validation) see
|
boot |
Character, "full" or "partial" see @details |
bootper |
Numeric, percentage of fire events randomly added or removed in the "partial" replication procedure (default 0.1) |
nbboot |
Numeric, number of bootstrap replicates |
alpha |
Numeric, confidence interval (default 0.01) |
pseudo |
Logical, apply (TRUE) or not (FALSE) Mann (2004) correction (default=FALSE) |
pseudo_per |
percentage of actual data used in reflection in the Mann (2004) correction |
By using boot="partial" option (beta!) fire dates are randomly removed or added within a defined percentage (by default between 1 and 10% of total number of events) in order to make new series that are then used to calculate ensemble members fire frequencies. This procedure differs slightly from the full bootstrapp where fire dates are randomly picked with replacement. Theoretically classic bootstrap could result in a sample where a single fire event date is replicated n times which makes no sense for fires. By randomly removing or adding fire dates the confidence intervals are narrower and likely better reflect the long term fire regime variablility.
ff data.frame, with fire frequency, bandwidth and CIs
O. Blarquez
Mann, M. E. (2004). On smoothing potentially non-stationary
climate time series. Geophysical Research Letters, 31(7).
Mudelsee, M., B<c3><b6>rngen, M., Tetzlaff, G., & Gr<c3><bc>newald, U. (2004). Extreme floods
in central Europe over the past 500 years: Role of cyclone pathway <e2><80><9c>Zugstrasse Vb<e2><80><9d>.
Journal of Geophysical Research: Atmospheres (1984<e2><80><93>2012), 109(D23).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ## Not run:
set.seed(123)
fevent=c(round(abs(rnorm(20,mean=7,sd=5))*1000),round(abs(rnorm(10,mean=8,sd=1))*1000))
ff=kdffreq(fevent,bandwidth = 1000, nbboot=10)
# Estimate the frequency of armed conflicts from 1946 to 2014
# Data from the The Uppsala Conflict Data Program (UCDP) available at: https://www.prio.org
dat=read.csv('http://ucdp.uu.se/downloads/ucdpprio/ucdp-prio-acd-4-2016.csv')
res=kdffreq(dat$Year,bandwidth = "bw.ucv", nbboot=1000, up = 1946, lo = 2014, interval=1, pseudo=T)
plot(res, ylab="# armed conflict/year")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.