How to get uncertainty in peak location from HENzsearch

Open
#80 17 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Assessment

Difficulty
4/5
Estimated time
3-5 days
Newbie friendliness
25/100
Issue type
Feature
Clarity
Mostly clear
Activity status
Stale
Tech stack
jupyter-notebook
Domain
data

Research direction

Start with HENzsearch --fit-candidates and inspect how SincSquareModel results are produced and exposed through best_fits. Define how uncertainty for the fitted peak frequency should be returned alongside the existing parameters, then verify that the result is accessible in the shown workflow.

Written by the indexing model from the issue text.

Description

When I run HENzsearch with --fit-candidates, I see that I can get the parameters of the SincSquareModel by using

In [15]: efperiod = load_folding('cleanfilt_bary_nicer_xti_Z2n.nc')  
In [16]: efperiod.best_fits                                                                                                                               
Out[16]: [<SincSquareModel(amplitude=116321.52283363, mean=0.10777897, width=0.0000043)>]

This give the best fit frequency and the width of the sinc function. It would be nice if there were a way to easily get the uncertainty in the frequency determination. Can that be included somehow? Is that one of the outputs the fitter could return?

Dominant language
Jupyter Notebook
Stars
26
Forks
18
Avg merge
1d 20h
Merged PRs (30d)
8

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from StingraySoftware/HENDRICS

All issues in StingraySoftware/HENDRICS

Similar issues

More Data Engineering issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.