一、实验数据
需要考察多个变量的边缘正态性
二、考察边缘正态性
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