ContingencyTableGenerator.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: The Bio++ Development Group
//
// SPDX-License-Identifier: CECILL-2.1

#ifndef BPP_NUMERIC_RANDOM_CONTINGENCYTABLEGENERATOR_H
#define BPP_NUMERIC_RANDOM_CONTINGENCYTABLEGENERATOR_H


#include "../Matrix/Matrix.h"
#include "RandomTools.h"

// From the STL:
#include <cmath>
#include <vector>

namespace bpp
{
class ContingencyTableGenerator
{
private:
  std::vector<size_t> nrowt_;
  std::vector<size_t> ncolt_;
  size_t nrow_;
  size_t ncol_;
  size_t nrowm_;
  size_t ncolm_;
  std::vector<size_t> jwork_; // workspace
  size_t ntot_; // total number of observations
  std::vector<double> fact_; // log factorial

public:
  ContingencyTableGenerator(const std::vector<size_t>& nrowt, const std::vector<size_t>& ncolt);

public:
  RowMatrix<size_t> rcont2();
};
} // end of namespace bpp.
#endif // BPP_NUMERIC_RANDOM_CONTINGENCYTABLEGENERATOR_H