phenology (version 3.40)

likelihood_phenology: Estimate the likelihood of timeseries based on a set of parameters.

Description

This function is used to estimate the likelihood based on a set of parameters.

Usage

likelihood_phenology(data = NULL, parametersfit = NULL,
    parametersfixed = NULL, zero_counts = NULL,
    method_incertitude = NULL, result = NULL)

Arguments

data

Dataset generated with add_format

parametersfixed

Set of fixed parameters

parametersfit

Set of parameters to be fitted

method_incertitude

2 [default] is the correct one from a statistical point of view; 0 is an aproximate method more rapid; 1 is an alternative more rapid but biased.

zero_counts

example c(TRUE, TRUE, FALSE) indicates whether the zeros have been recorder for each of these timeseries. Defaut is TRUE for all.

result

An object obtained after fit_phenology()

Value

• The likelihood of the data with the parameters

Details

likelihood_phenology estimate likelihood for a set of parameters.

Examples

# Read a file with data
# Gratiot<-read.delim("http://max2.ese.u-psud.fr/epc/conservation/BI/Complete.txt", , header=FALSE)
data(Gratiot)
# Generate a formatted list nammed data_Gratiot
data_Gratiot<-add_format(origin=NULL, add=Gratiot, name="Complete", reference=as.Date("2001-01-01"), format="%d/%m/%Y")
# Generate initial points for the optimisation
parg<-par_init(data_Gratiot, parametersfixed=NULL)
# Estimate likelihood with this initial set of parameters
likelihood_phenology(data=data_Gratiot, parametersfit=parg, parametersfixed=NULL)