nag_correg_linregm_coeffs_noconst (g02ch) performs a multiple linear regression with no constant on a set of variables whose sums of squares and cross-products about zero and correlation-like coefficients are given.

Syntax

[result, coef, rznv, cz, ifail] = g02ch(n, sspz, rz, 'k1', k1)

[result, coef, rznv, cz, ifail] = nag_correg_linregm_coeffs_noconst(n, sspz, rz, 'k1', k1)

Note: the interface to this routine has changed since earlier releases of the toolbox:

At Mark 23:

k was removed from the interface

Description

nag_correg_linregm_coeffs_noconst (g02ch) fits a curve of the form

y = b_{1} x_{1} + b_{2} x_{2} + \dots + b_{k} x_{k}

to the data points

\begin{matrix} (x_{11}, x_{21}, \dots, x_{k 1}, y_{1}) \\ (x_{12}, x_{22}, \dots, x_{k 2}, y_{2}) \\ ⋮ \\ (x_{1 n}, x_{2 n}, \dots, x_{k n}, y_{n}) \end{matrix}

such that

y_{i} = b_{1} x_{1 i} + b_{2} x_{2 i} + \dots + b_{k} x_{k i} + e_{i}, i = 1, 2, \dots, n .

The function calculates the regression coefficients,

b_{1}, b_{2}, \dots, b_{k}

, (and various other statistical quantities) by minimizing

\sum_{i = 1}^{n} e_{i}^{2} .

The actual data values

(x_{1 i}, x_{2 i}, \dots, x_{k i}, y_{i})

are not provided as input to the function. Instead, input to the function consists of:

(i)	The number of cases, $n$ , on which the regression is based.
(ii)	The total number of variables, dependent and independent, in the regression, $(k + 1)$ .
(iii)	The number of independent variables in the regression, $k$ .
(iv)	The $(k + 1)$ by $(k + 1)$ matrix $[{\tilde{S}}_{i j}]$ of sums of squares and cross-products about zero of all the variables in the regression; the terms involving the dependent variable, $y$ , appear in the $(k + 1)$ th row and column.
(v)	The $(k + 1)$ by $(k + 1)$ matrix $[{\tilde{R}}_{i j}]$ of correlation-like coefficients for all the variables in the regression; the correlations involving the dependent variable, $y$ , appear in the $(k + 1)$ th row and column.

The quantities calculated are:

(a)

The inverse of the

k

k

partition of the matrix of correlation-like coefficients,

[{\tilde{R}}_{i j}]

, involving only the independent variables. The inverse is obtained using an accurate method which assumes that this sub-matrix is positive definite (see Further Comments).

(b)

The modified matrix,

C = [c_{i j}]

, where

c_{i j} = \frac{{\tilde{R}}_{i j} {\tilde{r}}^{i j}}{{\tilde{S}}_{i j}}, i, j = 1, 2, \dots, k,

where

{\tilde{r}}^{i j}

is the

(i, j)

th element of the inverse matrix of

[{\tilde{R}}_{i j}]

as described in (a) above. Each element of

C

is thus the corresponding element of the matrix of correlation-like coefficients multiplied by the corresponding element of the inverse of this matrix, divided by the corresponding element of the matrix of sums of squares and cross-products about zero.

(c)

The regression coefficients:

b_{i} = \sum_{j = 1}^{k} c_{i j} {\tilde{S}}_{j (k + 1)}, i = 1, 2, \dots, k,

where

{\tilde{S}}_{j (k + 1)}

is the sum of cross-products about zero for the independent variable

x_{j}

and the dependent variable

y

(d)

The sum of squares attributable to the regression,

S S R

, the sum of squares of deviations about the regression,

S S D

, and the total sum of squares,

S S T

$S S T = {\tilde{S}}_{(k + 1) (k + 1)}$ , the sum of squares about zero for the dependent variable, $y$ ;
$S S R = \sum_{j = 1}^{k} b_{j} {\tilde{S}}_{j (k + 1)}; S S D = S S T - S S R$ .

