r - グループ/クラスに基づく PCA の機能の貢献

Question

各タイプのワイン クラス (バローロ、グリニョリーノ、バルベーラ) の機能の寄与を取得/プロットしたいと思います。fviz_contribを使用すると、以下の MWE に示すように、すべてのクラスの貢献度を取得できます。しかし、クラス/グループで個別にフィルタリングして計算/プロットできるかどうか、またどのようにできるか疑問に思っていました。

library(ggbiplot)
library(factoextra)

data(wine)

wine.pca <- prcomp(wine, scale. = TRUE)

# plot the PCA 
print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, groups = wine.class, ellipse = TRUE, circle = TRUE))

# plot the contributions of the features for all wine classes
g.contr <- fviz_contrib(wine.pca, choice = "var", axes = 1:2, fill = "lightblue", color = "darkblue", top = 45)
print(g.contr)