從軸標簽的摘要（rda）中提取“解釋的方差比例”

Question

我正在尋找 plot 一個 RDA，並讓軸標簽自動填充方差比例，該方差比例由摘要 output 解釋的方差比例解釋。 不幸的是，我似乎無法弄清楚如何調用特定的 object。 使用示例：

a<-as.matrix(varespec)
b<-as.matrix(varechem)
R<-rda(b,a)
summary(R)
R2<-RsquareAdj(R)
P2<-anova(R)
title<-paste("P=", P2$`Pr(>F)`, "R2=",R2$r.squared)
plot(R,type="n",bty="n",main= title,
     xlab="{var RDA1}% variance explained", 
     ylab="{var RDA2}% variance explained", 
     col.main="black",col.lab="black", col.axis="white",
     xaxt="n",yaxt="n")
abline(h=0,v=0,col="black",lwd=2)
points(R,display="species",col="gray",pch=20) 
points(R,display="cn",lwd=3)
text(R,display="cn",cex=0.5) 

在該示例中，我想將 {var RDA1} 替換為分別來自 RDA1 和 RDA2 列下的“組件的重要性：”部分的“比例解釋”的值（分別為第 1 列和第 2 行，第 2 行）。 我似乎無法弄清楚如何引用這些特定值。

Answer 1

這些值由summary(R)計算，因此您必須將它們保存為 object：

R.sum <- summary(R)
R.sum$cont   # Prints the "Importance of components" table
R.sum$cont$importance[2, "RDA1"]
# 0.74785
R.sum$cont$importance[2, "RDA2"]
# [1] 0.19804