R: Set Correlation and Multiple Regression from matrix or raw...

setCor {psych} | R Documentation

## Set Correlation and Multiple Regression from matrix or raw input

### Description

Finds Cohen's Set Correlation between a predictor set of variables (x) and a criterion set (y). Also finds multiple correlations between x variables and each of the y variables. Will work with either raw data or a correlation matrix. A set of covariates (z) can be partialled from the x and y sets. Regression diagrams are automatically included.

### Usage

    setCor(y,x,data,z=NULL,n.obs=NULL,use="pairwise",std=TRUE,square=FALSE,
           main="Regression Models",plot=TRUE)
    setCor.diagram(sc,main="Regression model",digits=2,show=FALSE,...)
    set.cor(y,x,data,z=NULL,n.obs=NULL,use="pairwise",std=TRUE,square=FALSE,
           main="Regression Models",plot=TRUE)    #an alias to setCor
    mat.regress(y, x,data, z=NULL,n.obs=NULL,use="pairwise",square=FALSE)
    matReg(x,y,C,n.obs=0)

### Arguments

`y` | either the column numbers of the y set (e.g., c(2,4,6) or the column names of the y set (e.g., c("Flags","Addition")
`x` | either the column numbers of the x set (e.g., c(1,3,5) or the column names of the x set (e.g. c("Cubes","PaperFormBoard")
`data` | a matrix or data.frame of correlations or, if not square, of raw data
`C` | A variance/covariance matrix, or a correlation matrix
`z` | the column names or numbers of the set of covariates
`n.obs` | If specified, then confidence intervals, etc. are calculated, not needed if raw data are given
`use` | find the correlations "pairwise" (default) or just use "complete" cases (to match the lm function)
`std` | Report standardized betas (based upon the correlations) or raw (based upon covariances)
`main` | The title for setCor.diagram
`square` | if FALSE, then square matrices are treated as correlation matrices not as data matrices. In the rare case that one has as many cases as variables, then set square=TRUE.
`sc` | The output of setCor may be used for drawing diagrams
`digits` | How many digits should be displayed in the setCor.diagram?
`show` | Show the unweighted matrix correlation between the x and y sets?
`plot` | By default, setCor makes a plot of the results, set to FALSE to suppress the plot
`...` | Additional graphical parameters for setCor.diagram

### Details

Although it is more common to calculate multiple regression and canonical correlations from the raw data, it is, of course, possible to do so from a matrix of correlations or covariances. In this case, the input to the function is a square covariance or correlation matrix, as well as the column numbers (or names) of the x (predictor), y (criterion) variables, and if desired z (covariates). The function will find the correlations if given raw data.

The output is a set of multiple correlations, one for each dependent variable in the y set, as well as the set of canonical correlations.

An additional output is the R2 found using Cohen's set correlation (Cohen, 1982). This is a measure of how much variance and the x and y set share.

Cohen (1982) introduced the set correlation, a multivariate generalization of the multiple correlation to measure the overall relationship between two sets of variables. It is an application of canoncial correlation (Hotelling, 1936) and 1 - ∏(1-ρ_i^2) where ρ_i^2 is the squared canonical correlation. Set correlation is the amount of shared variance (R2) between two sets of variables. With the addition of a third, covariate set, set correlation will find multivariate R2, as well as partial and semi partial R2. (The semi and bipartial options are not yet implemented.) Details on set correlation may be found in Cohen (1982), Cohen (1988) and Cohen, Cohen, Aiken and West (2003).

R2 between two sets is just

R2 = 1- |R| /(|Ry| * |Rx|)

where R is the complete correlation matrix of the x and y variables and Rx and Ry are the two sets involved.

Unfortunately, the R2 is sensitive to one of the canonical correlations being very high. An alternative, T2, is the proportion of additive variance and is the average of the squared canonicals. (Cohen et al., 2003), see also Cramer and Nicewander (1979). This average, because it includes some very small canonical correlations, will tend to be too small. Cohen et al. admonition is appropriate: "In the final analysis, however, analysts must be guided by their substantive and methodological conceptions of the problem at hand in their choice of a measure of association." ( p613).

Yet another measure of the association between two sets is just the simple, unweighted correlation between the two sets. That is,

Ruw=1Rxy1' / (sqrt(1Ryy1'* 1Rxx1'))

where Rxy is the matrix of correlations between the two sets. This is just the simple (unweighted) sums of the correlations in each matrix. This technique exemplifies the robust beauty of linear models and is particularly appropriate in the case of one dimension in both x and y, and will be a drastic underestimate in the case of items where the betas differ in sign.

When finding the unweighted correlations, as is done in `alpha`, items are flipped so that they all are positively signed.

A typical use in the SAPA project is to form item composites by clustering or factoring (see `fa`,`ICLUST`, `principal`), extract the clusters from these results (`factor2cluster`), and then form the composite correlation matrix using `cluster.cor`. The variables in this reduced matrix may then be used in multiple R procedures using `set.cor`.

Although the overall matrix can have missing correlations, the correlations in the subset of the matrix used for prediction must exist.

If the number of observations is entered, then the conventional confidence intervals, statistical significance, and shrinkage estimates are reported.

If the input is rectangular (not square), correlations or covariances are found from the data.

The print function reports t and p values for the beta weights, the summary function just reports the beta weights.

The Variance Inflation Factor is reported but should be taken with the normal cautions of interpretation discussed by Guide and Ketokivm. That is to say, VIF > 10 is not a magic cuttoff to define colinearity. It is merely 1/(1-smc(R(x)).

`matReg` is primarily a helper function for `mediate` but is a general multiple regression function given a covariance matrix and the specified x, and y variables. Its output includes betas, se, t, p and R2. It does not work on data matrices, nor does it take formula input.

### Value

`beta` | the beta weights for each variable in X for each variable in Y
`R` | The multiple R for each equation (the amount of change a unit in the predictor set leads to in the criterion set).
`R2` | The multiple R2 (% variance acounted for) for each equation
`VIF` | The Variance Inflation Factor which is just 1/(1-smc(x))
`se` | Standard errors of beta weights (if n.obs is specified)
`t` | t value of beta weights (if n.obs is specified)
`Probability` | Probability of beta = 0 (if n.obs is specified)
`shrunkenR2` | Estimated shrunken R2 (if n.obs is specified)
`setR2` | The multiple R2 of the set correlation between the x and y sets

itemresidualThe residual correlation matrix of Y with x and z removed

`ruw` | The unit weighted multiple correlation
`Ruw` | The unit weighted set correlation

### Note

As of April 30, 2011, the order of x and y was swapped in the call to be consistent with the general y ~ x syntax of the lm and aov functions. In addition, the primary name of the function was switched to setCor from mat.regress to reflect the estimation of the set correlation.

The denominator degrees of freedom for the set correlation does not match that reported by Cohen et al., 2003 in the example on page 621 but does match the formula on page 615, except for the typo in the estimation of F (see Cohen 1982). The difference seems to be that they are adding in a correction factor of df 2 = df2 + df1.

### Author(s)

William Revelle

Maintainer: William Revelle <revelle@northwestern.edu>

### References

J. Cohen (1982) Set correlation as a general multivariate data-analytic method. Multivariate Behavioral Research, 17(3):301-341.

J. Cohen, P. Cohen, S.G. West, and L.S. Aiken. (2003) Applied multiple regression/correlation analysis for the behavioral sciences. L. Erlbaum Associates, Mahwah, N.J., 3rd ed edition.

H. Hotelling. (1936) Relations between two sets of variates. Biometrika 28(3/4):321-377.

E.Cramer and W. A. Nicewander (1979) Some symmetric, invariant measures of multivariate association. Psychometrika, 44:43-54.

V. D. R. Guide Jr. and M. Ketokivim (2015) Notes from the Editors: Redefining some methodological criteria for the journal. Journal of Operations Management. 37. v-viii.

### See Also

`cluster.cor`, `factor2cluster`,`principal`,`ICLUST`, `link{cancor}` and cca in the yacca package.

### Examples

    #the Kelly data from Hoteling
    kelly <- structure(list(speed = c(1, 0.4248, 0.042, 0.0215, 0.0573), power = c(0.4248,
    1, 0.1487, 0.2489, 0.2843), words = c(0.042, 0.1487, 1, 0.6693,
    0.4662), symbols = c(0.0215, 0.2489, 0.6693, 1, 0.6915), meaningless = c(0.0573,
    0.2843, 0.4662, 0.6915, 1)), .Names = c("speed", "power", "words",
    "symbols", "meaningless"), class = "data.frame", row.names = c("speed",
    "power", "words", "symbols", "meaningless"))

    kelly

    setCor(1:2,3:5,kelly)
    #Hotelling reports canonical correlations of .3073 and .0583  or squared correlations of
    # 0.09443329 and 0.00339889 vs. our values of  0.0946 0.0035,

    setCor(y=c(7:9),x=c(1:6),data=Thurstone,n.obs=213)
    #now try partialling out some variables
    set.cor(y=c(7:9),x=c(1:3),z=c(4:6),data=Thurstone) #compare with the previous
    #compare complete print out with summary printing
    sc <- setCor(x=c("gender","education"),y=c("SATV","SATQ"),data=sat.act) # regression from raw data
    sc
    summary(sc)

[Package psych version 1.7.3.21 Index]