R里面主成分分析的双变量图怎么画比较好看

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有哪位知道主成分分析中关于第1主成分和第2主成分所画出来的双变量图怎么画比较好看的吗？比如说，上面的数字变成样本的名称，坐标轴以0为界分成四个象限这样子的。（我把代码贴在下面了），谢谢~~~麻烦各位了！
【相关程序如下：

>a=read.csv("shanghai2010.csv",header=T)

> attach(a)

> a[,3:8]    #提取a中第3-8列

> PCA=princomp(a[,3:8],cor=T)  #主成分分析

> PCA   #特征值

> PCA$loadings      #主成分载荷

> screeplot(PCA,type="lines")  #主成分碎石图

> PCA$scores   #主成分得分

>biplot(PCA)   #画双变量图】

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：主成分分析 主成分 双变量 Loadings princomp scores 坐标轴 程序 样本 主成分分析法 spss主成分分析 逐步回归分析 多元回归分析 因子分析法 应用时间序列分析

R好像就是这样的。其他软件也许不是。比如Minitab

tmdxyz 发表于 2012-4-7 19:29
R好像就是这样的。其他软件也许不是。比如Minitab

谢谢你！不过这样看了半天，感觉太丑了

丑不丑，那是决定于你的数据呀。

昨天无意间看到了：你在画完biplot之后，再加上以下两个语句，运行即可。

lines(c(-10,10),c(0,0))

lines(c(0,0),c(-10,10))

数字要变成样本的名称得重新对行和列命名，用rownames(a)=("","") colnames(a)=("")

抠脚小仙 发表于 2014-7-10 21:03
数字要变成样本的名称得重新对行和列命名，用rownames(a)=("","") colnames(a)=("")

具体怎么弄能说一下么？不胜感激

可以试试ggbiplot包，调用ggbiplot()函数