Correlation analysis in R Programming - Step-by-Step Execution

Concept Flow - Correlation analysis

Start with two numeric vectors

↓

Calculate means of each vector

↓

Calculate deviations from means

↓

Multiply deviations element-wise

↓

Sum multiplied deviations

↓

Calculate standard deviations

↓

Divide sum by product of std devs and n-1

↓

End

Correlation analysis calculates a number showing how two sets of numbers move together, step-by-step computing means, deviations, and combining them into a coefficient.

Execution Sample

R Programming

x <- c(1, 2, 3, 4, 5)
y <- c(2, 4, 6, 8, 10)
correlation <- cor(x, y)
print(correlation)

This code calculates the correlation coefficient between two numeric vectors x and y and prints the result.

Execution Table

Step	Action	Value/Calculation	Result
1	Input vectors	x = (1,2,3,4,5), y = (2,4,6,8,10)	Vectors ready
2	Calculate means	mean(x) = 3, mean(y) = 6	Means computed
3	Calculate deviations	x_dev = (-2,-1,0,1,2), y_dev = (-4,-2,0,2,4)	Deviations computed
4	Multiply deviations element-wise	x_dev * y_dev = (8,2,0,2,8)	Products computed
5	Sum multiplied deviations	sum = 20	Sum = 20
6	Calculate standard deviations	sd(x) = 1.58, sd(y) = 3.16	Std devs computed
7	Calculate correlation	r = 20 / ((5-1)*1.58*3.16) = 1	Correlation coefficient = 1
8	Print result	print(1)	Output: 1
9	End	Correlation calculated successfully	Stop

💡 All steps completed; correlation coefficient calculated and printed.

Variable Tracker

Variable	Start	After Step 2	After Step 3	After Step 4	After Step 5	After Step 6	After Step 7	Final
x	(1,2,3,4,5)	(1,2,3,4,5)	(-2,-1,0,1,2)	(-2,-1,0,1,2)	(-2,-1,0,1,2)	(1,2,3,4,5)	(1,2,3,4,5)	(1,2,3,4,5)
y	(2,4,6,8,10)	(2,4,6,8,10)	(-4,-2,0,2,4)	(-4,-2,0,2,4)	(-4,-2,0,2,4)	(2,4,6,8,10)	(2,4,6,8,10)	(2,4,6,8,10)
mean_x	NA	3	3	3	3	3	3	3
mean_y	NA	6	6	6	6	6	6	6
x_dev	NA	NA	(-2,-1,0,1,2)	(-2,-1,0,1,2)	(-2,-1,0,1,2)	NA	NA	NA
y_dev	NA	NA	(-4,-2,0,2,4)	(-4,-2,0,2,4)	(-4,-2,0,2,4)	NA	NA	NA
prod_dev	NA	NA	NA	(8,2,0,2,8)	(8,2,0,2,8)	NA	NA	NA
sum_prod_dev	NA	NA	NA	NA	20	NA	NA	NA
sd_x	NA	NA	NA	NA	NA	1.58	1.58	1.58
sd_y	NA	NA	NA	NA	NA	3.16	3.16	3.16
correlation	NA	NA	NA	NA	NA	NA	1	1

Key Moments - 3 Insights

Why do we subtract the mean from each value before multiplying?

Subtracting the mean centers the data around zero, so we measure how values move relative to their average, as shown in step 3 of the execution_table.

Why do we divide by (n-1) times the product of standard deviations?

Dividing by (n-1) and the product of standard deviations normalizes the sum of products to get a value between -1 and 1, explained in step 7 of the execution_table.

What does a correlation of 1 mean in this example?

A correlation of 1 means perfect positive linear relationship; as x increases, y increases proportionally, as seen in the final result in step 8.

Visual Quiz - 3 Questions

Look at the execution_table at step 3. What are the deviations of x from its mean?

A (2, 1, 0, -1, -2)

B (1, 2, 3, 4, 5)

C (-2, -1, 0, 1, 2)

D (0, 0, 0, 0, 0)

💡 Hint

Check the 'Value/Calculation' column at step 3 in execution_table.

At which step does the correlation coefficient get calculated?

A Step 7

B Step 3

C Step 5

D Step 2

💡 Hint

Look for the step where division by standard deviations and n-1 happens in execution_table.

If the vectors x and y were identical, what would the correlation coefficient be according to variable_tracker?

A -1

B 1

C 0

D Undefined

💡 Hint

Correlation of identical vectors is shown in the final correlation value in variable_tracker.

Concept Snapshot

