apriori.h

Go to the documentation of this file.

/**
 *
 *   Copyright (c) 2005-2021 by Pierre-Henri WUILLEMIN(@LIP6) & Christophe GONZALES(@AMU)
 *   info_at_agrum_dot_org
 *
 *  This library is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  This library is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with this library.  If not, see <http://www.gnu.org/licenses/>.
 *
 */


/** @file
 * @brief the base class for all a priori
 *
 * @author Christophe GONZALES(@AMU) and Pierre-Henri WUILLEMIN(@LIP6)
 */
#ifndef GUM_LEARNING_A_PRIORI_H
#define GUM_LEARNING_A_PRIORI_H

#include <string>
#include <vector>

#include <agrum/agrum.h>
#include <agrum/tools/core/bijection.h>
#include <agrum/BN/learning/aprioris/aprioriTypes.h>
#include <agrum/tools/database/databaseTable.h>
#include <agrum/tools/stattests/idCondSet.h>

namespace gum {

  namespace learning {

    /** @class Apriori
     * @brief the base class for all a priori
     * @headerfile apriori.h <agrum/tools/database/apriori.h>
     * @ingroup learning_apriori
     */
    template < template < typename > class ALLOC = std::allocator >
    class Apriori: private ALLOC< NodeId > {
      public:
      /// type for the allocators passed in arguments of methods
      using allocator_type = ALLOC< NodeId >;

      // ##########################################################################
      /// @name Constructors / Destructors
      // ##########################################################################
      /// @{

      /// default constructor
      /** @param database the database from which learning is performed. This is
       * useful to get access to the random variables
       * @param nodeId2Columns a mapping from the ids of the nodes in the
       * graphical model to the corresponding column in the DatabaseTable.
       * This enables estimating from a database in which variable A corresponds
       * to the 2nd column the parameters of a BN in which variable A has a
       * NodeId of 5. An empty nodeId2Columns bijection means that the mapping
       * is an identity, i.e., the value of a NodeId is equal to the index of
       * the column in the DatabaseTable.
       * @param alloc the allocator used to allocate the structures within the
       * RecordCounter.*/
      Apriori(const DatabaseTable< ALLOC >&                                 database,
              const Bijection< NodeId, std::size_t, ALLOC< std::size_t > >& nodeId2columns
              = Bijection< NodeId, std::size_t, ALLOC< std::size_t > >(),
              const allocator_type& alloc = allocator_type());

      /// virtual copy constructor
      virtual Apriori< ALLOC >* clone() const = 0;

      /// virtual copy constructor with a given allocator
      virtual Apriori< ALLOC >* clone(const allocator_type& alloc) const = 0;

      /// destructor
      virtual ~Apriori();

      /// @}

      // ##########################################################################
      /// @name Accessors / Modifiers
      // ##########################################################################
      /// @{

      /// sets the weight of the a priori (kind of effective sample size)
      virtual void setWeight(const double weight);

      /// returns the weight assigned to the apriori
      double weight() const;

      /// indicates whether an apriori is of a certain type
      virtual bool isOfType(const std::string& type) = 0;

      /// returns the type of the apriori
      virtual const std::string& getType() const = 0;

      /// indicates whether the apriori is potentially informative
      /** Basically, only the NoApriori is uninformative. However, it may happen
       * that, under some circonstances, an apriori, which is usually not equal
       * to the NoApriori, becomes equal to it (e.g., when the weight is equal
       * to zero). In this case, if the apriori can detect this case, it shall
       * inform the classes that use it that it is temporarily uninformative.
       * These classes will then be able to speed-up their code by avoiding to
       * take into account the apriori in their computations. */
      virtual bool isInformative() const = 0;

      /// adds the apriori to a counting vector corresponding to the idset
      /** adds the apriori to an already created counting vector defined over
       * the union of the variables on both the left and right hand side of the
       * conditioning bar of the idset.
       * @warning the method assumes that the size of the vector is exactly
       * the domain size of the joint variables set. */
      virtual void addAllApriori(const IdCondSet< ALLOC >&               idset,
                                 std::vector< double, ALLOC< double > >& counts)
         = 0;

      /** @brief adds the apriori to a counting vectordefined over the right
       * hand side of the idset
       *
       * @warning the method assumes that the size of the vector is exactly
       * the domain size of the joint RHS variables of the idset. */
      virtual void addConditioningApriori(const IdCondSet< ALLOC >&               idset,
                                          std::vector< double, ALLOC< double > >& counts)
         = 0;

      /// returns the allocator used by the internal apriori
      allocator_type getAllocator() const;

      /// @}


      protected:
      /// the weight of the apriori
      double weight_{1.0};

      /// a reference to the database in order to have access to its variables
      const DatabaseTable< ALLOC >* database_;

      /** @brief a mapping from the NodeIds of the variables to the indices of
       * the columns in the database */
      Bijection< NodeId, std::size_t, ALLOC< std::size_t > > nodeId2columns_;


      /// copy constructor
      Apriori(const Apriori< ALLOC >& from);

      /// copy constructor with a given allocator
      Apriori(const Apriori< ALLOC >& from, const allocator_type& alloc);

      /// move constructor
      Apriori(Apriori< ALLOC >&& from);

      /// move constructor with a given allocator
      Apriori(Apriori< ALLOC >&& from, const allocator_type& alloc);

      /// copy operator
      Apriori< ALLOC >& operator=(const Apriori< ALLOC >& from);

      /// move operator
      Apriori< ALLOC >& operator=(Apriori< ALLOC >&& from);
    };

  } /* namespace learning */

} /* namespace gum */

/// include the template implementation
#include <agrum/BN/learning/aprioris/apriori_tpl.h>

#endif /* GUM_LEARNING_A_PRIORI_H */