Bonjour,
Je suis actuellement en train d'essayer d'apprendre le langage R seul et je travaille en ce moment sur la fonction gpd du package VGAM.
Je n'arrive cependant pas à comprendre la syntaxe de celle ci.
Le but de mon petit exercice est qu'après avoir simulé une GPD avec rgpd je dois utiliser la fonction gpd pour calculer les estimateurs sigma et gamma par moments pondérés et maximum de vraisemblance.
Voila ce que j'ai écrit pour l'instant mais je ne vois pas comment avancer.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12 # packages-charger packages evir et evd library(evir) library(evd) library(VGAM) GPD<-rgpd(30) GPDc<-sort(GPD) seuil<-min(GPDc) gpd(??????)
Merci de vos conseils.
Partager