[REDMINE1D-308] [RM-8275] Probleme prior N(z) ? Created: 13/Sep/23  Updated: 24/Oct/23  Resolved: 14/Oct/23

Status: Done
Project: 1D Redmine
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Task Priority: Normal
Reporter: Redmine-Jira Migtation Assignee: Redmine-Jira Migtation
Resolution: Done Votes: 0
Labels: None
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Attachments: PNG File Capture d’écran du 2023-10-12 14-59-13.png     PNG File clipboard-202310121628-znerx.png     PNG File clipboard-202310121651-zuhi5.png     PNG File clipboard-202310121652-qxoa1.png    

 Description   

Created on 2023-09-12 06:27:15 by Vincent Le Brun. % Done: 100

j'ai fait deux runs avec la 0.44 sur les 200000 spectres de EL-COSMOS, prior H-alpha, et avec 3 valeurs du prior euclidnhaemittersStrength différentes : -1, 1. et 0.001 et les résultats sont identiques (à quelques fluctuations près) et les PDF ne montrent aucune évolution...



 Comments   
Comment by Redmine-Jira Migtation [ 24/Oct/23 ]

Comment by Didier Vibert on 2023-09-12 06:46:03:
à regarder avant que @pychabau et @gdaste ne livre la 0.44 à EUCLID ?

Comment by Redmine-Jira Migtation [ 24/Oct/23 ]

Comment by Vincent Le Brun on 2023-09-12 06:50:49:
les répertoires :
/net/CESAM/amazed/vlebrun/EL-COSMOS2020/output_0.44
/net/CESAM/amazed/vlebrun/EL-COSMOS2020/output_0.44_Halpha_N_of_z_1
/net/CESAM/amazed/vlebrun/EL-COSMOS2020/output_0.44_Halpha_N_of_z_0.001

Comment by Redmine-Jira Migtation [ 24/Oct/23 ]

Comment by Didier Vibert on 2023-10-12 14:56:49:
si je prends un spectre au hazard (spec_1000043)

je vois une diff sur les pdf:
output_0.44 (prior Ha, pas de N(z)):

output_0.44_Halpha_N_of_z_1 (prior Ha, avec N(z)):

on voit clairement le prior N(z) pozzetti ajouté sur la logPDF ("fond" basse fréquence). Sur ce spectre ça ne change pas le best z, mais par exemple le candidat 4 à z~1.5 de la solution sans N(z) disparait au profit d'un candidat plus centré vers le max du prior.

un autre exemple: sans prior N(z):

avec N(z):

Les candidats sont aux mêmes z, mais leur ordre est modifié par le prior.

Du coup, où est le problème ? si sur les stats ça change pas, il faut sans doute recalibrer la force du prior... En tout cas, on n'a pas changé dans le code la façon dont le prior est ajouté, cependant les modifs des différentes release ont pu conduire à une variation du chi2 qui induise une recalibration du prior pour avoir le même effet.

@vlebrun , je mets en resolved ?

Comment by Redmine-Jira Migtation [ 24/Oct/23 ]

Comment by Vincent Le Brun on 2023-10-13 13:58:32:
Ok je ne savais pas etre très frais, les PDF sont bien modifiées mais on voit que le prior N(z) n'est pas assez fort pour contrer les solutions données par le prior H-alpha meme quand elles ne sont pas très fortes
on peut fermer

Generated at Sat Feb 10 15:31:24 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.