View source: R/summary.phenologymap.R
summary.phenologymap | R Documentation |
summary.phenologymap print information on a phenologymap object
## S3 method for class 'phenologymap'
summary(object, ...)
object |
A map generated with map_phenology. |
... |
Not used |
Return None
Marc Girondot
Other Phenology model:
AutoFitPhenology()
,
BE_to_LBLE()
,
Gratiot
,
LBLE_to_BE()
,
LBLE_to_L()
,
L_to_LBLE()
,
MarineTurtles_2002
,
MinBMinE_to_Min()
,
adapt_parameters()
,
add_SE()
,
add_phenology()
,
extract_result()
,
fit_phenology()
,
likelihood_phenology()
,
logLik.phenology()
,
map_Gratiot
,
map_phenology()
,
par_init()
,
phenology()
,
phenology2fitRMU()
,
phenology_MHmcmc()
,
phenology_MHmcmc_p()
,
plot.phenology()
,
plot.phenologymap()
,
plot_delta()
,
plot_phi()
,
print.phenology()
,
print.phenologymap()
,
print.phenologyout()
,
remove_site()
,
result_Gratiot
,
result_Gratiot1
,
result_Gratiot2
,
result_Gratiot_Flat
,
result_Gratiot_mcmc
,
summary.phenology()
,
summary.phenologyout()
library("phenology")
# Read a file with data
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")
# Generate initial points for the optimisation
parg<-par_init(data_Gratiot, fixed.parameters=NULL)
# Run the optimisation
## Not run:
result_Gratiot<-fit_phenology(data=data_Gratiot,
fitted.parameters=parg, fixed.parameters=NULL)
## End(Not run)
data(result_Gratiot)
# Extract the fitted parameters
parg1<-extract_result(result_Gratiot)
# Add constant Alpha and Tau values
# [day d amplitude=(Alpha+Nd*Beta)^Tau with Nd being the number of counts for day d]
pfixed<-c(parg1, Alpha=0, Tau=1)
pfixed<-pfixed[-which(names(pfixed)=="Theta")]
# The only fitted parameter will be Beta
parg2<-c(Beta=0.5, parg1["Theta"])
# Generate a likelihood map
# [default Phi=seq(from=0.1, to=20, length.out=100) but it is very long]
# Take care, it takes 20 hours ! The data map_Gratiot has the result
## Not run:
map_Gratiot<-map_phenology(data=data_Gratiot,
Phi=seq(from=0.1, to=20, length.out=100),
fitted.parameters=parg2, fixed.parameters=pfixed)
## End(Not run)
data(map_Gratiot)
# Print the information on a map
summary(map_Gratiot)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.