(e)

The degrees of freedom attributable to the regression,

D F R

, the degrees of freedom of deviations about the regression,

D F D

, and the total degrees of freedom,

D F T

D F R = k; D F D = n - k; D F T = n .

(f)

The mean square attributable to the regression,

M S R

, and the mean square of deviations about the regression,

M S D

M S R = S S R / D F R; M S D = S S D / D F D .

(g)

The

F

value for the analysis of variance:

F = M S R / M S D .

(h)

The standard error estimate:

s = \sqrt{M S D} .

(i)

The coefficient of multiple correlation,

R

, the coefficient of multiple determination,

R^{2}

, and the coefficient of multiple determination corrected for the degrees of freedom,

{\bar{R}}^{2}

R = \sqrt{1 - \frac{S S D}{S S T}}; R^{2} = 1 - \frac{S S D}{S S T}; {\bar{R}}^{2} = 1 - \frac{S S D \times D F T}{S S T \times D F D} .

(j)

The standard error of the regression coefficients:

s e (b_{i}) = \sqrt{M S D \times c_{i i}}, i = 1, 2, \dots, k .

(k)

The

t

values for the regression coefficients:

t (b_{i}) = \frac{b_{i}}{s e (b_{i})}, i = 1, 2, \dots, k .

References

Draper N R and Smith H (1985) Applied Regression Analysis (2nd Edition) Wiley

Parameters

Compulsory Input Parameters

1: $n$ – int64int32nag_int scalar: $n$ , the number of cases used in calculating the sums of squares and cross-products and correlation-like coefficients.
2: $sspz (ldsspz, k1)$ – double array: ldsspz, the first dimension of the array, must satisfy the constraint $ldsspz \geq k1$ .
$sspz (i, j)$ must be set to ${\tilde{S}}_{i j}$ , the sum of cross-products about zero for the $i$ th and $j$ th variables, for $i = 1, 2, \dots, k + 1$ and $j = 1, 2, \dots, k + 1$ ; terms involving the dependent variable appear in row $k + 1$ and column $k + 1$ .
3: $rz (ldrz, k1)$ – double array: ldrz, the first dimension of the array, must satisfy the constraint $ldrz \geq k1$ .
$rz (i, j)$ must be set to ${\tilde{R}}_{i j}$ , the correlation-like coefficient for the $i$ th and $j$ th variables, for $i = 1, 2, \dots, k + 1$ and $j = 1, 2, \dots, k + 1$ ; coefficients involving the dependent variable appear in row $k + 1$ and column $k + 1$ .

Optional Input Parameters

1: $k1$ – int64int32nag_int scalar: Default: the first dimension of the arrays sspz, rz and the second dimension of the arrays sspz, rz. (An error is raised if these dimensions are not equal.)
The total number of variables, independent and dependent $(k + 1)$ , in the regression.

Constraint: $2 \leq k1 \leq n$ .

Output Parameters

1: $result (13)$ – double array

The following information:

$result (1)$	$S S R$ , the sum of squares attributable to the regression;
$result (2)$	$D F R$ , the degrees of freedom attributable to the regression;
$result (3)$	$M S R$ , the mean square attributable to the regression;
$result (4)$	$F$ , the $F$ value for the analysis of variance;
$result (5)$	$S S D$ , the sum of squares of deviations about the regression;
$result (6)$	$D F D$ , the degrees of freedom of deviations about the regression;
$result (7)$	$M S D$ , the mean square of deviations about the regression;
$result (8)$	$S S T$ , the total sum of squares;
$result (9)$	$D F T$ , the total degrees of freedom;
$result (10)$	$s$ , the standard error estimate;
$result (11)$	$R$ , the coefficient of multiple correlation;
$result (12)$	$R^{2}$ , the coefficient of multiple determination;
$result (13)$	${\bar{R}}^{2}$ , the coefficient of multiple determination corrected for the degrees of freedom.

