回歸分析程式第1頁，課件共51頁，創(chuàng)作于2023年2月回歸分析程式概述常用的回歸分析程式過程有:(1)CORR過程、（2）REG過程、(3)NLIN過程、（4）RSREG過程?；貧w分析第2頁，課件共51頁，創(chuàng)作于2023年2月CORR過程1REG過程2NLIN過程3RSREG過程4ＳＡＳ的應(yīng)用基礎(chǔ)回歸分析第3頁，課件共51頁，創(chuàng)作于2023年2月CORR過程CORR過程可用來分析變量之間的相關(guān)關(guān)系。

CORR過程的格式如下：

PROCCORR[選擇項(xiàng)]；

VAR變量表；

PARTIAL變量表；在PROCCORR過程中可能出現(xiàn)的選擇項(xiàng)是DATA=數(shù)據(jù)集名，用來指定CORR的數(shù)據(jù)集名，如果省略這一選擇，則指定最新建立的數(shù)據(jù)集。

PARTIAL語句用來求偏相關(guān)系數(shù)，指定保持不變的偏變量名。回歸分析第4頁，課件共51頁，創(chuàng)作于2023年2月CORR過程例題6.1，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputx1-x4@@；Cards；47.254.7416.747.253.3754.58.56.791.745.74.4166.7692250.54.252250.542.57526350253.7166.70.54350142.575275300；proccorr；varx1-x4；proccorrnosimple；varx1x2；partialx3x4；proccorrnosimple；varx1x3；partialx2x4；run；回歸分析第5頁，課件共51頁，創(chuàng)作于2023年2月CORR過程程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：甘薯實(shí)生苗栽培試驗(yàn)中薯塊重x1(g)、塊根粗x2(mm)、單株結(jié)薯數(shù)x3及單株產(chǎn)量y(g)的12組觀測(cè)值如下表（引自溫奎課程論文，1987），試計(jì)算三種相關(guān)系數(shù)并作相關(guān)系數(shù)顯著性檢驗(yàn)？程式中的nosimple將抑制簡(jiǎn)單統(tǒng)計(jì)量的輸出?；貧w分析i123456789101112x1ix2iX3iyi4.04.04.54.06.00.50.52.02.00.51.02.07.257.258.55.79.04.254.06.05.04.04.07.04.73.36.74.42.02.02.53.03.73.02.55.0416.775.091.7166.725.025.075.050.0166.750.075.0300.0第6頁，課件共51頁，創(chuàng)作于2023年2月CORR過程本程式運(yùn)行的結(jié)果如下:

(圖一)回歸分析第7頁，課件共51頁，創(chuàng)作于2023年2月CORR過程

(圖二)回歸分析第8頁，課件共51頁，創(chuàng)作于2023年2月CORR過程

(圖三)回歸分析第9頁，課件共51頁，創(chuàng)作于2023年2月REG過程

REG過程可用來建立因變量與自變量之間的線性回歸方程。

REG過程的格式如下：

PROCREG[選擇項(xiàng)]；

MODEL因變量=自變量/[選擇項(xiàng)]；

在PROCREG過程中可能出現(xiàn)的選擇項(xiàng)是DATA=數(shù)據(jù)集名，用來指定要做REG的數(shù)據(jù)集名，如果省略這一選擇，則指定最新建立的數(shù)據(jù)集。

回歸分析第10頁，課件共51頁，創(chuàng)作于2023年2月REG過程例題6.2，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputxy@@；cards；14.346.31430.769.3144.622.769.27.316812.31.32.77.926.3；procreg；modely=x/xpxi；run；回歸分析第11頁，課件共51頁，創(chuàng)作于2023年2月REG過程程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：棉花紅鈴蟲第一代產(chǎn)卵高峰日百株卵量x(粒)與百株累計(jì)卵量y(粒)的8組觀測(cè)數(shù)據(jù)如下表，試建立二者間的一元線性回歸方程并作回歸方程的顯著性檢驗(yàn)？程式中的xpx及i指定輸出離均差矩陣及其逆矩陣。回歸分析i

