偏相关分析

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问大侠，做偏相关分析时，控制变量的个数多了，为什么SAS或SPSS都没有结果输出？左偏相关有什么要注意的？谢谢啦！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：偏相关分析 相关分析 偏相关 控制变量 SPSS 偏相关

把SAS程序贴出来看看

data a;
input y x1-x14;
datalines;
0.0142         0.45         2.26         6.05         3.63         13.15         74.45         23.39         7.71         11.21         2.06         9.86         19.24         21.26         5.27
0.0132         1.02         5.82         24.92         11.44         31.60         25.20         38.80         15.03         23.25         6.36         3.94         6.70         3.89         2.02
0.0410         0.76         5.34         24.53         11.82         21.90         35.65         35.39         13.25         27.29         13.93         2.46         3.56         3.04         1.08
0.0121         2.26         8.55         32.95         14.04         20.75         21.45         34.23         15.02         25.53         6.80         3.63         7.40         5.50         1.90
0.0392         1.62         6.50         24.35         20.33         18.70         28.50         33.01         13.43         33.99         16.70         2.32         0.29         0.24         0.02
0.0129         4.76         16.11         40.13         12.00         10.05         16.95         36.45         15.24         40.51         6.51         0.51         0.39         0.40         0.00
0.0104         1.71         9.08         34.58         14.58         15.90         24.15         32.72         14.83         40.91         10.52         0.35         0.31         0.34         0.03
0.0091         2.79         9.68         37.61         13.27         9.20         27.45         32.68         13.87         40.03         12.94         0.18         0.17         0.11         0.00
0.0148         1.05         8.74         39.18         23.28         18.20         9.55         31.49         16.54         40.27         10.82         0.24         0.32         0.32         0.00
0.0150         2.68         9.37         44.18         18.61         12.65         12.50         38.20         15.15         36.33         8.15         0.43         0.46         0.77         0.51

run;

proc corr;

with y;
partial x1-x9;

run;

我用SAS不多，都是用的时候在书上找到相应代码。以上是代码，不知有什么问题。控制前x1-x7，SAS可以输出y与x8-x14的偏相关系数；控制x1-x8时，y与x9-x14的偏相关系数都是1，没有给显著水平值；像上面代码里控制x1-x9，SAS就不能给出y与x10-x14的偏相关系数及显著水平了。

不知是不是可以控制任意个变量，是我对偏相关分析的理解有误还是SAS代码有问题，向各位请教，谢谢！

1. proc corr data=kk;
   
2. var a b;
   
3. with c;
   
4. particl d1 d2 d3;
   
5. run;

复制代码

var 语句漏掉了

data a;
input y x1-x14;
datalines;
0.0142         0.45         2.26         6.05         3.63         13.15         74.45         23.39         7.71         11.21         2.06         9.86         19.24         21.26         5.27
0.0132         1.02         5.82         24.92         11.44         31.60         25.20         38.80         15.03         23.25         6.36         3.94         6.70         3.89         2.02
0.0410         0.76         5.34         24.53         11.82         21.90         35.65         35.39         13.25         27.29         13.93         2.46         3.56         3.04         1.08
0.0121         2.26         8.55         32.95         14.04         20.75         21.45         34.23         15.02         25.53         6.80         3.63         7.40         5.50         1.90
0.0392         1.62         6.50         24.35         20.33         18.70         28.50         33.01         13.43         33.99         16.70         2.32         0.29         0.24         0.02
0.0129         4.76         16.11         40.13         12.00         10.05         16.95         36.45         15.24         40.51         6.51         0.51         0.39         0.40         0.00
0.0104         1.71         9.08         34.58         14.58         15.90         24.15         32.72         14.83         40.91         10.52         0.35         0.31         0.34         0.03
0.0091         2.79         9.68         37.61         13.27         9.20         27.45         32.68         13.87         40.03         12.94         0.18         0.17         0.11         0.00
0.0148         1.05         8.74         39.18         23.28         18.20         9.55         31.49         16.54         40.27         10.82         0.24         0.32         0.32         0.00
0.0150         2.68         9.37         44.18         18.61         12.65         12.50         38.20         15.15         36.33         8.15         0.43         0.46         0.77         0.51

run;

proc corr;
         
var x10-x14;
with y;
partial x1-x9;

run;

加上了var还是没有偏相关系数和显著水平结果，log窗口里有4行警告：
WARNING: The variable x10 is singular in calculating the partial PEARSON correlations.
WARNING: The variable x11 is singular in calculating the partial PEARSON correlations.
WARNING: The variable x12 is singular in calculating the partial PEARSON correlations.
WARNING: The variable x13 is singular in calculating the partial PEARSON correlations.
WARNING: The variable x14 is singular in calculating the partial PEARSON correlations.
不知是哪出问题了

可以先控制1个变量试试

我从x1开始逐渐增加控制的个数，一直增加到x7时正常，到x8时y与其余变量偏相关系数都是1，没有显著水平结果，到x9及以上就算不出偏相关系数和显著水平了

数据观测次数太少啦，直接就奇异阵出来了。
偏相关系数是控制其它变量对相关系数考虑变量的影响后的一个相关系数，你总共就只有10次观测，然后基于9个变量去控制，对于每十个变量来讲，一个10×10的矩阵很容易就奇异了，最后，所有的考察变量的变异都被9个控制变量解释了，最后压根就没有什么能用来做相关系数的计算。

样本量问题