utilisaton de persp

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Olivier Orcel
Messages : 29
Enregistré le : 28 Nov 2014, 09:51

utilisaton de persp

Messagepar Olivier Orcel » 06 Avr 2016, 12:45

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 ?

Par avance, je vous remercie.

PS : je vous transmets copie partiel du prgramme.


Copule<-function(){
F10 = 1/7060
H10 = 2.4806*(-log(F10) - 2.8117)^{0.32}
S10 = - 0.1095*log(F10) - 0.96
#=====================================================
#================LECTURE DES DONNES===================
#=====================================================
tab<-read.table("TRII.txt",h=T)
tab1<-as.matrix(tab)
#=====================================================
#
for (S in seq(1,10,by=1)){
if (S > 0.3){fs = exp(-((S+0.96)/0.1095))}
else{fs=1exp(0.0711/(S0.029))}
x = 0.56
fHT = abs(1 - fs^{-x} + (F10)^{-x})^{-1/x}
#=====================================================
#
for (H in seq(1,10,by=1)){
if (H < 3){fH = exp(-((H-0.55)/0.53))}
else {fH = exp(-(((H/2.4806)^{1/0.32} + 2.8117)))}
x = 2.37
fT = abs(1 - fH^{-x} + (fHT)^{-x})^{-1/x}
#=====================================================
#===================DETERMINATION DE T================
#=====================================================
if (fT > 6){T = exp((fT - 5)/2)}
else {T = -0.11*log(fT)}
#=====================================================
#=========REPRESENTATION GRAPHIQUE ISOVALEURS 3D======
#=====================================================
#
f<-function(S,H){ )
# S<-seq(1,10,by=1)
# H<-S
# z<-outer(S,H,f)
# persp(S,H,z,theta=30,phi=30,expand=0.5)
}
}
plot(S,H,T)
}
Copule()
#====================================================
Bonjour,

Pouvez vous me dire comment utiliser les copules de Marshall et de Pareto ?

Ca devrait etre ( comme pour les autres copules ) :

pareto.cop<-paretoCopula(400)
fit.mpl<-fitCopula(pareto.cop,...)
fit.mpl

Je vous remercie.

O ORCEL

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité