从因子分析中消除R中的箭头

Question

在下面的代码生成的biplot中，我试图去掉红线。如果有人能帮忙我会很感激的。

library(psych)
data(bfi)
fa2 <- fa(bfi[16:25],2)  #factor analysis
fa2$scores <- fa2$scores[1:100,]  #just take the first 100
biplot(fa2,pch=c(24,21)[bfi[1:100,"gender"]],bg=c("blue","red")[bfi[1:100,"gender"]],
main="Biplot of Conscientiousness and Neuroticism by gender")

​

​

Answer 1

无论出于什么原因，psych库决定重写它自己的双图，因此它忽略了许多标准参数。您可以创建自己的版本，只需删除箭头绘图。这个方法有点麻烦，但是用psych_1.4.5进行了测试。只要确认一下

body(biplot.psych)[[c(11,3,12)]]

返回

arrows(0, 0, x$loadings[, 1L] * 0.8, x$loadings[, 2L] * 0.8, 
    col = col[2L], length = arrow.len)

以确保我们正在更改正确的行。那你就可以

biplot.psych2<-biplot.psych
body(biplot.psych2)[[11]][[3]][[12]]<-NULL

然后调用我们的新函数

biplot.psych2(fa2,pch=c(24,21)[bfi[1:100,"gender"]],
    bg=c("blue","red")[bfi[1:100,"gender"]],
    main="Biplot of Conscientiousness and Neuroticism by gender")

要获得

​

​

Answer 2

在你贴出问题几个月后，这里有一个简短的答案。实际上，您要求的不是双情节(包括因素分数和因素负载)，而是要素分数的图表。用你的例子，只是

Bg=c(“蓝色”，“红色”)[bfi1:100，“性别”]，pch=c(24,21)[bfi1:100，“性别”]，ylim=c(-3,3)，xlim=c(-3,3))。