Correlation analysis in R:
- Use cor(x, y) to find correlation coefficient r
- r ranges from -1 (perfect negative) to 1 (perfect positive)
- Steps: compute means, deviations, multiply, sum, divide by std dev product
- Measures linear relationship strength and direction
- Perfect correlation means one variable predicts the other exactly

Full Transcript

Correlation analysis calculates how two numeric vectors relate linearly. We start with two vectors, find their means, then find how each value deviates from its mean. We multiply these deviations element-wise and sum them. Then we calculate the standard deviations of each vector. Finally, we divide the sum of multiplied deviations by the product of standard deviations and n-1 to get the correlation coefficient. This coefficient tells us if the vectors move together positively, negatively, or not at all. In the example, vectors x and y have a perfect positive correlation of 1.

Practice

1. What does the cor() function in R calculate?

easy

A. The correlation coefficient between two numeric vectors

B. The sum of two numeric vectors

C. The difference between two numeric vectors

D. The product of two numeric vectors

Solution

1. Step 1: Understand the purpose of cor()

   The cor() function is designed to calculate how two numeric vectors move together, which is called correlation.

2. Step 2: Identify the correct description

   Correlation coefficient measures the strength and direction of a linear relationship between two numeric vectors.

3. Final Answer:

   The correlation coefficient between two numeric vectors -> Option A

4. Quick Check:

   cor() calculates correlation [OK]

Hint: Remember: cor() finds how two numbers move together [OK]

Common Mistakes:

• Confusing correlation with sum or difference
• Thinking cor() calculates causation
• Using cor() on non-numeric data

2. Which of the following is the correct syntax to calculate correlation between vectors x and y in R?

easy

A. correlation_coefficient(x, y)

B. cor(x, y)

C. correlate(x, y)

D. correlation(x, y)

Solution

1. Step 1: Recall the R function name

   The correct function to calculate correlation in R is cor().

2. Step 2: Check the syntax

   The function takes two numeric vectors as arguments: cor(x, y).

3. Final Answer:

   cor(x, y) -> Option B

4. Quick Check:

   Use cor() with two vectors [OK]

Hint: Use exactly cor(x, y) to get correlation [OK]

Common Mistakes:

• Using incorrect function names like correlation()
• Adding extra arguments unnecessarily
• Using wrong argument order

3. What is the output of the following R code?

x <- c(1, 2, 3, 4, 5)
y <- c(2, 4, 6, 8, 10)
cor(x, y)

medium

A. 1

B. -1

C. 0

D. NA

Solution

1. Step 1: Analyze the vectors

   Vector y is exactly twice x for each element, so they have a perfect positive linear relationship.

2. Step 2: Calculate correlation

   Perfect positive linear relationship means correlation coefficient is 1.

3. Final Answer:

   1 -> Option A

4. Quick Check:

   Perfect positive correlation = 1 [OK]

Hint: If one vector is a perfect multiple of another, cor = 1 [OK]

Common Mistakes:

• Confusing correlation with sum or mean
• Expecting 0 when vectors are related
• Misreading negative correlation

4. Identify the error in this R code for correlation:

x <- c(1, 2, 3)
y <- c(4, 5)
cor(x, y)

medium

A. cor() function is misspelled

B. Vectors contain non-numeric data

C. Vectors have different lengths causing an error

D. No error, code runs fine

Solution

1. Step 1: Check vector lengths

   Vector x has length 3, vector y has length 2, so lengths differ.

2. Step 2: Understand cor() requirements

   cor() requires vectors of equal length to compute correlation.

3. Final Answer:

   Vectors have different lengths causing an error -> Option C

4. Quick Check:

   cor() needs equal length vectors [OK]

Hint: Vectors must be same length for cor() [OK]

Common Mistakes:

• Ignoring vector length mismatch
• Assuming cor() auto-adjusts lengths
• Thinking cor() works with non-numeric data

5. You have two numeric vectors a and b with some missing values (NA). Which is the correct way to calculate correlation ignoring missing values in R?

hard

A. cor(a, b, removeNA = TRUE)

B. cor(a, b, na.rm = TRUE)

C. cor(a, b, method = "ignoreNA")

D. cor(a, b, use = "complete.obs")

Solution

1. Step 1: Understand handling of missing values in cor()

   The cor() function uses the use argument to specify how to handle missing values.

2. Step 2: Identify correct argument value

   Setting use = "complete.obs" tells R to use only pairs without any NA values.

3. Final Answer:

   cor(a, b, use = "complete.obs") -> Option D

4. Quick Check:

   Use use = "complete.obs" to ignore NA [OK]

Hint: Use use="complete.obs" to skip NA pairs [OK]

Common Mistakes:

• Using na.rm which is not an argument for cor()
• Using incorrect argument names like removeNA
• Assuming method argument handles NA