#### Correlation

KEY:

ds = dataset you are currently using.

Var1, Var2, etc. = Variables in the dataset

XYXY = dummy name for a variable, matrix, or data frame into which you are moving information.

## Topics

## *Correlation between two variables in a data set

*Correlations between two variables in a data set...with p-value

*Correlations among all variables in a data set

*Correlations among all variables in a data set...with p-values

*Partial Correlation

*Semipartial Correlation

*Bayesian correlation

## Correlation between two variables in a data set

Note that IV = Variable that will go on X axis. DV = variable that will go on Y axis.

_____________________________

cor(ds$Var1, ds$Var2)

_____________________________

## Correlations between two variables in a data set...with p-value

Note that IV = Variable that will go on X axis. DV = variable that will go on Y axis.

_____________________________

cor.test(ds$Var1, ds$Var2, method="pearson")

_____________________________

## Correlations among all variables in a data set

Note [-1] removes the first column (i.e., subject number) from the analysis.

_____________________________

cor(ds[-1], method="pearson")

_____________________________

## Correlations among all variables in a data set...with p-values

You will need to download the "Hmisc" package in order to obtain the "rcorr" function.

Note [-c(1)] removes the first column (i.e., subject number) from the analysis.

_____________________________

library(Hmisc)

rcorr(as.matrix(ds[-c(1)]), type= "pearson")

_____________________________

## Partial Correlation

You will need to download the "ppcor" package in order to obtain the "pcor.test" function.

Var1 and Var2 is the correlation of interest. VarZ is the factor that is partialled out.

_____________________________

library(ppcor)

pcor.test(x=ds$Var1, y=ds$Var2, z=ds$VarZ, method="pearson")

_____________________________

## Semipartial Correlation

You will need to download the "ppcor" package in order to obtain the "spcor.test" function.

Var1 and Var2 is the correlation of interest. VarZ is the factor that is partialled out of Var2 only.

_____________________________

library(ppcor)

spcor.test(x=ds$Var1, y=ds$Var2, z=ds$VarZ, method="pearson")

_____________________________

## Bayesian correlation

Bayesian correlations can be helpful when you want to know how much trust you can place in an observed correlation value. This can be particularly helpful in the case of non-significant relationships.

You will need to download the "BayesFactor" (to create the model) and "bayestestR" (to expand the descriptive capabilities) packages.

These commands are intended to get you up and running. For extensive discussion of these packages, please visit https://richarddmorey.github.io/BayesFactor/ and https://easystats.github.io/bayestestR/

_____________________________

library(BayesFactor)

XYXY <- correlationBF(DS$Var1, DS$Var2)

library(bayestestR)

describe_posterior(XYXY)

bayesfactor(XYXY)

_____________________________