MODULE II: CORRELATION AND REGRESSION ANALYSIS
1. FOUNDATIONAL CONCEPTS OF BIVARIATE RELATIONSHIPS
In corporate and economic systems, variables rarely operate in isolation. Managerial decisions depend critically upon understanding how alterations in one market factor correspond with variations in another. When an investigation simultaneously examines two quantifiable characteristics across identical statistical entities, the distribution is characterized as a Bivariate Frequency Distribution.
Concept of Co-Variation
Causation vs. Covariation
Taxonomy of Correlation
Statistical theory categorizes correlation across three distinct dimensions based on direction, linearity, and number of variables under examination:
| Classification Basis | Sub-Categories | Mathematical Characterization & Commercial Illustration |
|---|---|---|
| 1. Direction of Change |
Positive (Direct) Negative (Inverse) |
• Positive: Both variables move in the identical direction ((Delta X > 0 implies Delta Y > 0)). E.g., Advertising expenditure and product sales revenue; Disposable income and consumer expenditure. • Negative: Variables move in opposite directions ((Delta X > 0 implies Delta Y < 0)). E.g., Unit price and market quantity demanded; Production volume and per-unit fixed overhead cost. |
| 2. Ratio of Variation |
Linear Correlation Curvilinear (Non-linear) |
• Linear: The amount of change in one variable bears a strictly constant ratio to the change in the other variable. Plotted on a Cartesian graph, points cluster along a straight line ((Y = a + bX)). • Curvilinear: The ratio of variation is not constant, altering at varying levels of the independent variable. Plotted points trace a curve (e.g., diminishing marginal returns curve, production cost curves). |
| 3. Number of Variables |
Simple Correlation Partial Correlation Multiple Correlation |
• Simple: Studies mutual co-variation between exactly two variables ((X) and (Y)). • Partial: Measures relationship between two variables while holding one or more interrelated variables mathematically constant (denoted (r_{12.3})). • Multiple: Investigates the joint relationship between one dependent variable and the simultaneous combined effect of two or more independent variables (denoted (R_{1.23})). |
2. SCATTER DIAGRAM METHOD
The Scatter Diagram is a visual, diagrammatic device for inspecting bivariate data. Values of the independent or predictor variable ((X)) are plotted along the horizontal (X)-axis, while the corresponding values of the dependent or criterion variable ((Y)) are plotted along the vertical (Y)-axis. Each pair of observations ((X_i, Y_i)) is plotted as a discrete coordinate dot.
1. Perfect Positive Correlation ((r = +1.0))
2. Perfect Negative Correlation ((r = -1.0))
3. High Degree Positive / Negative
4. Zero / Absence of Correlation ((r = 0))
3. KARL PEARSON'S COEFFICIENT OF CORRELATION ((r))
Introduced by the English biometrician Karl Pearson, the Product-Moment Correlation Coefficient ((r)) is the most widely utilized parametric measure of linear association between two quantitative variables. It is defined mathematically as the ratio of the covariance between (X) and (Y) to the product of their respective standard deviations.
Alternative Operational Formulas for Computation
In applied business statistics, direct computation from actual arithmetic means becomes cumbersome when or results in fractional values. Consequently, statisticians employ specific computational transformations:
Letting deviations from actual means be (x = ) and (y = ):
Formula 2: Shortcut Method (Using Assumed Means (A_x) and (A_y))
Where (dx = (X - A_x)) and (dy = (Y - A_y)):
Formula 3: Step-Deviation Method (Grouped or Scaled Series)
Where and , where (c_1) and (c_2) represent class intervals:
Crucial Mathematical Properties of Pearson's (r)
Karl Pearson's correlation coefficient possesses vital mathematical characteristics that govern its analytical behavior:
Property 1: Strictly Bounded Limits
Property 2: Invariance to Origin and Scale
Property 3: Symmetry of Relationship
4. STATISTICAL SIGNIFICANCE: PROBABLE ERROR & DETERMINATION
When a correlation coefficient (r) is calculated from a random sample of (N) pairs of observations drawn from a bivariate population, it represents a sample statistic subject to sampling fluctuations. Assessing whether the calculated (r) reflects a genuine population relationship requires statistical significance testing.
Standard Error and Probable Error of (r)
The standard error of the correlation coefficient is given by:
Probable Error of (r) ((P.E._r)):
Probable Error is a classical measure defining the limits within which the population correlation coefficient ( ho) is expected to lie with a 50% probability:
Rules for Interpreting Correlation via Probable Error
- Condition 1 (Absence of Significance): If (|r| < P.E._r), the correlation is not significant. There is no empirical evidence of any association between the variables.
- Condition 2 (Definite Significance): If , the correlation is definitely significant. The sample statistic reliably demonstrates true population correlation.
- Condition 3 (Inconclusive Zone): If , the relationship is weak or inconclusive; no definitive decision can be drawn without increasing sample size.
- Confidence Band: The limits for the true population correlation coefficient (
ho) are:
Coefficient of Determination ((r^2)) and Alienation
The correlation coefficient (r) measures the degree of linear association, but it does not represent a proportion. To evaluate the explanatory power of the relationship, statisticians use the Coefficient of Determination ((r^2)).
| Statistical Metric | Algebraic Formula | Managerial & Analytical Interpretation |
|---|---|---|
| Coefficient of Determination | Represents the proportion (or percentage) of the total variation in the dependent variable that is directly explained by its linear relationship with the independent variable. For instance, if (r = 0.80), then (r^2 = 0.64), signifying that 64% of the variance in (Y) is explained by (X), while 36% remains unexplained. | |
| Coefficient of Non-Determination | Measures the proportion of variation in the dependent variable attributable to chance, omitted factors, or random disturbances. | |
| Coefficient of Alienation | Measures the lack of association or the degree of departure from perfect linear relationship between the two variables. |
5. SPEARMAN'S RANK CORRELATION COEFFICIENT ((R) OR ( ho))
Developed in 1904 by British psychologist Charles Edward Spearman, the Rank Correlation Coefficient is a non-parametric technique applied when variables cannot be measured on an exact quantitative scale, but can be ordered or ranked based on qualitative attributes (e.g., brand loyalty, employee leadership, aesthetics, consumer taste preferences).
Let (R_1) and (R_2) be the ranks assigned to an individual on attributes (X) and (Y). Defining rank difference (D = R_1 - R_2):
Handling Tied / Repeated Ranks (Correction Factor)
When two or more items possess identical numerical values, they cannot be awarded distinct integer ranks. Under the bracket rank procedure, each tied observation is assigned the arithmetic mean of the ranks they would have occupied if they were distinct.
When ties occur, the term ((m(m^2 - 1)) / (12)) is added to for each tied group, where (m) denotes the number of items sharing that particular rank:
6. CONCURRENT DEVIATION METHOD
The Method of Concurrent Deviations is a rapid, non-parametric approach focusing purely on the direction of change from one time period or observation to the next, ignoring the magnitude of changes. It is particularly useful for short-term financial time series and volatile economic data.
- (N) = Total number of paired observations.
- (n = (N - 1)) = Total number of paired directional changes.
- (c) = Number of concurrent deviations (i.e., pairs where both signs are (+), or both signs are (-)).
- The algebraic sign is determined by the sign of ((2c - n)). If ((2c - n) > 0), the positive sign is retained both inside and outside the radical. If ((2c - n) < 0), the negative sign is retained outside, and the inner term is multiplied by (-1) to prevent imaginary numbers.
7. REGRESSION ANALYSIS: PRINCIPLES & LEAST SQUARES ARCHITECTURE
While correlation measures the degree and direction of association between two variables, it does not imply functional dependency or enable forecasting. Regression Analysis is the statistical methodology used to estimate or predict the unknown value of one variable (dependent or explained variable) from known values of another variable (independent or predictor variable).
Fundamental Distinction: Correlation vs. Regression
| Analytical Dimension | Correlation Analysis | Regression Analysis |
|---|---|---|
| 1. Objective | Measures the degree, strength, and direction of linear association. | Estimates the quantitative functional relationship to enable predictive forecasting. |
| 2. Dependency Structure | Symmetrical: Variables (X) and (Y) are treated co-equally ((r_{xy} = r_{yx})). | Asymmetrical: Differentiates strictly between cause (independent) and effect (dependent). |
| 3. Number of Equations | A single numerical coefficient ((r)) summarizes the relationship. | Yields two distinct regression lines: and . |
| 4. Dimensionality | Relative, unit-free measure (). | Absolute measures expressed in the original physical units of the dependent variable. |
| 5. Impact of Origin/Scale | Invariant to changes of both origin and scale. | Invariant to change of origin, but sensitive to change of scale. |
| 6. Error Terms | No assumption of stochastic error or disturbance term. | Explicitly minimizes the sum of squared residual errors (). |
Principle of Ordinary Least Squares (OLS)
The standard criterion for fitting a line of best fit through scattered bivariate data is the Principle of Least Squares. The method determines parameter coefficients such that the sum of the squared vertical deviations of the observed data points from the fitted line is minimized:
Normal Equations for Line of (X) on (Y) ((X = a' + b'Y)):
[
8. THE TWO BIVARIATE REGRESSION EQUATIONS & COEFFICIENTS
In simple linear bivariate analysis, there are always two distinct regression lines:
1. Regression Equation of (Y) on (X)
2. Regression Equation of (X) on (Y)
Formulas for Regression Coefficients
The regression coefficients (b_{yx}) and (b_{xy}) can be computed using several equivalent formulations depending on the available data summary:
| Method | Coefficient (b_{yx}) ((Y) on (X)) | Coefficient (b_{xy}) ((X) on (Y)) |
|---|---|---|
| From Correlation & Std Dev | (b_{yx} = r · (σ_y) / (σ_x)) | (b_{xy} = r · (σ_x) / (σ_y)) |
| From Actual Mean Deviations ((x = X - bar{X}, y = Y - bar{Y})) |
(b_{yx} = (∑ xy) / (∑ x^2)) | (b_{xy} = (∑ xy) / (∑ y^2)) |
| From Assumed Means ((dx = X - A_x, dy = Y - A_y)) |
(b_{yx} = (N ∑ dx dy - (∑ dx)(∑ dy)) / (N ∑ dx^2 - (∑ dx)^2)) | (b_{xy} = (N ∑ dx dy - (∑ dx)(∑ dy)) / (N ∑ dy^2 - (∑ dy)^2)) |
| From Raw Scores | (b_{yx} = (N ∑ XY - (∑ X)(∑ Y)) / (N ∑ X^2 - (∑ X)^2)) | (b_{xy} = (N ∑ XY - (∑ X)(∑ Y)) / (N ∑ Y^2 - (∑ Y)^2)) |
Fundamental Mathematical Properties of Regression Coefficients
The regression coefficients (b_{yx}) and (b_{xy}) exhibit critical mathematical relationships:
1. Common Sign Property
2. Geometric Mean Property
3. Complementary Magnitude
4. Arithmetic Mean Bound
5. Point of Intersection
6. Origin vs. Scale Invariance
9. ANGLE BETWEEN REGRESSION LINES & STANDARD ERROR OF ESTIMATE
The acute angle (theta) formed between the two regression lines reflects the degree of correlation between the variables. The formula for the tangent of this angle is:
[ tan theta = ≤ft| (1 - r^2) / (r) ≤ft( (σ_x · σ_y) / (σ_x^2 + σ_y^2) right) right|
Two Limiting Boundary Conditions:
• Case 1: When (Perfect Correlation):
( an heta = 0 implies heta = 0^circ). The two regression lines coincide completely into a single straight line.
• Case 2: When (r = 0) (Complete Absence of Linear Correlation):
. The two regression lines are perpendicular to each other, running parallel to the Cartesian axes ( and ).
Standard Error of Estimate ((S_{yx}) and (S_{xy}))
The Standard Error of Estimate measures the dispersion of actual observed values around the fitted regression line, providing an index of forecasting accuracy:
Standard Error of (X) on (Y) ((S_{xy})):
10. STEP-BY-STEP SOLVED NUMERICAL ILLUSTRATIONS
A manufacturing firm tracked its advertising expenditure ((X), in lakhs of ₹) and resulting sales revenue ((Y), in crores of ₹) over 8 consecutive quarterly operating periods:
Data:
Advertising Expenditure ((X)): 12, 14, 15, 18, 20, 22, 24, 27
Sales Revenue ((Y)): 20, 22, 24, 28, 32, 35, 38, 41
Required:
1. Calculate Karl Pearson's Coefficient of Correlation ((r)) using the direct mean deviation method.
2. Compute the Standard Error ((S.E._r)) and Probable Error ((P.E._r)) of the correlation coefficient.
3. Test the statistical significance of the correlation coefficient.
4. Calculate the Coefficient of Determination ((r^2)) and interpret the findings.
Solution & Step-by-Step Working Notes:
Step 1: Compute Arithmetic Means:
| Quarter | (X) | (x = X - 19) | (x^2) | (Y) | (y = Y - 30) | (y^2) | (xy) |
|---|---|---|---|---|---|---|---|
| 1 | 12 | -7 | 49 | 20 | -10 | 100 | +70 |
| 2 | 14 | -5 | 25 | 22 | -8 | 64 | +40 |
| 3 | 15 | -4 | 16 | 24 | -6 | 36 | +24 |
| 4 | 18 | -1 | 1 | 28 | -2 | 4 | +2 |
| 5 | 20 | +1 | 1 | 32 | +2 | 4 | +2 |
| 6 | 22 | +3 | 9 | 35 | +5 | 25 | +15 |
| 7 | 24 | +5 | 25 | 38 | +8 | 64 | +40 |
| 8 | 27 | +8 | 64 | 41 | +11 | 121 | +88 |
| Total | 152 | 0 | 190 | 240 | 0 | 418 | 281 |
Step 2: Calculate Karl Pearson's Correlation Coefficient ((r)):
Step 3: Standard Error and Probable Error:
Step 4: Significance Test:
Step 5: Coefficient of Determination ((r^2)):
Ten corporate acquisition targets were evaluated independently by two investment analysts (Senior Analyst (X) and Risk Officer (Y)). Their evaluations resulted in the following performance scores:
Data:
Target Company: A, B, C, D, E, F, G, H, I, J
Analyst (X) Scores: 68, 75, 90, 75, 50, 62, 80, 75, 88, 92
Risk Officer (Y) Scores: 60, 68, 85, 70, 45, 55, 80, 68, 82, 90
Required: Compute Spearman's Rank Correlation Coefficient, applying the correction factor for repeated values.
Solution & Step-by-Step Working Notes:
1. Ranking for Analyst (X):
Sorting (X) in descending order: 92 (Rank 1), 90 (Rank 2), 88 (Rank 3), 80 (Rank 4).
The value 75 appears 3 times (positions 5, 6, 7). Average rank = ((5 + 6 + 7) / (3) = 6.0). Here (m_1 = 3).
Next values: 68 (Rank 8), 62 (Rank 9), 50 (Rank 10).
2. Ranking for Risk Officer (Y):
Sorting (Y) in descending order: 90 (Rank 1), 85 (Rank 2), 82 (Rank 3), 80 (Rank 4).
Next: 70 (Rank 5).
The value 68 appears twice (positions 6 and 7). Average rank = ((6 + 7) / (2) = 6.5). Here (m_2 = 2).
Next values: 60 (Rank 8), 55 (Rank 9), 45 (Rank 10).
| Target | (X) Score | Rank (R_1) | (Y) Score | Rank (R_2) | (D = R_1 - R_2) | (D^2) |
|---|---|---|---|---|---|---|
| A | 68 | 8.0 | 60 | 8.0 | 0.0 | 0.00 |
| B | 75 | 6.0 | 68 | 6.5 | -0.5 | 0.25 |
| C | 90 | 2.0 | 85 | 2.0 | 0.0 | 0.00 |
| D | 75 | 6.0 | 70 | 5.0 | +1.0 | 1.00 |
| E | 50 | 10.0 | 45 | 10.0 | 0.0 | 0.00 |
| F | 62 | 9.0 | 55 | 9.0 | 0.0 | 0.00 |
| G | 80 | 4.0 | 80 | 4.0 | 0.0 | 0.00 |
| H | 75 | 6.0 | 68 | 6.5 | -0.5 | 0.25 |
| I | 88 | 3.0 | 82 | 3.0 | 0.0 | 0.00 |
| J | 92 | 1.0 | 90 | 1.0 | 0.0 | 0.00 |
| Total | - | - | - | - | (∑ D = 0) | (∑ D^2 = 1.50) |
3. Compute Correction Factors:
For (X) series: Value 75 repeated (m_1 = 3) times:
[
CF_1 = (m_1(m_1^2 - 1)) / (12) = (3(3^2 - 1)) / (12) = (3 × 8) / (12) = 2.0
4. Calculate Corrected Spearman's Rank Correlation:
An industrial plant operations director examined the relationship between factory capacity utilization rate ((X), in %) and total factory monthly maintenance expenditure ((Y), in thousands of ₹) across 10 operating months:
Data:
Capacity Utilization ((X)): 60, 62, 65, 70, 72, 75, 78, 80, 82, 86
Maintenance Expenditure ((Y)): 25, 28, 30, 36, 38, 40, 42, 45, 48, 52
Required:
1. Calculate the arithmetic means and , and standard deviations and .
2. Compute the two regression coefficients (b_{yx}) and (b_{xy}), and deduce the correlation coefficient (r).
3. Construct the regression equation of (Y) on (X) and forecast maintenance expenditure for a projected capacity utilization rate of 85%.
4. Construct the regression equation of (X) on (Y) and estimate capacity utilization if maintenance budget is restricted to ₹ 50,000 ((Y = 50)).
5. Compute the Standard Error of Estimate (S_{yx}).
Solution & Step-by-Step Working Notes:
Let Assumed Mean for (X) be (A_x = 72), and for (Y) be (A_y = 38).
Then (dx = X - 72) and (dy = Y - 38).
| Month | (X) | (dx = X - 72) | (dx^2) | (Y) | (dy = Y - 38) | (dy^2) | |
|---|---|---|---|---|---|---|---|
| 1 | 60 | -12 | 144 | 25 | -13 | 169 | +156 |
| 2 | 62 | -10 | 100 | 28 | -10 | 100 | +100 |
| 3 | 65 | -7 | 49 | 30 | -8 | 64 | +56 |
| 4 | 70 | -2 | 4 | 36 | -2 | 4 | +4 |
| 5 | 72 | 0 | 0 | 38 | 0 | 0 | 0 |
| 6 | 75 | +3 | 9 | 40 | +2 | 4 | +6 |
| 7 | 78 | +6 | 36 | 42 | +4 | 16 | +24 |
| 8 | 80 | +8 | 64 | 45 | +7 | 49 | +56 |
| 9 | 82 | +10 | 100 | 48 | +10 | 100 | +100 |
| 10 | 86 | +14 | 196 | 52 | +14 | 196 | +196 |
| Total | 730 | +10 | 702 | 384 | +4 | 702 | 698 |
Step 1: Compute Means and Standard Deviations:
Step 2: Compute Regression Coefficients (b_{yx}) and (b_{xy}):
Step 3: Regression Equation of (Y) on (X) (Predicting Maintenance Cost):
Step 5: Standard Error of Estimate ((S_{yx})):
11. COMPREHENSIVE SELF-ASSESSMENT MODULE
Select the most appropriate option. Complete solution key with mathematical rationale follows.
1. If Karl Pearson's correlation coefficient (r = +0.80), the percentage of variation in the dependent variable left unexplained by the independent variable is:
(a) 80% (b) 20% (c) 64% (d) 36%
2. If the two regression coefficients are (b_{yx} = -1.6) and (b_{xy} = -0.4), the value of the correlation coefficient (r) is:
(a) +0.80 (b) -0.80 (c) -0.64 (d) +0.64
3. Which of the following statements regarding regression coefficients is mathematically FALSE?
(a) Both regression coefficients must share the same algebraic sign.
(b) If (b_{yx} > 1), then (b_{xy}) must be less than 1.
(c) Regression coefficients are independent of change of scale.
(d) The arithmetic mean of (b_{yx}) and (b_{xy}) is greater than or equal to (r).
4. If the acute angle between the two regression lines is ( heta = 0^circ), the coefficient of correlation (r) is equal to:
(a) 0 (b) (c) (d)
5. The two lines of regression intersect at the point defined by:
(a) ((0, 0)) (b) () (c) () (d) ((b_{yx}, b_{xy}))
7. If in a tied rank correlation problem, a score repeats 4 times, the correction factor ((m(m^2 - 1)) / (12)) added to is:
(a) 4 (b) 5 (c) 16 (d) 2.5
8. If (u = 2X - 5) and (v = -3Y + 4), and the correlation between (X) and (Y) is (r_{xy} = 0.75), the correlation between (u) and (v) is:
(a) +0.75 (b) -0.75 (c) +1.50 (d) -0.50
9. The Standard Error of Estimate (S_{yx}) becomes equal to the standard deviation when:
(a) (r = 1.0) (b) (r = -1.0) (c) (r = 0) (d) (b_{yx} = 1.0)
10. In the Method of Concurrent Deviations, if (N = 11) and the number of concurrent deviations (c = 2), then ((2c - n)) evaluates to:
(a) -6 (b) +6 (c) -7 (d) +4
1. (d) 36%: (r^2 = (0.80)^2 = 0.64) (explained variance). Unexplained variance is (1 - r^2 = 1 - 0.64 = 0.36 = 36%).
2. (b) -0.80: (r = -√((-1.6) × (-0.4)) = -√(0.64) = -0.80). The sign of (r) must match the signs of the coefficients.
3. (c): Regression coefficients are dependent on change of scale (); only correlation is independent of scale.
4. (c) : When , ( an heta = 0 implies heta = 0^circ), meaning both lines coincide.
5. (c) (): The lines (Y - Ȳ = b_{yx}) and (X - X̄ = b_{xy}) simultaneously satisfy .
6. (c) : By Horace Secrist and classical statistical doctrine, significance is assured when (r) exceeds 6 times its probable error.
7. (b) 5: ((4(4^2 - 1)) / (12) = (4(15)) / (12) = (60) / (12) = 5.0).
8. (b) -0.75: The scale transformation on (Y) involves multiplying by (-3) (a negative scalar), while (X) is multiplied by (+2). An odd number of negative scale changes flips the sign of (r).
9. (c) (r = 0): . When no correlation exists, regression line prediction error equals total standard deviation.
10. (a) -6: (n = N - 1 = 11 - 1 = 10). Thus (2c - n = 2(2) - 10 = 4 - 10 = -6).
- Why are there two regression lines in bivariate analysis? Under what operational circumstances will the two regression lines coincide, and under what circumstances will they be mutually perpendicular? Explain with mathematical proofs.
- Evaluate the properties of Regression Coefficients. Prove mathematically that the arithmetic mean of the two regression coefficients is always greater than or equal to the geometric mean (correlation coefficient).
- Contrast Karl Pearson's Correlation with Spearman's Rank Correlation. Delineate specific commercial and behavioral circumstances where Spearman's rank method is statistically superior to product-moment correlation.
- Define the Coefficient of Determination ((r^2)). How does it distinguish genuine explanatory power from simple co-variation? Explain its role in multiple regression analysis.
- Explain the role of Probable Error in Statistical Inference. State the precise conditions under which Probable Error can be legitimately employed to assess the reliability of a sample correlation coefficient.
Problem 1: Bivariate Sales & Production Regression Analysis
A consumer goods corporation records the following summary statistics from 100 operating branch stores across South India:
• Mean Annual Branch Promotion Outlay ((X)) = ₹ 50,000, Standard Deviation () = ₹ 10,000
• Mean Annual Branch Sales Revenue ((Y)) = ₹ 400,000, Standard Deviation () = ₹ 50,000
• Karl Pearson's Correlation Coefficient between Promotion and Sales ((r)) = +0.80
Required:
1. Compute the two regression coefficients (b_{yx}) and (b_{xy}).
2. Construct the regression equation of Sales ((Y)) on Promotion Outlay ((X)).
3. Forecast the expected sales revenue of a branch that spends ₹ 70,000 on promotion.
4. Construct the regression equation of Promotion Outlay ((X)) on Sales ((Y)) and determine the promotion expenditure recommended to achieve a target branch sales revenue of ₹ 550,000.
5. Calculate the Standard Error of Estimate (S_{yx}) and construct a 95% confidence interval () for the sales forecast in (3).
Problem 2: Spearman's Rank Correlation with Severe Ties
Two market research firms ranked 8 competing mobile smartphone brands based on customer satisfaction metrics. The assigned scores were as follows:
Brand: A, B, C, D, E, F, G, H
Agency 1: 85, 70, 70, 90, 60, 70, 85, 55
Agency 2: 78, 65, 72, 85, 50, 65, 78, 48
Required:
1. Assign average fractional ranks to both agencies' scores.
2. Identify the tied groups and compute the individual and aggregate correction factors.
3. Calculate the corrected Spearman's Rank Correlation Coefficient ((R)) and comment on the degree of agreement between the two rating agencies.