KernelDensityEstimation.h

Go to the documentation of this file.

// --------------------------------------------------------------------------
//                   OpenMS -- Open-Source Mass Spectrometry
// --------------------------------------------------------------------------
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Justin Sing $
// $Authors: Justin Sing $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <cstddef>
#include <utility>
#include <vector>
 
#include <OpenMS/config.h>
 
namespace OpenMS
{
  namespace Math
  {
    struct OPENMS_DLLAPI KernelDensityEstimation
    {
      static double bwNrd0(const std::vector<double>& x);
 
      static std::vector<double> linBin(const std::vector<double>& x,
                                        double xmin,
                                        double xmax,
                                        std::size_t nbins,
                                        const std::vector<double>* weights);
 
      static std::vector<double> linBin(const std::vector<double>& x,
                                        double xmin,
                                        double xmax,
                                        std::size_t nbins);
 
      static std::vector<double> forRt(const std::vector<double>& X, std::size_t M = 0);
 
      static std::vector<double> revRt(const std::vector<double>& Xp, std::size_t M = 0);
 
      static std::vector<double> silvermanKernelFFT(double bw, std::size_t M, double RANGE);
 
      static std::pair<std::vector<double>, std::vector<double>> gridKdeFFT(const std::vector<double>& x,
                                                                            double bw,
                                                                            std::size_t gridsize = 512,
                                                                            double cut = 3.0);
 
      static std::vector<double> kdeFFTEval(const std::vector<double>& x,
                                            double bw,
                                            std::size_t gridsize = 512,
                                            double cut = 3.0);
    };
 
  } // namespace Math
} // namespace OpenMS