2.6 相关与回归分析的SAS程序
在阅读以下内容之前,请先阅读第一章“SAS软件基本操作”。
2.6.1 一元线性回归分析
例 2.16 以课本上表10-1中的数据为例,求出一元回归方程、检验回归显著性并求出回归及预测值的0.95置信区间。
解: 使用PROC REG过程进行分析,SAS程序如下:
options linesize = 76;
data soil;
input salt dw @@;
cards;
0 80 0.8 90 1.6 95 2.4 115 3.2 130 4.0 115 4.8 135
proc reg;
model dw = salt;
run;
输出结果见表2-20。
表2-20 例2.16所输出的结果
The SAS System
Model: MODEL1
Dependent Variable: DW
Analysis of Variance
|
|
Sum of |
Mean |
|
|
Source |
DF |
Squares |
Square |
F Value |
Prob>F |
|
|
|
|
|
|
Model |
1 |
2232.14286 |
2232.14286 |
31.566 |
0.0025 |
Error |
5 |
353.57143 |
70.71429 |
|
|
C Total |
6 |
2585.71429 |
|
|
|
Root MSE |
8.40918 |
R-square |
0.8633 |
Dep Mean |
108.57143 |
Adj R-sq |
0.8359 |
C.V. |
7.74530 |
|
|
Parameter Estimates
|
|
Parameter |
Standard |
T for H0: |
|
Variable |
DF |
Estimate |
Error |
Parameter=0 |
Prob>|T| |
|
|
|
|
|
|
INTERCEP |
1 |
81.785714 |
5.72988941 |
14.274 |
0.0001 |
SALT |
1 |
11.160714 |
1.98648174 |
5.618 |
0.0025 |
表2-20的第一部分是对回归所做的方差分析,第二部分给出了截距(表中的INTERCEP),即回归方程中的常数项a,和回归系数b(表中的SALT)。可以得出回归方程:
并给出在H0: a(b) = 0下,对a和b所做的t检验。
为了得到残差和置信区间,可以将过程步做以下补充:
proc sort out = sorted;
by salt;
run;
proc reg data = sorted;
model dw = salt / r clm;
id salt;
来源:不详