SAS實驗-相關(guān)與回歸分析

學生實驗報告學生實驗報告學生姓名學號同組人實驗工程相關(guān)與回歸分析必修操作性實驗實驗地點0505實驗儀器臺號指導教師董成武實驗日期及節(jié)次一、實驗目的及要求：1、目的學會和掌握參數(shù)估計與假設(shè)檢驗2、內(nèi)容及要求學會用SAS作總體均值、方差估計，以及單樣本均值、方差比例檢驗，雙樣本均值、方差、比例的比擬和配對樣本均值比擬檢驗等操作二、儀器用具：儀器名稱規(guī)格/型號數(shù)量備注計算機1SAS軟件1三、實驗方法與步驟:步驟一、運行SAS軟件，新建各題要用到的數(shù)據(jù)集；步驟二、新建數(shù)據(jù)集后，輸入analyst命令或在點擊解決方案—分析—分析家，進入分析員應用環(huán)境；步驟三、在分析員應用環(huán)境翻開新建的數(shù)據(jù)集；步驟四、再根據(jù)各題要求，點擊統(tǒng)計選項來分析。四、實驗結(jié)果與數(shù)據(jù)處理：1．下面是某地一年級的12位女生的體重(kg)和肺活量(L)的數(shù)據(jù)試求肺活量對體重的回歸方程。體重424246464650505052525858肺活量2.552.22.752.42.82.813.413.13.462.853.53〔1〕對肺活量和體重進行相關(guān)分析。繪制體重與肺活量的散點圖，計算其相關(guān)系數(shù)。并對輸出結(jié)果進行解釋。①將數(shù)據(jù)保存至EXCEL文件，導入數(shù)據(jù)，得到數(shù)據(jù)集；②在INSIGHT環(huán)境中調(diào)入數(shù)據(jù)，繪制體重與肺活量的散點圖，如下圖。③在INSIGHT環(huán)境中利用多元進行分析，得出體重與肺活量的相關(guān)系數(shù)?！皢巫兞拷y(tǒng)計量〞中提供了每個變量的5個統(tǒng)計量，分別為數(shù)量、均值、標準差、最小值和最大值?！跋嚓P(guān)系數(shù)舉證〞中那么提供了變量間的相關(guān)系數(shù)，體重〔weight〕和肺活量(pulmonary)的相關(guān)系數(shù)為0.7495〔2〕判斷體重和肺活量兩變量的關(guān)系，擬合肺活量對體重的方程，對回歸系數(shù)進行解釋和對回歸系數(shù)進行檢驗，繪制殘差圖。①用insight視圖下的分析來擬合肺活量對體重的方程；得到簡單線性回歸方程為對回歸系數(shù)進行解釋和對回歸系數(shù)進行檢驗：上表提供模型擬合的匯總度量，肺活量的均值是2.9025，模型的,說明肺活量變異的大約56%可以由體重來說明。從方差分析表中可以看到，P值小于，所以拒絕原假設(shè)接受備擇假設(shè)，所以這一模型擬合數(shù)據(jù)較好。繪制殘差圖：〔3〕對體重為57kg的女生的肺活量進行估計，翻開數(shù)據(jù)集直接在weight(體重)最下方輸入57,得到預測值。所以體重為57kg的女生的肺活量的估計值是3.3535。〔4〕并求出肺活量期望值的95%的可信區(qū)間。由上圖表可知，肺活量期望值的95%的區(qū)間為[2.6392,3.1658]。2、利用dst.fitness數(shù)據(jù)集，擬合多元線性方程〔以O(shè)XYGEN變量為因變量〕?！?〕擬合多元線性方程的方程。①用編程擬合多元線性回歸，編程如下：Procregdata=dst.fitness; ModelOXYGEN=ageweightruntimerstpulserunpulsemaxpulsegroup;run;②提交這一程序后可得到如下的輸出：這一輸出與簡單線性回歸的輸出時完全相仿的，只是進入回歸的自變量有7個。從參數(shù)的估計值容易得到擬合的回歸為：OXYGEN=-0.21043AGE–0.073700WEIGHT–2.680931RUNTIME–0.001829RSTPULSE–0.372725RUNPULSE+0.305364MAXPULSE-0.073612GROUP+101.824147〔2〕判斷方程是否存在共線性和異方差？如果方程存在共線性，怎樣用逐步回歸法進行變量篩選擬合回歸方程〔保存模式的顯著性水平為0.15〕。從參數(shù)估計檢驗局部可以看出，變量AGE、WEIGHT、RSTPULSE和GROUP的回歸系數(shù)，不能拒絕它們?yōu)?的原假設(shè)。在這里須小心地看待這些檢驗，因為它們都是在其他變量參加回歸的前提下進行顯著性檢驗的，完全可能因為自變量間存在較強的相關(guān)而掩蓋它們對回歸的奉獻，所以在剔除不顯著的回歸變量時必須逐個進行。①用編程方法判斷方程是否存在共線性Procregdata=dst.fitness;ModelOXYGEN=ageweightruntimerstpulserunpulsemaxpulsegroup/collinvif;run;②提交這一程序后可得到如下的輸出：由檢測結(jié)果可知，方差膨脹因子最大只有8.83077，并沒有超過10，這說明不存在共線性問題。③用編程的方法做逐步回歸法進行變量篩選擬合回歸方程。Procregdata=dst.fitness;ModelOXYGEN=ageweightruntimerstpulserunpulsemaxpulsegroup/selection=stepwise;run;④提交這一程序后可得到如下的輸出：在輸出報告中，提供了進入回歸變量逐次改變后回歸方差分析和擬合的信息。在報告的最后，利用逐步回歸法挑選，變量RUNTIME、AGE、RUNPULSE和MAXPULSE四個變量進入回歸。所有進入回歸的變量在0.15的水平下是顯著的?！?〕檢驗樣本中是否存在強影響點。①用編程的方法作回歸CookD診斷：Procregdata=dir.fitness;ModelOXYGEN=ageweightruntimerstpulserunpulsemaxpulsegroup/r;run;②提交這一程序后可得到如下的輸出：CookD統(tǒng)計量大于4/n說明是一個有影響的觀測。在這里，n=31，所以當CookD統(tǒng)計量大于4/31=0.129時，說明存在強影響點。由上圖可知，存在大于0.129的點，即檢驗樣本中存在強影響點。五、討論與結(jié)論花費了一定的時間完場此份實驗報告，證明課堂上所學的與實際應用之間存在著差距。在相關(guān)與回歸