bootstrap.phenologyFit {chillR}

bootstrap.phenologyFit

Description

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

Usage

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.

Details

bootstrap an object of S3 class 'phenologyFit'

Value

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.

Author(s)

Carsten Urbach <[e-mail]>