在R包geomorph中,如何解释$ rotation返回的PC负载矩阵?

问题描述 投票:0回答:1

我使用R中的地貌进行了二维几何形态分析。32个标本中有13个地标。

在地标上进行过程分析并执行PCA之后

lands<-readland.tps(file="2d_pelvis_GM.TPS", specID = "ID", negNA=TRUE)
lands<-estimate.missing(lands, method = "TPS")
gpa<-gpagen(lands)
gm2<-gm.prcomp(gpa$coords)

在gm.prcomp()文档中,它说$ rotation将给我“变量加载的矩阵,即分解后的矩阵的特征向量。”

但是,当我尝试此操作时,我得到:

gm2$rotation
              [,1]         [,2]         [,3]        [,4]
 [1,] -0.292637442 -0.247819659 -0.106858931  0.01109767
 [2,]  0.308917172 -0.393705793 -0.059492315 -0.07341983
 [3,]  0.030000517 -0.081826909 -0.065424375  0.09800728
 [4,]  0.052069888 -0.165765161 -0.182545759 -0.01158679
 [5,]  0.548194257 -0.381569139 -0.018172830 -0.01715879
 [6,] -0.006640823  0.509101670 -0.110753022  0.07658526
 [7,] -0.010698931 -0.016356517 -0.003072061  0.04752315
 [8,]  0.055072685  0.043019180 -0.011659494  0.10461193
 [9,] -0.104866710 -0.047445617 -0.016123797 -0.06791883
[10,] -0.022524580 -0.050072314  0.103334864  0.16119988
[11,] -0.229010887  0.016002306  0.348741743 -0.15581785
[12,] -0.037487586 -0.124647558  0.147662540  0.07513417
[13,]  0.181659233  0.257690024 -0.159775115 -0.34637649
[14,] -0.401315680 -0.177654634 -0.267710400 -0.13799427
[15,]  0.120250466  0.287037481 -0.160233065 -0.23500784
[16,] -0.308168182 -0.098020742 -0.316002315 -0.11907409
[17,] -0.201954619  0.002635231  0.103542233 -0.06862048
[18,]  0.014172165 -0.003263227  0.167454235  0.22778472
[19,] -0.139870351 -0.038388240 -0.043993696 -0.03608711
[20,] -0.060611865  0.111322129  0.389121041  0.18020552
[21,]  0.084694959  0.083795575  0.048581716  0.40593457
[22,]  0.084609928  0.040567815  0.521324332 -0.52082230
[23,]  0.025200160  0.191674102 -0.008007330  0.34023543
[24,]  0.237283438  0.264741813 -0.268730192 -0.14407308
[25,] -0.010960652 -0.025428638  0.080795508  0.02418930
[26,]  0.084623439  0.044376822 -0.112003514  0.18144889

...依此类推,一直到[,26]。自然地,我注意到该矩阵中的行数和列数是地标数量的两倍,但我不知所措,无法用它来解释哪些地标对每个PCA的贡献最大。当我在其他R PCA分析中查找相似的组件时,似乎给出了更直观的矩阵(行是地标,列是PC),因此它对我没有太大帮助。

非常感谢您的帮助

r pca
1个回答
0
投票

我正在经历同样的谜语。我所知道的是,此矩阵中的行和列数是地标数的两倍,因为其中一个(行或列)代表变量,即tps中地标的x坐标和y坐标格式。解释此矩阵的主要问题是弄清楚R如何布置变量。我猜这将是1x,1y,2x,2y,3x,3y,...

而且这也不能解释为什么列和行都是地标数的两倍。我(再次)猜测,列(可能)表示变量,行(可能)表示按解释的变化百分比排列的许多主成分,默认情况下将其设置为与变量相同。

我已复制您的问题,以在ResearchGate中问同样的问题。

我希望这是有用的。

© www.soinside.com 2019 - 2024. All rights reserved.