"глвные компоненты" = "главные оси эллипсоида инерции" (вроде как)
Разве это противоречит тому, что я написал?
Я бы сказал "не очень похоже"... особенно если не нормализовывать данные каждого месяца...
Первая компонента очень похожа (нормализация входит в понятие "по сути"):
Код:
> cor(rowMeans(x[,1:12]),pr.out$x[,1])
[1] 0.9828019
> cor(rowMeans(x[,13:24]),pr.out$x[,1])
[1] 0.9796792
А вот вторую я назвал неточно.
Код:
> cor(rowMeans(x[,5:11])-rowMeans(x[,c(1:4,12)]),pr.out$x[,2])
[1] 0.9219097
То есть это разница между летом и зимой. Со знаком.
И спасибо за совет почистить от нулей.
Код:
r=r[r[[12]]!=0 & r[[13]]!=0,]