多元统计分析-----例8.1:今有14名学生的身高和体重数据,做相关图以显示相关变量间的关系。

例8.1:今有14名学生的身高和体重数据,做相关图以显示相关变量间的关系。

x1=c(147,171,175,159,155,152,158,154,164,168,166,159,164,177) 

x2=c(32,57,64,1,38,35,44,41,54,57,49,47,46,63)

plot(x1,x2,xlim = c(145,180),ylim=c(25,75))

出错:

lines(c(146,178),c(30,66),lwd=2)

lines(c(163,166),c(54,47))

library(shape)

解决方法如下:

1.安装该包:install.packages("shape")

library(shape);

lines(getellipse(24,3,mid=c(162,48),angle = 48),lty=3)

        正如二维椭圆有两个主轴,三维桶球有三个主轴一样,有几个变量,就有几个主成分。当然,选择越少的主成分,降维就越好。那么,什么是标准呢?那就是这些被选的主成分所代表的主轴的长度之和占了主轴长度总和的大部分。有些学者建议,所选的主轴总长度占所有主轴长度之和的80%左右即可,其实,这只是一个大体的说法;具体选几个,要根据实际情况而定。但如果所有涉及的变量都不那么相关,就很难降维。

X=data.frame(x1,x2) S=cov(X);S #协方差阵

x1 x2 x1 77.45604 95.31868 x2 95.31868 257.36264

R=cor(X);R #相关系数阵

x1 x2 x1 1.0000000 0.6751148 x2 0.6751148 1.0000000

主成分的性质:

主成分分析函数

princomp(x,cor = FAL.SE,scores = TRUE,-…)

princomp(X,cor = FALSE,scores = TRUE) #主成分分析函数

x为数据矩阵或数据框;

cor表示是否用相关系数阵,默认为协方差阵;scores表示是否输出成分得分

pc=princomp(X);pc

names(pc) #主成分的对象

[1] "sdev" "loadings" "center" "scale" "n.obs" "scores" "call"

pc $ sdev^2 #主成分方差

Comp.1 Comp.2 277.15195 33.75111

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

墨染盛夏呀

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值