j'ai refais tes simus en mettant le meme nombre de modalités à tous les facteurs 5. Normalement on ne devrait pas avoir à saisir ce parametre pour chaque facteur puisque Morris demande à ce que ce nombre soit le meme pour tous les facteurs. Ca c'est un bug que les Lutins doivent corriger. Le fait que tu aies eu des NA pourrait etre expliqué par le 0 qui trainait car la seule division faite pour calculer les IA de morris est la division par 1/(p-1) avec p le nombre de modalités. Pourtant quand je corrige 0 par 5, avec une duree de simu de un an, j'obtiens des results dans les csv tous égaux à 1 (Y6 et Y4) alors que si je regarde les series temporelles il y a bien une difference entre le debut et la fin de la simu (à t=12 Y=7.56E9 et à t=0 Y=8.11E9) Est-ce un probleme d'arrondis? Eric comment peut-on avoir 7.56E9/8.11E9=1? Ensuite j'ai des NA pour les IA de Y6 mais pas de Y4. du coup je ne comprends pas le pourquoi des NA. des 0 c'est possible, des NAs non.... stephanie Bastien PREUSS a écrit :

J'ai 8 paramètres et j'ai un r = 4 ce qui fait 36 simus. Je viens d'essayer avec r = 20 ce qui fait 180 simus, mais j'ai toujours des valeurs égales à NA... Il ne me semble pas que le calcul de ces indices soit contraint par un nombre de simu, par exemple sur l'exemple golf de Gascogne en juin on avait testé r=4 pour 7 paramètres (soit 32 simu) et ça marchait sans soucis.

Du coup je ne vois pas trop d'où cela pourrait venir? Quelqu'un aurait-il une idée?

Merci.

Bastien

-----Message d'origine----- De : isis-fish-users-bounces@list.isis-fish.org [mailto:isis-fish-users-bounces@list.isis-fish.org] De la part de Jean Couteau Envoyé : Wednesday, December 08, 2010 2:47 PM À : isis-fish-users@list.isis-fish.org Objet : Re: [Isis-fish-users] NA dans les indices de sensibilité de Morris

Le 08/12/2010 14:28, Bastien PREUSS a écrit :

...

Je lance un AS avec Morris, mon plan comprend 35 simu, ça tourne parfaitement, toutes les simu se terminent, mais quand je vais voir dans

les

...

résultats, mes indices de sensibilité ont la valeur NA partout. J’ai

regardé

...

dans les exports .csv et dans le RData, c’est pareil.

Y aurait-il une raison connue à cela ?

Oui, la fonction Morris n'arrive pas à calculer les indices de sensibilité avec les résultats que tu lui donnes (d'après le debug, tu as bien des résultats pour toutes tes simus, mais à la sortie de la fonction morris, tu as des NA pour non calculé en gros). A première vue je dirais qu'il te faut plus de simulations pour calculer les indices de sensibilité. Stéphanie tu confirmes ?

Un moyen de tester avant de lancer les simus ISIS (qui sont longues quand même :D) : Dans R, tu crées une fonction qui a n paramètres, n étant le nombre de paramètres de ton AS (exemple : X1*X2*X...*Xn) et tu lances l'AS que tu veux sur cette fonction (qui simule ton modèle ISIS). Si tu as des résultats à ton AS, c'est bon tu peux lancer la même AS dans ISIS. (Note, j'ai jamais tester cette méthode, mais ça devrais le faire normalement).

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

_______________________________________________ Isis-fish-users mailing list Isis-fish-users@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-users

-- ...................................................................... Stephanie MAHEVAS (Stephanie.Mahevas@ifremer.fr) IFREMER/EMH (Ecologie et Modèles pour l'Halieutique) Tel: (33) 2 40 37 41 81 Fax: (33) 2 40 37 40 75 o \ o / _ o __| \ / |__ o _ \ o / o /|\ | /\ ___\o \o | o/ o/__ /\ | /|\ / \ / \ | \ /) | ( \ /o\ / ) | (\ / | / \ / \ ......................................................................