ラスターのスタックで PCA を実行した後 (これと2014 Raster Package のドキュメントに類似)、固有値、固有ベクトル、および負荷を確認したいと思います...
スクリー プロット、変動の比率、累積比率 (summary()、print()、hist()) を返す princomp の一般的な呼び出しは、RasterBrick 出力から情報を取得していないようです。コード例は次のとおりです。
#bring in rasterBrick
logo <- brick(system.file("external/rlogo.grd", package="raster"))
#select random samples & do PCA
sr <- sampleRandom(logo, 100)
pca <- princomp(sr)
# to visualize pcs as rasters
x <- predict(logo, pca, index=1:3)
plot(x)
##ANSWERED QUESTION
summary(pca) # importance of components
plot (pca) # scree plot
loadings (pca) #eigens
summary()
解析からの値ではなく、ラスター レイヤーの要約統計と思われるものを返します。print()
ラスターなどの最小値と最大値を示します。
特にPCAに関連する固有値と固有ベクトルを見つける方法についての明確さについて、ご意見をお寄せいただきありがとうございます。