Principal Component Analysis (PCA) in MATLAB

Introduction

Principal Component Analysis (PCA) helps to simplify data by finding the main directions where the data changes the most. It reduces many numbers into fewer important ones.

You have many measurements and want to find the most important patterns.

You want to reduce data size to make it easier to analyze or visualize.

You want to remove noise or less important details from your data.

You want to prepare data before using machine learning.

You want to understand relationships between variables in your data.

Syntax

MATLAB

coeff = pca(X)
[coeff, score, latent] = pca(X)

% X is a data matrix where rows are observations and columns are variables.

pca function returns principal components in coeff.

score gives the data projected onto principal components.

Examples

Calculate principal components for a small 4x2 data matrix.

MATLAB

X = [2.5 2.4; 0.5 0.7; 2.2 2.9; 1.9 2.2];
coeff = pca(X);

Get principal components and transformed data.

MATLAB

[coeff, score] = pca(X);

Also get variance explained by each principal component.

MATLAB

[coeff, score, latent] = pca(X);

Sample Program

This program runs PCA on a 10x2 data matrix. It shows the main directions (coeff), the data in the new system (score), and how much each direction explains the data (latent).

MATLAB

X = [2.5 2.4; 0.5 0.7; 2.2 2.9; 1.9 2.2; 3.1 3.0; 2.3 2.7; 2 1.6; 1 1.1; 1.5 1.6; 1.1 0.9];
[coeff, score, latent] = pca(X);

fprintf('Principal Components (coeff):\n');
disp(coeff);

fprintf('Data projected on principal components (score):\n');
disp(score);

fprintf('Variance explained by each component (latent):\n');
disp(latent);

Important Notes

Make sure your data matrix X has rows as samples and columns as features.

PCA centers the data by default before calculating components.

Use latent to understand how much each component matters.

Summary

PCA finds main directions to simplify data.

Use pca function in MATLAB to get components and transformed data.

Check variance latent to see importance of each component.

Practice

1.

What is the main purpose of Principal Component Analysis (PCA) in MATLAB?

easy

A. To reduce the number of variables by finding main directions of data

B. To sort data in ascending order

C. To calculate the sum of all data points

D. To generate random numbers for data simulation

Solution

1. Step 1: Understand PCA goal

   PCA finds new axes (principal components) that capture most data variation.

2. Step 2: Identify main purpose

   It reduces data dimensions by focusing on important directions, simplifying analysis.

3. Final Answer:

   To reduce the number of variables by finding main directions of data -> Option A

4. Quick Check:

   PCA = Dimensionality reduction [OK]

Hint: PCA simplifies data by finding main directions [OK]

Common Mistakes:

• Thinking PCA sorts data
• Confusing PCA with summing data
• Assuming PCA generates random data

2.

Which MATLAB command correctly performs PCA on a data matrix X?

easy

A. score = pca(X, 'score');

B. pca = X.coeff();

C. [coeff, score, latent] = pca(X);

D. coeff = pca(X, 'coeff');

Solution

1. Step 1: Recall MATLAB PCA syntax

   The function pca returns coefficients, scores, and variances as outputs.

2. Step 2: Match correct syntax

   Only [coeff, score, latent] = pca(X); correctly calls PCA and stores outputs.

3. Final Answer:

   [coeff, score, latent] = pca(X); -> Option C

4. Quick Check:

   Correct PCA syntax = [coeff, score, latent] = pca(X); [OK]

Hint: Use brackets to capture multiple PCA outputs [OK]

Common Mistakes:

• Using dot notation on matrix
• Passing output names as strings
• Assigning only one output without brackets

3.

Given the MATLAB code below, what is the size of score if X is a 100x5 matrix?

 [coeff, score, latent] = pca(X); 

medium

A. 5x5

B. 100x5

C. 100x100

D. 5x100

Solution

1. Step 1: Understand PCA output sizes

   score contains transformed data points in principal component space, same rows as X.

2. Step 2: Match dimensions

   Since X is 100x5, score is 100x5 (100 samples, 5 components).

3. Final Answer:

   100x5 -> Option B

4. Quick Check:

   score size = samples x components = 100x5 [OK]

Hint: Score matches rows of X and number of components [OK]

Common Mistakes:

• Confusing score size with coeff size
• Assuming score is square matrix
• Swapping rows and columns

4.

Identify the error in this MATLAB PCA code snippet:

 X = rand(50,3);
 [coeff, score, latent] = pca(X);
 disp(score(1,4)); 

medium

A. Indexing score at column 4 causes an error

B. Missing semicolon after rand function

C. Incorrect function name, should be pca_analysis

D. Matrix X size is invalid for PCA

Solution

1. Step 1: Check size of PCA outputs

   Input X has 3 columns, so score has 3 columns (components).

2. Step 2: Identify indexing error

   Accessing score(1,4) tries to get 4th column which does not exist, causing an error.

3. Final Answer:

   Indexing score at column 4 causes an error -> Option A

4. Quick Check:

   score columns = 3, index 4 invalid [OK]

Hint: Check matrix size before indexing columns [OK]

Common Mistakes:

• Assuming more components than input columns
• Ignoring MATLAB indexing errors
• Thinking semicolon omission causes error here

5.

You want to reduce a dataset X (200x10) to 2 principal components using MATLAB PCA. Which code correctly achieves this?

 % Options:
 A) [coeff, score] = pca(X, 'NumComponents', 2);
 B) [coeff, score] = pca(X, 'NumComponents', 2); score = score(:,3:4);
 C) [coeff, score] = pca(X); score = score(:,11:12);
 D) [coeff, score] = pca(X, 'Components', 2); 

hard

A. [coeff, score] = pca(X, 'Components', 2);

B. [coeff, score] = pca(X, 'NumComponents', 2); score = score(:,3:4);

C. [coeff, score] = pca(X); score = score(:,11:12);

D. [coeff, score] = pca(X, 'NumComponents', 2);

Solution

1. Step 1: Understand PCA dimension reduction parameter

   MATLAB's pca accepts 'NumComponents' to limit output components.

2. Step 2: Identify correct usage

   [coeff, score] = pca(X, 'NumComponents', 2); uses correct parameter 'NumComponents', 2 to get 2 components directly.

3. Step 3: Check other options

   B: [coeff, score] = pca(X, 'NumComponents', 2); score = score(:,3:4); tries to slice non-existent columns 3-4; C: [coeff, score] = pca(X); score = score(:,11:12); tries to slice non-existent columns 11-12; D: [coeff, score] = pca(X, 'Components', 2); uses invalid parameter.

4. Final Answer:

   [coeff, score] = pca(X, 'NumComponents', 2); -> Option D

5. Quick Check:

   Use 'NumComponents' to limit PCA output [OK]

Hint: Use 'NumComponents' to get desired PCA components [OK]

Common Mistakes:

• Using wrong parameter name 'Components'
• Using invalid slice indices on PCA score
• Assuming PCA auto-limits components without parameter