shift_sinusoid: Shift sinusoid information.

In phenology: Tools to Manage a Parametric Function that Describes Phenology and More

Shift sinusoid information.

Description

This function is used to shift sinusoid parameters from ”, '1' or '2'.

Usage

shift_sinusoid(parameters = NULL, from = "", to = "1")

Arguments

parameters	set of parameters
from	The number of series to change
to	The number of series to change

Details

shift_sinusoid shift sinusoid information.

Value

Return a set of modified parameters

Author(s)

Marc Girondot

Examples

# Read a file with data
library("phenology")
data(Gratiot)
# Generate a formatted list nammed data_Gratiot 
data_Gratiot<-add_phenology(Gratiot, name="Complete", 
		reference=as.Date("2001-01-01"), format="%d/%m/%Y")
# Fix parameter FLat to 0
pfixed=c(Flat=0)
# Generate initial points for the optimisation
parg<-par_init(data_Gratiot, fixed.parameters=pfixed)
# Fit is done
## Not run: 
result_Gratiot_Flat<-fit_phenology(data=data_Gratiot, 
		fitted.parameters=parg, fixed.parameters=pfixed)

## End(Not run)
data(result_Gratiot_Flat)
parg<-extract_result(result_Gratiot_Flat)
# Add data for one sinusoid superimposed 
# [day d amplitude=(Alpha+Nd*Beta)^Tau with Nd being the number of counts for day d]
parg<-c(parg, Alpha=0.5, Beta=0.8, Delta=3, Phi=15)
# Tau is fixed to 1
pfixed=c(Flat=0, Tau=1)
# Run the optimisation
## Not run: 
result_Gratiot1<-fit_phenology(data=data_Gratiot, 
		fitted.parameters=parg, fixed.parameters=pfixed)
# Plot the phenology
output1<-plot(result_Gratiot1, moon=TRUE)
#' 
## End(Not run)
data(result_Gratiot1)
# Extract the fitted parameters
parg1<-extract_result(result_Gratiot1)
# Shift sunusoid information to the '1'
parg2<-shift_sinusoid(parameters=parg1, from="", to="1")
# Tau is fixed to 1
pfixed=c(Flat=0, Tau1=1, Tau=1)
# Add data for another sinusoid superimposed 
# [day d amplitude=(Alpha+Nd*Beta)^Tau with Nd being the number of counts for day d]
parg<-c(parg2, Alpha=0.5, Beta=0.8, Delta=3, Phi=10)
# Run the optimisation
## Not run: 
result_Gratiot2<-fit_phenology(data=data_Gratiot, 
		fitted.parameters=parg, fixed.parameters=pfixed)
# Plot the phenology
output2<-plot(result_Gratiot2, moon=TRUE)

## End(Not run)
data(result_Gratiot2)

phenology documentation built on Oct. 16, 2023, 9:06 a.m.