Bonjour, Il doit en effet y avoir un pépin dans le calcul de SensitivityBiomassRelativeY6 dans la V4... Car que ce soit avec la V3 ou la V4, mes SensitivityBiomassY3 sont à première vue calculés correctement (même s'il est trop tôt pour dire s'ils sont justes ils sont au moins cohérents entre eux). SensitivityBiomassRelativeY6 étant le ratio entre SensitivityBiomassY3 (biomasse au dernier pas de temps)et la biomasse au début de la simulation, je ne vois pas d'où peuvent venir tous ces 1, à part d'une erreur dans le calcul de SensitivityBiomassRelativeY6 dans la V4. Loïc Jean Couteau a écrit :
-----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----- _______________________________________________ Isis-fish-devel mailing list Isis-fish-devel@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-devel