mt.benjamini_hochberg

benjamini_hochberg(maxima, alpha, plot=False)[source]

Select a subset of the maxima to be candidates to Point Source.

It applies the procedure in multiple testing called Benjamini-Hochberg procedure.

Parameters

  • maxima (pd.DataFrame) – Pandas DataFrame with the information of the maxima. Has to contain at least the column ‘pvalue’.

  • alpha (float) – Parameter alpha of the Benjamini-Hochberg procedure. Should be between 0 and 1.

  • plot (bool, optional) – If True, the function prints the expected number of sources for the intensities of the maxima versus the actual number of maxima at that intensity, along with the threshold (expected*alpha). Two plots are given: one with all the detections plus 5; and one with all the maxima.

Returns

The same Pandas DataFrame as the input maxima, but only with the candidates to be point sources.

Return type

pd.DataFrame