.. _program_listing_file_include_gwmodelpp_GWRBase.h:

Program Listing for File GWRBase.h
==================================

|exhale_lsh| :ref:`Return to documentation for file <file_include_gwmodelpp_GWRBase.h>` (``include/gwmodelpp/GWRBase.h``)

.. |exhale_lsh| unicode:: U+021B0 .. UPWARDS ARROW WITH TIP LEFTWARDS

.. code-block:: cpp

   #ifndef GWRBASE_H
   #define GWRBASE_H
   
   #include "SpatialMonoscaleAlgorithm.h"
   #include "IRegressionAnalysis.h"
   
   namespace gwm
   {
   
   class GWRBase : public SpatialMonoscaleAlgorithm, public IRegressionAnalysis
   {
   public:
   
       static arma::vec Fitted(const arma::mat& x, const arma::mat& betas)
       {
           return sum(betas % x, 1);
       }
   
       static double RSS(const arma::mat& x, const arma::mat& y, const arma::mat& betas)
       {
           arma::vec r = y - Fitted(x, betas);
           return sum(r % r);
       }
   
       static double AICc(const arma::mat& x, const arma::mat& y, const arma::mat& betas, const arma::vec& shat)
       {
           double ss = RSS(x, y, betas), n = (double)x.n_rows;
           return n * log(ss / n) + n * log(2 * arma::datum::pi) + n * ((n + shat(0)) / (n - 2 - shat(0)));
       }
   
   public:
   
       GWRBase() {}
   
       GWRBase(const arma::mat& x, const arma::vec& y, const SpatialWeight& spatialWeight, const arma::mat& coords) : SpatialMonoscaleAlgorithm(spatialWeight, coords)
       {
           mX = x;
           mY = y;
       }
   
       ~GWRBase() {}
   
   public:
   
       const arma::mat& betas() const { return mBetas; }
   
   public:     // Implement IRegressionAnalysis
       virtual const arma::vec& dependentVariable() const override { return mY; }
       virtual void setDependentVariable(const arma::vec& y) override { mY = y; }
   
       virtual const arma::mat& independentVariables() const override { return mX; }
       virtual void setIndependentVariables(const arma::mat& x) override { mX = x; }
   
       virtual bool hasIntercept() const override { return mHasIntercept; }
       virtual void setHasIntercept(const bool has) override { mHasIntercept = has; }
   
       virtual RegressionDiagnostic diagnostic() const override { return mDiagnostic; }
   
   public:
       virtual bool isValid() override;
   
   protected:
   
       arma::mat mX;   
       arma::vec mY;   
       arma::mat mBetas;   
       bool mHasIntercept = true;  
   
       RegressionDiagnostic mDiagnostic;   
   };
   
   }
   
   #endif  // GWRBASE_H