bootstrap.phenologyFit function

bootstrap.phenologyFit

This function bootstraps the residuals of a phenologyFit. It internally calls phenologyFitter on each bootstrap replicate.

bootstrap.phenologyFit(
  object,
  boot.R = 99,
  control = list(smooth = FALSE, verbose = FALSE, maxit = 1000, nb.stop.improvement =
    250),
  lower,
  upper,
  seed = 1766588
)

Arguments

• object: class phenologyFit, the object to bootstrap
• boot.R: integer. The number of bootstrap replicates
• control: control parameters to GenSA, see GenSA::GenSA
• lower: Vector with length of ‘par.guess’. Lower bounds for components.
• upper: Vector with length of ‘par.guess’. Upper bounds for components. If missing, upper in object is used.
• seed: integer seed for the random number generator used by GenSA. If missing, lower in object is used.

Returns

Invisibly returns a list with elements boot.R, object, seed, residuals, lower, upper, and res. The latter list res has boot.R elements, which are lists again. Each of these lists contains named elements par, value, bloomJDays, and pbloomJDays. par are the best fit parameters on the particular bootstrap replicate, value the corresponding RSS, bloomJDays the re-sampled data and pbloomJDays the predicted bloom JDays for this sample.

Details

bootstrap an object of S3 class phenologyFit

Author(s)

Carsten Urbach urbach@hiskp.uni-bonn.de