R语言 获取由双标图()创建的红色箭头的值

blmhpbnm  于 2023-06-19  发布在  其他
关注(0)|答案(1)|浏览(93)

我想通过使用prcomp()函数确定bipot()创建的红色箭头的值。我想确定箭头的长度和位置,以便更准确地比较它们。下面是我使用的代码:

df <- read.csv(file = "data/PCA_data_all.csv", sep = ",")
df <- df %>%
  separate(model_name, into = c("Make"), sep = " ", remove = FALSE)
df <- df %>%
  select(-model_name)
df[,2:31] <- scale(df[,2:31])
brand.mean <- aggregate(df[,2:31], list(df[,1]), mean)
row.names(brand.mean) <- brand.mean$Group.1
brand.mean <- brand.mean[,-1]
rating.pc <- prcomp(brand.mean, scale=TRUE)
plot(rating.pc, type="l")
biplot(rating.pc, main="Brand Positioning", cex=c(0.5, 0.5))

据我所知,biplot()本身并不能让你看到这些值,这是一个耻辱:/

pqwbnv8z

pqwbnv8z1#

箭头的坐标“只是”变量在双标图中所示的两个主成分上的载荷(加上一些比例缩放 *,以与同一图中的观察值得分相拟合)。
如果rating.pc是你的prcomp结果,你将从rating.pc$rotation得到载荷。例如,具有四个主成分的prcomp结果pc将返回以下载荷结构:

> pc$rotation
                     PC1         PC2         PC3        PC4
Sepal.Length  0.36138659 -0.65658877  0.58202985  0.3154872
Sepal.Width  -0.08452251 -0.73016143 -0.59791083 -0.3197231
Petal.Length  0.85667061  0.17337266 -0.07623608 -0.4798390
Petal.Width   0.35828920  0.07548102 -0.54583143  0.7536574

......从该Angular 可以容易地计算和比较例如PC1 × PC2平面中的Angular 和长度。

  • 有关详细信息,请输入stats:::biplot.prcomp

相关问题