# Does it make sense to standardize Principal Components after performing Principal Component Analysis?

I am attempting to emulate the following paper by Messer et al., using year 2000 decennial Census data in R to create an index known as the Neighborhood Deprivation Index(NDI): https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3261293/#CR73

The Component Extraction and Index Construction section of the paper outlines the steps it took to construct the index using Principal Components Analysis.

I am aware that the variables used in PCA are often standardized before performing the method, as PCA is sensitive to the variances of the initial variables. If there are large differences between the ranges of the initial variables, those variables with larger ranges will dominate over those with small ranges, leading to biased results. However, the final step in the process outlined in the paper states:

The deprivation index was then standardized to have a mean of 0 and standard deviation (SD) of 1 by dividing the index by the square of the eigenvalue.

I read this statement as the resulting principal components being standardized and not that the data used in the creation of the principal components were standardized.

There are a few problems I've faced in trying to emulate their step:

1. The paper is a bit vague
2. There is no attached code where this paper is available
3. I couldn't find anything on post PC construction standardization

So I have a few questions.

1. Am I right in understanding the last step as standardizing the principal components?
2. Does it make sense to standardize the principal components if you standardized the data before performing PCA?
3. From my understanding, there should be an eigenvalue produced for each variable used in the construction of the PCA. If the last step is stating I should divide the index (the PCs?) by the eigenvalue^2, how am I supposed to do that if there are 8 eigenvalues for this specific scenario?
• Welcome to Cross Validated! Aren't the eigenvalues the variance, so we would divide by their square roots to standardize? Did you mean to type that it is the square root?
– Dave
Commented Feb 9, 2022 at 15:32
• Thanks Dave. No I directly took the quote from the paper. It could be a typo in the paper though. Commented Feb 9, 2022 at 17:04

• Thanks for your input Dave. I am also starting to think this is a typo. I know the principal() function from the psych package can standardize the output using the square root of the eigenvalues, so that may just be a uncorrected mistake in the paper. I think I was confusing the loadings and the PCs. The paper calls for only using the first PC. So the standardized index is the first PC divided by the square root of the first eigen value? Thanks for your assistance, I think I am getting a better understanding of the process undertaken by Messer now. Commented Feb 9, 2022 at 17:25