Разные результаты при запуске PCA в R или Stata - PullRequest
Новая
       45

Разные результаты при запуске PCA в R или Stata

0 голосов
/ 10 июля 2020

Я запускаю анализ основных компонентов в R, но значения, которые я получаю для оценок основного компонента 1 (PC1), противоположны (т.е. умножены на -1) тому, что я получил работает в Stata.

В Stata я использую следующие команды (после загрузки фрейма данных ниже): 

pca X1 X5 X7

(где X1, X5 и X7 - столбцы моего данные ) И: 

predict pc1 pc2 pc3, score

В R я использую команду prcomp из пакета factoextra: 

daf <- prcomp(x=d, scale. = TRUE, retx = TRUE)
daf$x[,1]

Первый элемент оценки ПК1 является отрицательным при использовании R, но положительным при использовании Stata, и то же самое происходит с другими значениями. Я неправильно использую PCA в R? Я был бы благодарен за любую помощь!

Ниже я поместил фрейм данных, который использовал. Прошу прощения за то, что не опубликовал воспроизводимый пример меньшего размера. Прошу прощения за свой англи sh тоже ...

Данные: 

    d<-data.frame(
  X1 = c(1,0,1,0,1,0,1,0,0,0,0,0,0,1,1,0,0,0,1,1,1,0,0,1,0,1,1,0,1,1,0,0,1,1),
  X5 = c(0.723264853549338, 0.293095848783131, 0.249530098571494, 0.568117582946663, 0.893216842616187, 0.458473806968546, 0.540790893374034, 0.310274194098472, 0.439974636572571, 0.283240938166311, 0.305186090937202, 0.477321362384237, 1, 0.723986731593571, 0.720605502538762, 0.587581748710709, 0.387998819832435, 0.41558154454738, 0.606024668782451, 0.595836088744153, 0.421111284523995, 0.551228660524501, 0.594579624892775, 0.585079156569203, 0.173307415389415, 0.416460554371002, 0.775289506275296, 0.515752968727375, 0.432484109191689, 0.283780145509236, 0.520817800909547, 0.560461124118364, 0, 0.396384701157797),
  X7 = c(0.584523240443511, 0.990811451154001, 0.944206882738442, 0.676674812375392, 0.835035924324797, 1, 0.804072676452289, 0.983723077483278, 0.377087978288423, 0.544645005558396, 0.975878431927861, 0.758974140733694, 0.551475086402936, 0.472944754460553, 0.3681911072948, 0.433919200876731, 0.653951942097384, 0.873974253793443, 0.655810896249537, 0.550689284075584, 0.934784347180246, 0.911025962497841, 0.90141604811194, 0.912158435125721, 0.880753164580567, 0.878081585424683, 0.176583243843766, 0.914869501586647, 0.940970117722151, 0.675276216278826, 0.894845196341543, 0.88948449388328, 0, 0.922911045673136)
)
...