2: $coef (ldcoef, 3)$ – double array

For

i = 1, 2, \dots, k

, the following information:

$coef (i, 1)$: $b_{i}$ , the regression coefficient for the $i$ th variable.
$coef (i, 2)$: $s e (b_{i})$ , the standard error of the regression coefficient for the $i$ th variable.
$coef (i, 3)$: $t (b_{i})$ , the $t$ value of the regression coefficient for the $i$ th variable.

3: $rznv (ldrznv, k)$ – double array

k = k1 - 1

The inverse of the matrix of correlation-like coefficients for the independent variables; that is, the inverse of the matrix consisting of the first

k

rows and columns of rz.

4: $cz (ldcz, k)$ – double array

k = k1 - 1

The modified inverse matrix,

C

, where

cz (i, j) = \frac{rz (i, j) \times rznv (i, j)}{sspz (i, j)}, i, j = 1, 2, \dots, k .

5: $ifail$ – int64int32nag_int scalar

ifail = 0

unless the function detects an error (see Error Indicators and Warnings).

Error Indicators and Warnings

Errors or warnings detected by the function:

$ifail = 1$

On entry,

k1 < 2

$ifail = 2$

On entry,

k1 \neq (k + 1)

$ifail = 3$

On entry,

n < k1

$ifail = 4$

On entry,	$ldsspz < k1$ ,
or	$ldrz < k1$ ,
or	$ldcoef < k$ ,
or	$ldrznv < k$ ,
or	$ldcz < k$ ,
or	$ldwkz < k$ .

$ifail = 5$: This indicates that the $k$ by $k$ partition of the matrix held in rz, which is to be inverted, is not positive definite.

$ifail = 6$: This indicates that the refinement following the actual inversion fails, indicating that the $k$ by $k$ partition of the matrix held in rz, which is to be inverted, is ill-conditioned. The use of nag_correg_linregm_fit (g02da), which employs a different numerical technique, may avoid the difficulty.

$ifail = 7$: Unexpected error in nag_linsys_real_posdef_solve_ref (f04ab).

$ifail = - 99$: An unexpected error has been triggered by this routine. Please contact NAG.

$ifail = - 399$: Your licence key may have expired or may not have been installed correctly.

$ifail = - 999$: Dynamic memory allocation failed.

Accuracy

The accuracy of any regression function is almost entirely dependent on the accuracy of the matrix inversion method used. In nag_correg_linregm_coeffs_noconst (g02ch), it is the matrix of correlation-like coefficients rather than that of the sums of squares and cross-products about zero that is inverted; this means that all terms in the matrix for inversion are of a similar order, and reduces the scope for computational error. For details on absolute accuracy, the relevant section of the document describing the inversion function used, nag_linsys_real_posdef_solve_ref (f04ab), should be consulted. nag_correg_linregm_fit (g02da) uses a different method, based on nag_linsys_real_gen_lsqsol (f04am), and that function may well prove more reliable numerically. It does not handle missing values, nor does it provide the same output as this function.

If, in calculating

F

or any of the

t (b_{i})

(see Description), the numbers involved are such that the result would be outside the range of numbers which can be stored by the machine, then the answer is set to the largest quantity which can be stored as a double variable, by means of a call to nag_machine_real_largest (x02al).

Further Comments

The time taken by nag_correg_linregm_coeffs_noconst (g02ch) depends on

k

This function assumes that the matrix of correlation-like coefficients for the independent variables in the regression is positive definite; it fails if this is not the case.

This correlation matrix will in fact be positive definite whenever the correlation-like matrix and the sums of squares and cross-products (about zero) matrix have been formed either without regard to missing values, or by eliminating completely any cases involving missing values for any variable. If, however, these matrices are formed by eliminating cases with missing values from only those calculations involving the variables for which the values are missing, no such statement can be made, and the correlation-like matrix may or may not be positive definite. You should be aware of the possible dangers of using correlation matrices formed in this way (see the G02 Chapter Introduction), but if they nevertheless wish to carry out regressions using such matrices, this function is capable of handling the inversion of such matrices, provided they are positive definite.

