Bonjour Je cherche à obtenir le maximum de vraisemblance avec la copule de Clayton : je pensais que c'était LogLik Clayton.cop<-claytonCopula(5887,2) fit.mpl<-fitCopula(Clayton.cop,FHS,method="mpl") fit.mpl print(coef(fit.mpl)) (ce qui donne bien le paramètre de la copule) Merci par avance...
Bonjour je souhaiterais superposer pluisieurs courbes agec contour. Je fais : contour(NH,NS,Tretour1,xlab="Hauteur de houle en m",ylab="Surcote en m",levels=c(10,100,1000), labelc=1.,col="blue",drawlabels=TRUE) contour(NH,NS,Tretour2,xlab="Hauteur de houle en m&quo...
j'avais effectivement trouvé (il y a une petite heure) qu'il convenait d'écrire U = uniroot(... puis de faire print(U), parfois des problèmes de syntaxe font perdre des heures.
Bonjour, Je souhaiterais obtenir toutes valeurs de x qui annulent x+H, avec H variant de 1 à 10. Pour cela, je fais : for (H in 1:10){ FG<-function(x){x+H} uniroot(FG,lower=1,upper=10) } Je ne récupère que la valeur 10. Merci de m'aider pour ce problème basique (peut être trop fatigué !). Olivier
Je suis bloqué depuis qq temps sur l'utilisation de wireframe. Je dois superposer deux graphiques en 3d en R. Si vous pouviez me donner un exemple avec wireframe, je vous serais très reconnaissant.
Pour plus de renseignements : Copule<-function(){ S<-seq(1,10,by=1) H<-seq(1,10,by=1) PERIODE<-matrix(nrow=10,ncol=10,byrow=TRUE) PERIODE1<-matrix(nrow=10,ncol=10,byrow=TRUE) # PERIODE[S,H] = fT PERIODE1[S,H] = fT1 #==fT et fT1 sont calculées en fonctions de S et H et de vraiables # z<-PERIODE[S,H] ...
Je dois visualiser deux courbes avec perps. D'après ce que j'ai pu voir sur intenet, il faut faire persp(x,y,z,theta=30,phi=40,..) perps(x,y,z,theta=40,phi=30,..)
Bref, pouvez vous m'aider, je vous en serais très reconnaissant.
Bonjour, Je dois visualiser en 3D trois variables (S,H,T). Je ne peux pas exprimer T en fonction des variables S et H, donc il me semble difficile de recourir à la fonction persp. En effet, l'expression de T s'exprime par : T = exp((fT-5)/2) Que me suggérez vous pour représenter graphiquement en 3d ...
Bonjour, Afin d'ajuster une loi statistique à un échantillon, j'utilise la méthode POT. Avec le package Renext, je fais Renouv. J'ai besoin de déterminer le threshold. Pouvez vous m'aider ? Je vous remercie beaucoup. O. ORCEL PS : j'ai trouvé un exemple dans lequel il était possible d'utiliser clust...