我想通过使用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()本身并不能让你看到这些值,这是一个耻辱:/
1条答案
按热度按时间pqwbnv8z1#
箭头的坐标“只是”变量在双标图中所示的两个主成分上的载荷(加上一些比例缩放 *,以与同一图中的观察值得分相拟合)。
如果
rating.pc
是你的prcomp结果,你将从rating.pc$rotation
得到载荷。例如,具有四个主成分的prcomp结果pc
将返回以下载荷结构:......从该Angular 可以容易地计算和比较例如PC1 × PC2平面中的Angular 和长度。
stats:::biplot.prcomp