Hello Bastien, Oui oui t as bien compris !
Lors de la calibration de la capturabilité dans ISIS, il y a q1 et q2 comme facteurs (c’est donc un simplexe à 2 facteurs ; avec trois couples de facteurs, initialisés avec M1, M2 et M3), c’est bien cela ?
Yep ! comme tu l as vu dans le walters c est geometrique donc attention a former un joli triangle de depart avec une distance pertinente entre les points initiaux : si c est trop grand il risque de partir rapidement en negatif, si c est trop petit c est moins grave mais plus long.... L ideal c est de tester plusieurs simplex initiaux pour verifier qu on est dans tombé dans des minimums locaux.
Par contre, je ne suis pas sure de bien comprendre à quoi correspondent q1 et q2… ? C’est une décomposition de la capturabilité q ?
Oui c est ca. En fait, dans isis, la capturabilité peut dependre de la classe et de la saison. Ca fait potentiellement bcp de valeurs à estimer. pour des probleme d identifiabilité et de pertinence on a fait des hypotheses, on imagine que certaines valeurs sont identiques et on essaie de garder seulement quelques (2 ou 3) coeff a estimer. Il existe je crois des simplexes de plus grande dimenssion et aussi d autres methodes de calibration qui permettent d estimer un plus grands nombre de parametres (algo genetiques par ex) mais il faut les coder.
Dans le script, la fonction objectif, c’est bien le critère ? c'est-à-dire le carré de l’écart entre observé et simulé (ligne 328) ? Mais pourquoi crit + (obs-simu)^2 ?
Oui c est une betise de l appeler critere c est bien la fonction objectif. Premiere chose, il n y a aucune obligation a calibrer sur les captures et encore moins sur les captures par age, c est selon la pertinence les données dispo... ensuite on a choisit comme fonction la somme sur les dates et les ages de : captures observées - captures simulees on a donc d abord creer la matrice captures simulées par age et dates (les captures observees sont au meme format): ensuite avec le matrix iterator on se promene dans la matrice et a chaque case on calcule (obs-sim)^2 et on l ajoute au total deja calculé (crit). Comme ca on a la valeur de la fonction objectif sur toute la duree de simulation et toutes les classes ... on pourrait choisir de ponderer aussi des dates ou des classes... Ca repond a ta question ? A++
Merci.
Bastien
Bastien Preuss
Doctorant (PhD Student)
IFREMER - IRD - Université de Nouvelle-Calédonie
101 Promenade Roger Laroque
BP A5
98848 Noumea cedex
New-Caledonia
Tel: (+687) 260 723
P_ __afin de contribuer au respect de l'environnement, merci de n'imprimer ce mail qu'en cas de nécessité_
------------------------------------------------------------------------
_______________________________________________ Isis-fish-users mailing list Isis-fish-users@list.isis-fish.org http://list.isis-fish.org/cgi-bin/mailman/listinfo/isis-fish-users
-- 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 23 (interne : 8123)