-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Bonjour, Le 10/11/2011 09:42, Loic GASCHE a écrit :
Bonjour,
j'ai lancé la même AS sous isis V3 puis sous isis V4, en utilisant la méthode de Morris.
Avec la V3 tout se passe bien et les sorties de la méthode de Morris (mu, mu.star et sigma) sont bien calculées. Par contre avec la V4, j'obtiens des NaN pour mu, mu.star et sigma. A première vue, cela vient du fait que la matrice a$ee vaut NaN et que c'est cette dernière qui est utilisée pour les calculs de la méthode de Morris.
Mais je ne sais pas pourquoi cette matrice vaut NaN avec la V4, donc si quelqu'un a une idée...
A priori, le problème vient du résultat Y6 qui est toujours égal à 1 quelque soit la simulation. Je n'est pas le problème avec Y3. Exemple dans R, d'abord avec Y3 puis Y6
a
Call: morris(model = NULL, factors = c("Sole7D_capturability", "Sole7D_naturalDeathRate_result"), r = 5, design = list(type = "oat", levels = c(5, 5), grid.jump = c(2, 2)), binf = c(0, 0), bsup = c(1, 1)) Model runs: 15 mu mu.star sigma Sole7D_capturability 1.306054 1.306054 0.4458962 Sole7D_naturalDeathRate_result 0.000000 0.000000 0.0000000
try(tell(a,y=SensitivityBiomassRelativeY6_Sole7D),silent=TRUE) a
Call: morris(model = NULL, factors = c("Sole7D_capturability", "Sole7D_naturalDeathRate_result"), r = 5, design = list(type = "oat", levels = c(5, 5), grid.jump = c(2, 2)), binf = c(0, 0), bsup = c(1, 1)) Model runs: 15 mu mu.star sigma Sole7D_capturability NaN NaN NA Sole7D_naturalDeathRate_result NaN NaN NA Dans la v3, tu as les mêmes résultats pour tes simus ? Est-ce que Y6 est toujours égal à 1 ? Jean -----BEGIN PGP SIGNATURE----- Version: GnuPG v2.0.18 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/ iQEcBAEBAgAGBQJOwMhnAAoJEFOQdnjKiPj334YH/2mH5NoZyqEJzGNIo/po+6tW si1lUTGhCKOwYeV8EdsZeDA/KhI+SbGJb8VLGqN84tGcCFSk77SY9r0CHebNoLE5 yaoptAKlPNe5yccLNitOOIaFvEU7RbeKrWa+xViDwNyHwDXtdthQT6njT8R8MWsF i8I3fWMVZXMi1snxOGQFqBVvNdIRstjYNa1cg5BvTOxk8fxtDBqFigUGKCorUS+W LMzyIr8G5EUAEzpe8Wl+/IWlNFtOb1OWuE4hYMLgdvk7/Iw3PQ8yq/w7yKsL5smA boYJ7mI7cvKWKmMcSk8Nf3QAUUxeoxDsbFm4QW4PDQ1pheZB/GPeEgMkNNDgqOQ= =JR0W -----END PGP SIGNATURE-----