12345678xi

14.314.069.322.77.38.01.37.9yi

46.330.7144.669.216.012.32.726.3第12頁，課件共51頁，創(chuàng)作于2023年2月REG過程本程式運(yùn)行的結(jié)果如下:

(圖一)回歸分析第13頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖二)回歸分析第14頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖三)回歸分析第15頁，課件共51頁，創(chuàng)作于2023年2月REG過程回歸分析

(圖四)第16頁，課件共51頁，創(chuàng)作于2023年2月REG過程例題6.3，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputx1x2y@@；cards；1.379.084.9311.341.891.869.673.062.330.7610.25.7817.670.050.0615.910.730.4315.7410.30.875.416.253.86；procreg；modely=x1x2/xpxi；run；回歸分析第17頁，課件共51頁，創(chuàng)作于2023年2月REG過程程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：某品種水稻糙米含鎘量y(mg/kg)與地上部生物量x1(10g/盆)及土壤含鎘量x2(100mg/kg)的8組觀測(cè)數(shù)據(jù)如下表，試建立二元線性回歸方程并作回歸方程的顯著性檢驗(yàn)？程式中的xpx及i指定輸出離均差矩陣及其逆矩陣。回歸分析i

12345678x1i

1.3711.349.670.7617.6715.9115.745.41x2i

9.081.893.0610.200.050.731.036.25yi

4.931.862.335.780.060.430.873.86第18頁，課件共51頁，創(chuàng)作于2023年2月REG過程本程式運(yùn)行的結(jié)果如下:(圖一)回歸分析第19頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖二)回歸分析第20頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖三)回歸分析第21頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖四)回歸分析第22頁，課件共51頁，創(chuàng)作于2023年2月REG過程例題6.4，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputxy@@；x1=x；x2=x**2；x3=x**3；cards；1.23.9141.53.722.82.22.62.61.4303.5-2.22.80.8；procreg；modely=x1x2；procreg；modely=x1-x3；run；回歸分析第23頁，課件共51頁，創(chuàng)作于2023年2月REG過程程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：假設(shè)變量x與y的9組觀測(cè)值如下表，試建立線性回歸方程及非線性回歸方程。程式中的modely=x1x2及modely=x1-x3指定建立一元二次拋物線回歸及一元三次拋物線回歸方程。回歸分析i

1

2

3

4

5

6

7

8

9x1.2

1

1.5

2

2.22.6

33.52.8y3.9

4

3.7

2.82.61.4

03.52.8第24頁，課件共51頁，創(chuàng)作于2023年2月REG過程本程式運(yùn)行的結(jié)果如下:

(圖一)回歸分析第25頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖二)回歸分析第26頁，課件共51頁，創(chuàng)作于2023年2月REG過程例題6.5，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputx1-x5y@@；cards；9.22.7321.4710.3321.1381.1159.13.7321.820.1120.8281.1468.64.8821.8270.3832.1311.84110.2333.9681.5870.1811.3491.3565.63.7321.8410.2971.8150.8635.3674.2361.8730.0631.3520.9036.1333.1461.9870.281.6470.1148.24.6461.6150.3794.5650.898回歸分析第27頁，課件共51頁，創(chuàng)作于2023年2月REG過程8.84.3781.5430.7442.0731.937.63.8641.5990.3422.4231.1049.74.3781.6910.2251.5221.4038.3675.0951.8070.4392.2811.76312.1674.8941.7280.1261.5811.63610.2673.7321.6120.1091.1981.4678.94.4721.880.0790.7950.9198.2335.2781.730.3033.0691.515；procreg；modely=x1-x5/selection=stepwise；run；回歸分析第28頁，課件共51頁，創(chuàng)作于2023年2月REG過程程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：在農(nóng)作物害蟲發(fā)生趨勢(shì)的預(yù)報(bào)研究中，所涉及的5個(gè)自變量及因變量的16組干測(cè)數(shù)據(jù)如下表，試建立逐步回歸方程。程式中的selection=stepwise指定用逐步回歸法建立線性回歸方程?；貧w分析

