diff --git a/extreme_fit/distribution/gev/main_evgan.R b/extreme_fit/distribution/gev/main_evgan.R index 21b400794d5ceb8343e447a7c4d38d210502334b..1ba5828b14962ad9123d3dc3c79b57b003e36b7a 100644 --- a/extreme_fit/distribution/gev/main_evgan.R +++ b/extreme_fit/distribution/gev/main_evgan.R @@ -3,10 +3,8 @@ # Created by: erwan # Created on: 30/03/2021 -library(mgcv) -# library(evgam) +library(evgam) library(SpatialExtremes) -library(ismev) # Sample from a GEV set.seed(42) N <- 50 @@ -22,6 +20,8 @@ colnames(coord) = c("T") print(coord) coord = data.frame(coord, stringsAsFactors = TRUE) # res = fevd_fixed(x_gev, data=coord, method='MLE', verbose=TRUE, use.phi=FALSE) -res = fevd_fixed(x_gev, data=coord, location.fun= ~T, scale.fun= ~T, method='MLE', type="GEV", verbose=FALSE, use.phi=FALSE) +# fmla = list(miaxi) +evgam() +# res = evgam(x_gev, data=coord, location.fun= ~T, scale.fun= ~T, method='MLE', type="GEV", verbose=FALSE, use.phi=FALSE) # res = fevd_fixed(x_gev, data=coord, shape.fun= ~T, method='MLE', type="GEV", verbose=FALSE, use.phi=FALSE) print(res)