[Isis-fish-user] Calibration : valeurs de q1 négatives...et emplacement d'un individu faisant la taille limite de classe
Bonjour à tous j'ai lancé plusieurs simus de calibration, mais les résultats me semblent bien étranges. Pour les 3 premières simus, il me calcule le paramètre en utilisant les 3 valeurs de M1, M2 et M3 proposées dans le script. Par contre, pour les simus suivantes, il n'a pas vraiment l'air de prendre en compte le pas de calcul de 1e-05. Il propose des valeurs de q1 négatives....je ne comprends pas trop pourquoi (fichier Historic.csv ci-joint). Est-ce que qqun a déjà eu ces soucis et pourrait m'expliquer le pourquoi du comment ? Sinon je n'ai pas reçu d'avis concernant l'emplacement d'un individu faisant la taille limite de 2 classes de taille...est-ce qu'un individu de 60 cm va dans la classe 50-60 ou dans la 60-70 ?? Merci ! delphine
Delphine.Rocklin@ifremer.fr a écrit :
Bonjour à tous Bonjour, j'ai lancé plusieurs simus de calibration, mais les résultats me semblent bien étranges. Pour les 3 premières simus, il me calcule le paramètre en utilisant les 3 valeurs de M1, M2 et M3 proposées dans le script. Par contre, pour les simus suivantes, il n'a pas vraiment l'air de prendre en compte le pas de calcul de 1e-05. Il propose des valeurs de q1 négatives....je ne comprends pas trop pourquoi (fichier Historic.csv ci-joint).
Est-ce que qqun a déjà eu ces soucis et pourrait m'expliquer le pourquoi du comment ? Je ne vais pas trop pouvoir vous aider, mais j'ai regarder le plan d'expérience.
En fait, le paramètres "pas" n'est pas utilisé (on vérifiera la semaine prochaine pourquoi). A partir de la 4eme simulation, il utilise certaines valeurs des années précédentes. Au vu des opérations effectuées, il y a des soustraction, donc il est possible que cela deviennent négatif. Par contre, en effet si j'ai bien compris c'est valeurs représentent la capturabilité, et ne peuvent pas être négatives. -- Éric <chatellier@codelutin.com> Tel: 02 40 50 29 28 http://www.codelutin.com
Bonjour Mr Chatellier, ah oui en effet c'est étrange que le pas ne soit pas utilisé...et oui, des valeurs négatives pour la capturabilité, ce n'est pas très pratique... Si c'est arrivé à quelqu'un d'autre, je veux bien savoir comment vous avez règlé ce problème ! Merci, Delphine Eric Chatellier <chatellier@codelutin.com> a écrit :
Delphine.Rocklin@ifremer.fr a écrit :
Bonjour à tous Bonjour, j'ai lancé plusieurs simus de calibration, mais les résultats me semblent bien étranges. Pour les 3 premières simus, il me calcule le paramètre en utilisant les 3 valeurs de M1, M2 et M3 proposées dans le script. Par contre, pour les simus suivantes, il n'a pas vraiment l'air de prendre en compte le pas de calcul de 1e-05. Il propose des valeurs de q1 négatives....je ne comprends pas trop pourquoi (fichier Historic.csv ci-joint).
Est-ce que qqun a déjà eu ces soucis et pourrait m'expliquer le pourquoi du comment ? Je ne vais pas trop pouvoir vous aider, mais j'ai regarder le plan d'expérience.
En fait, le paramètres "pas" n'est pas utilisé (on vérifiera la semaine prochaine pourquoi).
A partir de la 4eme simulation, il utilise certaines valeurs des années précédentes. Au vu des opérations effectuées, il y a des soustraction, donc il est possible que cela deviennent négatif.
Par contre, en effet si j'ai bien compris c'est valeurs représentent la capturabilité, et ne peuvent pas être négatives.
-- Éric <chatellier@codelutin.com> Tel: 02 40 50 29 28 http://www.codelutin.com
_______________________________________________ Isis-fish-user mailing list Isis-fish-user@lists.labs.libre-entreprise.org http://lists.labs.libre-entreprise.org/mailman/listinfo/isis-fish-user
Salut à tous ! en effet le pas n est pas utilisé parce que la methode c est le simplex A PAS VARIABLE, il est donc calculé en fonction des resultats des simus precedantes. On avait laissé le parametre "param_pas" au cas ou un jour on voudrait faire un pas fixe. c est normal qu il tente des valeurs negatives, ce n est pas normal qu il converge vers des valeurs negatives. Pour l éviter il faut peut etre prendre des valeurs des points initiaux plus proches (car ca marche par symetrie...). Mais connaissant ton cas c est peut etre aussi qu il n a pas assez de points pour ajuster correctement... si ca marche mieux avec des valeurs negatives c est surement que la capturablité est tres faible : essaie des points initiaux vers 10^-10 ou 10^-7... sinon au pire tu rajoute une ligne "si q est negatif alors q = 0" tu verras bien s il converge vers 0. Pour les limites de classes j ai pas trop d avis là dessus... c est des hypotheses l important c est de savoir ce qu elles impliquent pour tes résultats. Bon courage ! A+ Sigrid Delphine.Rocklin@ifremer.fr a écrit :
Bonjour Mr Chatellier,
ah oui en effet c'est étrange que le pas ne soit pas utilisé...et oui, des valeurs négatives pour la capturabilité, ce n'est pas très pratique...
Si c'est arrivé à quelqu'un d'autre, je veux bien savoir comment vous avez règlé ce problème !
Merci, Delphine
Eric Chatellier <chatellier@codelutin.com> a écrit :
Delphine.Rocklin@ifremer.fr a écrit :
Bonjour à tous Bonjour, j'ai lancé plusieurs simus de calibration, mais les résultats me semblent bien étranges. Pour les 3 premières simus, il me calcule le paramètre en utilisant les 3 valeurs de M1, M2 et M3 proposées dans le script. Par contre, pour les simus suivantes, il n'a pas vraiment l'air de prendre en compte le pas de calcul de 1e-05. Il propose des valeurs de q1 négatives....je ne comprends pas trop pourquoi (fichier Historic.csv ci-joint).
Est-ce que qqun a déjà eu ces soucis et pourrait m'expliquer le pourquoi du comment ? Je ne vais pas trop pouvoir vous aider, mais j'ai regarder le plan d'expérience.
En fait, le paramètres "pas" n'est pas utilisé (on vérifiera la semaine prochaine pourquoi).
A partir de la 4eme simulation, il utilise certaines valeurs des années précédentes. Au vu des opérations effectuées, il y a des soustraction, donc il est possible que cela deviennent négatif.
Par contre, en effet si j'ai bien compris c'est valeurs représentent la capturabilité, et ne peuvent pas être négatives.
-- Éric <chatellier@codelutin.com> Tel: 02 40 50 29 28 http://www.codelutin.com
_______________________________________________ Isis-fish-user mailing list Isis-fish-user@lists.labs.libre-entreprise.org http://lists.labs.libre-entreprise.org/mailman/listinfo/isis-fish-user
_______________________________________________ Isis-fish-user mailing list Isis-fish-user@lists.labs.libre-entreprise.org http://lists.labs.libre-entreprise.org/mailman/listinfo/isis-fish-user
-- Sigrid LEHUTA ~ ><> ~ Doctorante Département Ecologie et Modèles pour l'Halieutique IFREMER, rue de l'ile d'Yeu BP 21105 44311 Nantes Cedex 03 Tél : +33 (0)2 40 37 41 65 (interne : 8165)
salut, pour éviter ce genre de problème, quand on a des valeurs très faibles et pas trop d'idée de la puissance de 10 (10^-10 ou 10^-7), on a parfois intérêt à faire une petite transformation du paramètre à estimer. Dans ce cas, estimer log(capturabilité), au lieu d'estimer directement la capturabilité, devrait faciliter le travail de l'algorithme d'optimisation. Reste à modifier un peu le script de calibration mais ça doit pas être trop compliquer: faut juste qu'à l'endroit où tu faisais une transformation inverse (exponentielle) au moment où tu fixes la valeur. En espérant que ça serve A+ Hilaire Le 24/08/2009 10:36, Sigrid LEHUTA a écrit :
Salut à tous !
en effet le pas n est pas utilisé parce que la methode c est le simplex A PAS VARIABLE, il est donc calculé en fonction des resultats des simus precedantes. On avait laissé le parametre "param_pas" au cas ou un jour on voudrait faire un pas fixe.
c est normal qu il tente des valeurs negatives, ce n est pas normal qu il converge vers des valeurs negatives. Pour l éviter il faut peut etre prendre des valeurs des points initiaux plus proches (car ca marche par symetrie...). Mais connaissant ton cas c est peut etre aussi qu il n a pas assez de points pour ajuster correctement... si ca marche mieux avec des valeurs negatives c est surement que la capturablité est tres faible : essaie des points initiaux vers 10^-10 ou 10^-7... sinon au pire tu rajoute une ligne "si q est negatif alors q = 0" tu verras bien s il converge vers 0.
Pour les limites de classes j ai pas trop d avis là dessus... c est des hypotheses l important c est de savoir ce qu elles impliquent pour tes résultats.
Bon courage ! A+ Sigrid
Delphine.Rocklin@ifremer.fr a écrit :
Bonjour Mr Chatellier,
ah oui en effet c'est étrange que le pas ne soit pas utilisé...et oui, des valeurs négatives pour la capturabilité, ce n'est pas très pratique...
Si c'est arrivé à quelqu'un d'autre, je veux bien savoir comment vous avez règlé ce problème !
Merci, Delphine
Eric Chatellier <chatellier@codelutin.com> a écrit :
Delphine.Rocklin@ifremer.fr a écrit :
Bonjour à tous Bonjour, j'ai lancé plusieurs simus de calibration, mais les résultats me semblent bien étranges. Pour les 3 premières simus, il me calcule le paramètre en utilisant les 3 valeurs de M1, M2 et M3 proposées dans le script. Par contre, pour les simus suivantes, il n'a pas vraiment l'air de prendre en compte le pas de calcul de 1e-05. Il propose des valeurs de q1 négatives....je ne comprends pas trop pourquoi (fichier Historic.csv ci-joint).
Est-ce que qqun a déjà eu ces soucis et pourrait m'expliquer le pourquoi du comment ? Je ne vais pas trop pouvoir vous aider, mais j'ai regarder le plan d'expérience.
En fait, le paramètres "pas" n'est pas utilisé (on vérifiera la semaine prochaine pourquoi).
A partir de la 4eme simulation, il utilise certaines valeurs des années précédentes. Au vu des opérations effectuées, il y a des soustraction, donc il est possible que cela deviennent négatif.
Par contre, en effet si j'ai bien compris c'est valeurs représentent la capturabilité, et ne peuvent pas être négatives.
-- Éric <chatellier@codelutin.com> Tel: 02 40 50 29 28 http://www.codelutin.com
_______________________________________________ Isis-fish-user mailing list Isis-fish-user@lists.labs.libre-entreprise.org http://lists.labs.libre-entreprise.org/mailman/listinfo/isis-fish-user
_______________________________________________ Isis-fish-user mailing list Isis-fish-user@lists.labs.libre-entreprise.org http://lists.labs.libre-entreprise.org/mailman/listinfo/isis-fish-user
participants (4)
-
Delphine.Rocklin＠ifremer.fr -
Eric Chatellier -
hilaire.drouineau -
Sigrid LEHUTA