Merci Jean pour cette réponse, A l'origine ces questions émergent d'une autre : est-il possible de trouver dans les sorties du modèle les résultats de l'analyse de sensibilité correspondant aux "vraies" valeurs, c'est à dire ici celles appartenant à l'intervalle 0.8-1.2 pour la mortalité naturelle et pas à 0-1 ? Car pour le moment je n'ai trouvé de sorties que sous la forme de 3 colonnes: la standardisation qui est bien entre 0.8 et 1.2, la biomasse de fin, mais uniquement des valeurs entre 0 et 1 pour Ktemp. Ou alors faut il se débrouiller et convertir les valeurs entre 0 et 1 pour les faire correspondre à la bonne valeur dans l'intervalle souhaité ? Merci ! Loïc Jean Couteau a écrit :
Le 17/06/2011 10:38, Loic GASCHE a écrit :
Bonjour,
Bonjour,
Après vérification, mes interrogations demeurent en ce qui concerne l'utilisation de la méthode de Morris dans ISIS... Je souhaitais réaliser une analyse de sensibilité sur 2 facteurs (standardisation et Ktemp), le facteur de standardisation variant de manière continue entre 0.8 et 1.2 et Ktemp variant de 20% autour de sa valeur de base de 1. Pourtant les bornes appelées lors de l'appel de morris dans R ne correspondent pas à celles demandées (Ktemp a pour bornes 0 et 1)...
Dans le cas des équations continues (Ktemp) on a toujours 0 et 1 lors de l'appel de morris dans R, après on gère le retour pour le remettre dans les vraies bornes.
call morris(model = NULL, factors = c("BeamTrawl.standardisationFactor", "Sole7D.naturalDeathRate.Ktemp"), r = 4, design = list(type = "oat", levels = c(5, 5), grid.jump = c(2, 2)), binf = c(0.8, 0), bsup = c(1.2, 1))
Et les bornes indiquées par $isis.factors et $isis.factor.distribution sont également assez perturbantes...
$isis.factors nomFacteur Nominal Continu Binf Bsup 1 BeamTrawl.standardisationFactor 1 TRUE 0.8 1.2 2 Sole7D.naturalDeathRate.Ktemp 1 TRUE 0.8 1.2
Pas tant que ça en fait. Si tu pars sur une valeur nominale de 1 et une variation de 20%, tu as donc des "vraies" bornes de 0,8 et 1,2
$isis.factor.distribution NomFacteur NomDistribution ParametreDistribution 1 BeamTrawl.standardisationFactor qunif [0.8;1.2] 2 Sole7D.naturalDeathRate.Ktemp qunif [0.0;1.0]
Ici c'est les bornes de l'appel à la méthode de morris.
Dans la nouvelle version on a unifié tout ça et ça devrait être beaucoup moins perturbant (en tout cas j'espère).
Jean _______________________________________________ Isis-fish-users mailing list Isis-fish-users@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-users