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