phenology (version 7.2)

shift_sinusoid: 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

Value

Return a set of modified parameters

Details

shift_sinusoid shift sinusoid information.

Examples

# NOT RUN {
# Read a file with data
library("phenology")
# }
# NOT RUN {
Gratiot<-read.delim("http://max2.ese.u-psud.fr/epc/conservation/BI/Complete.txt", header=FALSE)
# }
# NOT RUN {
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)
# }
# 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)
#' 
# }
# 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)
# }
# NOT RUN {
data(result_Gratiot2)
# }