If the set of data points are orthogonal

> library(mvtnorm)
> x <- rmnorm(1000, mean = c(0, 6), varcov = diag(6))
> round(eigen(cov(x))$values)
[1] 1 1 1 1 1 1

Eigen values are all one and eigen vectors are orthonormal