Stefan Fröse

Results 32 comments of Stefan Fröse

When you say sensitivity you mean in the gamma astronomy case? -> How large has to be the parameter x (could be flux) to reach $5\sigma$? This case is covered...

ok I see but I am not sure if this is correct for the $df>1$ case. You calculate the p value and than retroject this to a TS value on...

sry, I think this is also wrong

Thought about it again and came to this conclusion: We are interested in the median significance. So we need the median TS value $q_{med}$: 1. Solve $P(q_0\le q_{med})=0.5$ for q_med...

```python In [86]: ts_to_sigma(20, df=100) Out[86]: 10.919638828491388 In [87]: ts_to_sigma_old(20, df=100, asimov_dataset=True) Out[87]: -inf In [88]: ts_to_sigma(10, df=1) Out[88]: 3.162277660486695 In [89]: np.sqrt(10) Out[89]: 3.1622776601683795 In [90]: ts_to_sigma(10, df=100) Out[90]:...

the above ipython snippet is for: ```python def ts_to_sigma(ts, df=1): """Convert delta ts to number of sigma.""" q_med = ncx2.ppf(0.5, df=df, nc=ts) return np.sqrt(q_med) ``` and ```python def ts_to_sigma_old(ts, df=1,...

@QRemy sry I misunderstood :) I think the only way to invert with this function is to search for the root: ```python def sigma_to_ts(Z, df=1): """ Invert the Asimov‐median map...

@Astro-Kirsty @QRemy I think this is correct, I can open a PR if you also agree

Hi @lucascolley @henryiii 😊 I just stumbled on this since I was wondering if nox can do this 😄 Do you know if there have been any updates to pixi...

I am currently on my phone and can't check but maybe we need to check if the authors of CosmiXs use the same interpolation as we do. https://github.com/ajueid/CosmiXs/blob/main/Interpolate.py Also, they...