Package org.apache.commons.math4.legacy.stat.regression

Class AbstractMultipleLinearRegression

    • Method Detail

      • isNoIntercept

        public boolean isNoIntercept()
        Returns:
        true if the model has no intercept term; false otherwise
        Since:
        2.2

      • setNoIntercept

        public void setNoIntercept​(boolean noIntercept)
        Parameters:
        noIntercept - true means the model is to be estimated without an intercept term
        Since:
        2.2

      • newSampleData

        public void newSampleData​(double[] data,
                          int nobs,
                          int nvars)

        Loads model x and y sample data from a flat input array, overriding any previous sample.

        Assumes that rows are concatenated with y values first in each row. For example, an input data array containing the sequence of values (1, 2, 3, 4, 5, 6, 7, 8, 9) with nobs = 3 and nvars = 2 creates a regression dataset with two independent variables, as below: 

           y   x[0]  x[1]
   --------------
   1     2     3
   4     5     6
   7     8     9

        Note that there is no need to add an initial unitary column (column of 1's) when specifying a model including an intercept term. If isNoIntercept() is true, the X matrix will be created without an initial column of "1"s; otherwise this column will be added.

        Throws IllegalArgumentException if any of the following preconditions fail:

        • data cannot be null
        • data.length = nobs * (nvars + 1)
        • nobs > nvars
        Parameters:
        data - input data array
        nobs - number of observations (rows)
        nvars - number of independent variables (columns, not counting y)
        Throws:
        NullArgumentException - if the data array is null
        DimensionMismatchException - if the length of the data array is not equal to nobs * (nvars + 1)
        InsufficientDataException - if nobs is less than nvars + 1

      • newXSampleData

        protected void newXSampleData​(double[][] x)

        Loads new x sample data, overriding any previous data.

        The input x array should have one row for each sample observation, with columns corresponding to independent variables. For example, if 
          x = new double[][] {{1, 2}, {3, 4}, {5, 6}}
        then setXSampleData(x) results in a model with two independent variables and 3 observations: 
           x[0]  x[1]
   ----------
     1    2
     3    4
     5    6

        Note that there is no need to add an initial unitary column (column of 1's) when specifying a model including an intercept term.

        Parameters:
        x - the rectangular array representing the x sample
        Throws:
        NullArgumentException - if x is null
        NoDataException - if x is empty
        DimensionMismatchException - if x is not rectangular

      • validateSampleData

        protected void validateSampleData​(double[][] x,
                                  double[] y)
                           throws MathIllegalArgumentException
        Validates sample data. Checks that
        • Neither x nor y is null or empty;
        • The length (i.e. number of rows) of x equals the length of y
        • x has at least one more row than it has columns (i.e. there is sufficient data to estimate regression coefficients for each of the columns in x plus an intercept.
        Parameters:
        x - the [n,k] array representing the x data
        y - the [n,1] array representing the y data
        Throws:
        NullArgumentException - if x or y is null
        DimensionMismatchException - if x and y do not have the same length
        NoDataException - if x or y are zero-length
        MathIllegalArgumentException - if the number of rows of x is not larger than the number of columns + 1

      • validateCovarianceData

        protected void validateCovarianceData​(double[][] x,
                                      double[][] covariance)
        Validates that the x data and covariance matrix have the same number of rows and that the covariance matrix is square.
        Parameters:
        x - the [n,k] array representing the x sample
        covariance - the [n,n] array representing the covariance matrix
        Throws:
        DimensionMismatchException - if the number of rows in x is not equal to the number of rows in covariance
        NonSquareMatrixException - if the covariance matrix is not square

      • estimateErrorVariance

        public double estimateErrorVariance()
        Estimates the variance of the error.
        Returns:
        estimate of the error variance
        Since:
        2.2

      • estimateRegressionStandardError

        public double estimateRegressionStandardError()
        Estimates the standard error of the regression.
        Returns:
        regression standard error
        Since:
        2.2

      • calculateBeta

        protected abstract RealVector calculateBeta()
        Calculates the beta of multiple linear regression in matrix notation.
        Returns:
        beta

      • calculateBetaVariance

        protected abstract RealMatrix calculateBetaVariance()
        Calculates the beta variance of multiple linear regression in matrix notation.
        Returns:
        beta variance

      • calculateYVariance

        protected double calculateYVariance()
        Calculates the variance of the y values.
        Returns:
        Y variance

      • calculateErrorVariance

        protected double calculateErrorVariance()

        Calculates the variance of the error term.

        Uses the formula 
         var(u) = u · u / (n - k)
        where n and k are the row and column dimensions of the design matrix X.
        Returns:
        error variance estimate
        Since:
        2.2

      • calculateResiduals

        protected RealVector calculateResiduals()
        Calculates the residuals of multiple linear regression in matrix notation. 
         u = y - X * b
        Returns:
        The residuals [n,1] matrix