《非参数统计》与MATLAB编程 第八章 相关分析(3)

NaN NaN 1.0000 -0.9623 NaN NaN -0.9623 1.0000 rup =

NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN 1.0000 0.6619 NaN NaN 0.6619 1.0000 说明：NaN表示不存在数据，在x中，只有第四列和每五列每行的数据都存在，[r,p,rlo,rup]=corrcoef(x,'rows', 'all')

只计算第三列和第四列的相关系数及相应的区间。即： [r,p,rlo,rup]=corrcoef(x(:,3), x(:,4)) r =

1.0000 -0.5296 -0.5296 1.0000 p =

1.0000 0.3587 0.3587 1.0000 rlo =

1.0000 -0.9623 -0.9623 1.0000 rup =

1.0000 0.6619 0.6619 1.0000

命令：[r,p,rlo,rup]=corrcoef(x,'rows', 'complete') 如果某行含有N aN，则去掉所有含N aN的行，再计算相关系数和区间。

例：x =[ 3 9 2 34 5 NaN 6 66 NaN 4 9 7 8 7 33 9 9 14 11 22]

[r,p,rlo,rup]=corrcoef(x,'rows', 'complete') r =

1.0000 0.3883 0.6080 -0.7630 0.3883 1.0000 -0.4956 0.2995 0.6080 -0.4956 1.0000 -0.9771 -0.7630 0.2995 -0.9771 1.0000 p =

1.0000 0.7462 0.5840 0.4475 0.7462 1.0000 0.6699 0.8064 0.5840 0.6699 1.0000 0.1365

0.4475 0.8064 0.1365 1.0000 rlo =

1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 rup =

1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1

[r,p,rlo,rup]=corrcoef(x,'rows', 'complete')与去掉二、三行所得矩阵再求 corrcoef相同。即与： [r,p,rlo,rup]=corrcoef([x(1,:);x(4,:);x(5,:)]) r =

1.0000 0.3883 0.6080 -0.7630 0.3883 1.0000 -0.4956 0.2995 0.6080 -0.4956 1.0000 -0.9771 -0.7630 0.2995 -0.9771 1.0000 p =

1.0000 0.7462 0.5840 0.4475 0.7462 1.0000 0.6699 0.8064

0.5840 0.6699 1.0000 0.1365 0.4475 0.8064 0.1365 1.0000 rlo =

1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 rup =

1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1 NaN NaN NaN NaN 1

这里rlo,rup中显示NaN，是因为求区间要用到n-3,这里n-3=0

命令：[r,p,rlo,rup]=corrcoef(x,'rows', 'pairwise') r(i,j)的计算是把X中i或j列中含有NaN的行去掉，再计算i与j列的相关系数等指标。 例：x =[ 3 9 2 34 5 NaN 6 66 NaN 4 9 7 8 7 33 9

9 14 11 22]

[r,p,rlo,rup]=corrcoef(x,'rows', 'pairwise') r =

1.0000 0.3883 0.6476 -0.5791 0.3883 1.0000 -0.1687 0.5868 0.6476 -0.1687 1.0000 -0.5296 -0.5791 0.5868 -0.5296 1.0000 p =

1.0000 0.7462 0.3524 0.4209 0.7462 1.0000 0.8313 0.4132 0.3524 0.8313 1.0000 0.3587 0.4209 0.4132 0.3587 1.0000 rlo =

1.0000 NaN -0.8302 -0.9895 NaN 1.0000 -0.9722 -0.8584 -0.8302 -0.9722 1.0000 -0.9623 -0.9895 -0.8584 -0.9623 1.0000 rup =

1.0000 NaN 0.9915 0.8614 NaN 1.0000 0.9457 0.9897 0.9915 0.9457 1.0000 0.6619 0.8614 0.9897 0.6619 1.0000