- 1、本文档共27页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《非参数统计》与MATLAB编程相关分析
第八章 相关分析
协方差
命令:C = cov(X)
当X为行或列向量时,它等于 var(X) 样本标准差。
X=1:15;cov(X)
ans =
20
var(X)
ans =
20
当X为矩阵时,此时X的每行为一次观察值,每列为一个变量。cov(X)为协方差矩阵,它是对称矩阵。
例:x=rand(100,3);c= cov(x)
c=
0.089672 -0.012641 -0.0055434
-0.012641 0.07928 0.012326
-0.0055434 0.012326 0.082203
c的对角线为:diag(c)
ans =
0.0897
0.0793
0.0822
它等于:var(x)
ans =
0.0897 0.0793 0.0822
sqrt(diag(cov(x)))
ans =
0.2995
0.2816
0.2867
它等于:std(x)
ans =
0.2995 0.2816 0.2867
命令:c = cov(x,y)
其中x和y是等长度的列向量(不是行向量),它等于cov([x y])或cov([x,y])
例:x=[1;4;9];y=[5;8;6];
c=cov(x,y)
c =
16.3333 1.1667
1.1667 2.3333
cov([x,y])
ans =
16.3333 1.1667
1.1667 2.3333
COV(X)、 COV(X,0)[两者相等] 或COV(X,Y)、COV(X,Y,0) [两者相等],它们都是除以n-1,而COV(X,1) or COV(X,Y,1)是除以n
x=[1;4;9];y=[5;8;6];
cov(x,y,1)
ans =
10.8889 0.7778
0.7778 1.5556
它的对角线与var([x y],1) 相等
ans =
10.8889 1.5556
协差阵的代数计算:
[n,p] = size(X);
X = X - ones(n,1) * mean(X);
Y = X*X/(n-1); Y为X的协差阵
相关系数(一)
命令:r=corrcoef(x)
x为矩阵,此时x的每行为一次观察值,每列为一个变量。
r为相关系数矩阵。它称为Pearson相关系数
例:x=rand(18,3);r=corrcoef(x)
r =
1.0000 0.1509 -0.2008
0.1509 1.0000 0.1142
-0.2008 0.1142 1.0000
r为对称矩阵,主对角阵为1
命令:r=corrcoef(x,y)
其中x和y是等长度的列向量(不是行向量),它等于cov([x y])或cov([x,y]),或x和y是等长度的行向量,r=corrcoef(x,y)它则等于r=corrcoef(x’,y’), r=corrcoef([x’,y’])
例:x=[1;4;9];y=[5;8;6]; corrcoef(x,y)
ans =
1.0000 0.1890
0.1890 1.0000
corrcoef([x,y])
ans =
1.0000 0.1890
0.1890 1.0000
C = COV(X)
Rij =C(i,j)/SQRT(C(i,i)*C(j,j))
如:X=[1 2 7 4 ;5 12 7 8;9 17 11 17];
cov(X)
ans =
16.0000 30.0000 8.0000 26.0000
30.0000 58.3333 13.3333 46.6667
8.0000 13.3333 5.3333 14.6667
26.0000 46.6667 14.6667 44.3333
corrcoef(X)
ans =
1.0000 0.9820 0.8660 0.9762
0.9820 1.0000 0.7559 0.9177
0.8660 0.7559 1.0000 0.9538
0.9762 0.9177 0.9538 1.0000
则有:30/sqrt(16*58.3333)
ans =
0.9820
命令:[r,p]=cor
文档评论(0)