If a matrix is positive definite, its subsequent re-organisation by either of nag_correg_linregm_service_select (g02ce) or nag_correg_linregm_service_reorder (g02cf) will not affect this property and the new matrix can safely be used in this function. Thus correlation matrices produced by any of nag_correg_coeffs_zero (g02bd), nag_correg_coeffs_zero_miss_case (g02be), nag_correg_coeffs_zero_subset (g02bk) or nag_correg_coeffs_zero_subset_miss_case (g02bl), even if subsequently modified by either nag_correg_linregm_service_select (g02ce) or nag_correg_linregm_service_reorder (g02cf), can be handled by this function.

It should be noted that the function requires the dependent variable to be the last of the

k + 1

variables whose statistics are provided as input to the function. If this variable is not correctly positioned in the original data, the means, standard deviations, sums of squares and cross-products about zero, and correlation-like coefficients can be manipulated by using nag_correg_linregm_service_select (g02ce) or nag_correg_linregm_service_reorder (g02cf) to reorder the variables as necessary.

Example

This example reads in the sums of squares and cross-products about zero, and correlation-like coefficients for three variables. A multiple linear regression with no constant is then performed with the third and final variable as the dependent variable. Finally the results are printed.

Open in the MATLAB editor: g02ch_example

function g02ch_example


fprintf('g02ch example results\n\n');

n = int64(5);
k = 2;
sspz = [245,  99, 82;
         99, 271, 52;
         82,  52, 54];
rz   = [1,      0.3842, 0.7129;
        0.3842, 1,      0.4299;
        0.7129, 0.4299, 1     ];

disp('Sums of squares and cross-products about zero:');
disp(sspz);
disp('Correlation-like coefficients:');
disp(rz);

% Fit multiple linear regression model
[result, coeff, rzinv, cz, ifail] = ...
  g02ch(n, sspz, rz);

% Display results
fprintf('  Variable    Coef      Std err  t-value\n');
disp([[1:k]' coeff]);
  
fprintf('Analysis of regression table :-\n\n');
fprintf('      Source        Sum of squares  DF   Mean square   F-value\n');
fprintf('Due to regression %11.3f%8d%13.3f%12.3f\n', result(1:4));
fprintf('About  regression %11.3f%8d%13.3f\n', result(5:7));
fprintf('Total             %11.3f%8d\n\n', result(8:9));
fprintf('Standard error of estimate = %8.4f\n', result(10));
fprintf('Multiple correlation (R)   = %8.4f\n', result(11));
fprintf('Determination (R squared)  = %8.4f\n', result(12));
fprintf('Corrected R squared        = %8.4f\n\n', result(13));

disp('Inverse of correlation matrix of independent variables:');
disp(rzinv);
disp('Modified inverse matrix:');
disp(cz);

g02ch example results

Sums of squares and cross-products about zero:
   245    99    82
    99   271    52
    82    52    54

Correlation-like coefficients:
    1.0000    0.3842    0.7129
    0.3842    1.0000    0.4299
    0.7129    0.4299    1.0000

  Variable    Coef      Std err  t-value
    1.0000    0.3017    0.1998    1.5098
    2.0000    0.0817    0.1900    0.4299

Analysis of regression table :-

      Source        Sum of squares  DF   Mean square   F-value
Due to regression      28.986       2       14.493       1.738
About  regression      25.014       3        8.338
Total                  54.000       5

Standard error of estimate =   2.8876
Multiple correlation (R)   =   0.7326
Determination (R squared)  =   0.5368
Corrected R squared        =   0.2280

Inverse of correlation matrix of independent variables:
    1.1732   -0.4507
   -0.4507    1.1732

Modified inverse matrix:
    0.0048   -0.0017
   -0.0017    0.0043

PDF version (NAG web site, 64-bit version, 64-bit version)

Chapter Contents

Chapter Introduction

NAG Toolbox