2014-01-04 2 views
0

Я пытаюсь применить анализ основных компонентов в ковариационной матрице, оцененной по отношениям между всеми индивидами, см. Mm в следующем примере. я оценю его так много, если кто может показать мне, как это сделатьПрименение PCA на определенной ковариационной матрице пользователя

Пример:

library(BLR) 
library (rrBLUP) 
data(wheat) 

Мм ковариационная матрица, что я хочу, чтобы вычислить PCA

Mm <- A.mat(X)     
+1

Обратитесь к нам за помощью?? Princomp'. – Fernando

+0

и? Prcomp. Вот несколько [примеров] (http://manuals.bioinformatics.ucr.edu/home/R_BioCondManual#TOC-Principal-Component-Analysis-PCA-). – marbel

ответ

0

Попробуйте следующее используя ваш набор данных. Я не уверен, что это вы после

pca <- prcomp(Mm, scale=T) 

# Check out whats in it 
str(pca) 

# Prints variance summary for all principal components. 
summary(pca) 

# Accesses subset of components. 
summary(pca)$importance[, 1:20] 

require(GGally) 
ggpairs(pca$x[,1:4]) 
+0

Martín Bel: Благодарю вас за ответ. Но меня беспокоит то, что обработка ковариационной матрицы является регулярной. – hema

+0

@hema Почему это проблема? Что не хватает для этого, чтобы быть действительным ответом? – marbel

Смежные вопросы