#Modélisation 3D du modèle Smith Kern d'apparition du dollar spot library(plotly) library(geoR) vtemp=seq(0,40) vhum=20:100 probab=matrix(nrow=length(vtemp),ncol=length(vhum)) for (i in 1:length(vtemp)){ for (j in 1:length(vhum)){ if(vtemp[i]<=10|vtemp[i]>=35){ probab[i,j]=0 } else{ logit=-11.4041+(0.0894*vhum[j])+(0.1932*vtemp[i]) probab[i,j]=(exp(logit)/(1+exp(logit)))*100 } } } #Affichage selon différents angles (theta et phi) par(mfrow=c(2,1)) persp(vtemp,vhum,probab,theta=30,phi=15,expand=0.5,col="lightblue",ticktype="detailed",xlab="température (%)", ylab="humid. relative (%)", zlab="proba. appar. dollar spot (%)") persp(vtemp,vhum,probab,theta=-40,phi=15,expand=0.5,col="lightblue",ticktype="detailed",xlab="température (%)", ylab="humid. relative (%)",zlab="proba. appar. dollar spot (%)")