ix1ix2ix3ix4ix5iyi

19.2002.7321.4710.3321.1381.155

29.1003.7321.8200.1120.8281.146

38.6004.8821.8270.3832.1311.841

410.2333.9681.5870.1811.3491.356

55.6003.7321.8410.2971.8150.863

65.3674.2361.8730.0631.3520.903第29頁，課件共51頁，創(chuàng)作于2023年2月REG過程回歸分析

76.1333.1461.9870.2801.6470.114

88.2004.6461.6150.3794.5650.898

98.8004.3781.5430.7442.0731.930

107.6003.8641.5990.3422.4231.104

119.7004.3781.6910.2251.5221.403

128.3675.0951.8070.4392.2811.763

1312.1674.8941.7280.1261.5811.636

1410.2673.7321.6120.1091.1981.467

158.9004.4721.8800.0790.7950.919

168.2335.2781.7300.3033.0691.515第30頁，課件共51頁，創(chuàng)作于2023年2月REG過程本程式運(yùn)行的結(jié)果如下:(圖一)回歸分析第31頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖二)回歸分析第32頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖三)回歸分析第33頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖四)回歸分析第34頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖五)回歸分析第35頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖六)回歸分析第36頁，課件共51頁，創(chuàng)作于2023年2月REG過程

(圖七)回歸分析第37頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程N(yùn)LIN過程可用來建立因變量與自變量之間的非線性回歸方程。

NLIN過程的格式如下：

PROCNLIN[選擇項(xiàng)]；

MODEL因變量=表達(dá)式；

PARMS參數(shù)名=數(shù)值；

DER.參數(shù)名=表達(dá)式；在PROCNLIN過程中可能出現(xiàn)的選擇項(xiàng)是DATA=數(shù)據(jù)集名，用來指定要做NLIN的數(shù)據(jù)集名，如果省略這一選擇，則指定最新建立的數(shù)據(jù)集?；貧w分析第38頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程例題6.6，進(jìn)入SAS的顯示管理系統(tǒng)，在編輯窗口輸入程式：dataex；inputxy@@；cards；11.8521.3731.0240.7540.5660.4160.3180.2380.17；procnlin；modely=a*exp(b*x)；parmsa=2.5168b=-0.3221；der.a=exp(b*x)；der.b=a*x*exp(b*x)；run；回歸分析第39頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程

程式中的數(shù)據(jù)來自以下試驗(yàn)結(jié)果及問題：假設(shè)變量x與y的9組觀測(cè)值如下表，試選用多個(gè)一元非線性回歸方程進(jìn)行擬合，并且比較各個(gè)回歸方程的你和情況。

回歸分析i

1

2

3

4

5

6

7

8

9xi

1

2

3

4

4

6

6

8

8yi1.851.371.020.750.560.410.310.230.17第40頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程本程式運(yùn)行的結(jié)果如下:

(圖一)回歸分析第41頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程

(圖二)回歸分析第42頁，課件共51頁，創(chuàng)作于2023年2月NLIN過程(圖三)回歸分析第43頁，課件共51頁，創(chuàng)作于2023年2月RSREG過程RSREG過程可用來建立二次響應(yīng)面的回歸方程。

RSREG過程的格式如下：

PROCRSREG[選擇項(xiàng)]；

MODEL因變量=自變量；在PROCRSREG過程中可能出現(xiàn)的選擇項(xiàng)是DATA=數(shù)據(jù)集名，用來指定要做RSREG的數(shù)據(jù)集名，如果省略這一選擇，則指定最新建立的數(shù)據(jù)集?；貧w分析第44頁，課件共51頁，創(chuàng)作于2023年2月RSREG過程例題6.7，進(jìn)入SAS的顯示管理系