利用Matlab考察数据的边缘正态性 作Q_Q图

一、实验数据

需要考察多个变量的边缘正态性

二、考察边缘正态性

1.原理

可以通过对每个变量的样本观测值的边缘分布作Q_Q图和计算相关系数 r_Q 来考察。
Matlab自带有作Q_Q图的函数qqplot() ，但是并没有计算相关系数 r_Q 的函数。参考实用多元统计分析的4.6节，用其中作Q_Q图的步骤和 r_Q 的公式自编Matlab程序进行求解，并且封装成函数。自编程序中用到了Matlab的norminv（probability,mean,standard_dev）函数，令mean = 0，standard_dev = 1 可以得到对应概率probability的标准正态分布的下侧分位数。

function[r]= Q_Q (x)
x1=sort(x); %对原矩阵每列进行升序排列
[n,m]=size(x1);  %返回维数
p=ones(n