j'ai fait une analyse de sensibilité de Morris avec le package "sensitivity" de R. J'ai réaliser des scatter plots afin d'identifier les paramètres (7 au total) qui impact le plus les sorties A,B et C du modèle.
Seulement, j'aimerais rajouter sur chaque graphique trois droites de pente sigma/mu*=0.1,0.5 et 1 (voici un exemple en page 6 de cet article: https://www.sciencedirect.com/science/a ... 8812004768).
J'ai calculé la valeur de mu* et sigma de chaque paramètre mais je ne vois pas comment obtenir et tracer ces droites sur les graphiques.
J'ai mis une partie de mon script en bas, je ne sais pas si cela va aider.
En vous remerciant
Code : Tout sélectionner
etude.morris <- morris( model = NULL, factors = facteurs$name, r = 60,
design = list(type = "oat", levels = 10, grid.jump = 10/2),
scale=T, binf=facteurs$binf, bsup=facteurs$bsup)
Sorties.morris <- model.func(etude.morris$X) #AS
Sorties.morris <- as.data.frame(Sorties.morris)
Nsorties <- c("A", "B", "C") # Output names
colnames(Sorties.morris) <- Nsorties
# Value of mu, mu* and sigma for each parameters:
mu <- apply(etude.morris$ee, 2, mean)
mu.star <- apply(etude.morris$ee, 2, function(etude.morris) mean(abs(etude.morris)))
sigma <- apply(etude.morris$ee, 2, sd)
# Morris scatter plot mu* vs sigma :
windows()
op = par(mfrow=c(2,2))
for(i in 1:3){
etude.morris = tell(x=etude.morris, y=Sorties.morris[,i])
plot(etude.morris, col=rainbow(7))
title(Nsorties[i])
# abline(v=0.01,col="red")
}