公衛醫師醫學統計學輔導：直線回歸
發布時間：2011-08-26 共2頁
計算出相關系數后，如果r顯著，且又需要進一步了解兩變量中一個變量依另一個變量而變動的規律時，則可進行回歸分析。
“回歸”是個借用已久因而相沿成習的名稱。若某一變量(Y)隨另一變量(X)的變動而變動，則稱X為自變量，Y為應變量。這種關系在數學上被稱為Y是X的函數，但在醫學領域里，自變量與應變量的關系和數學上的函數關系有所不同。例如成年人年齡和血壓的關系，通過大量調查，看出平均收縮壓隨年齡的增長而增高，并且呈直線趨，但各點并非恰好都在直線上。為強調這一區別，統計上稱這是血壓在年齡上的回歸。
直線回歸分析的任務就是建立一個描述應變量依自變量而變化的直線方程，并要求各點與該直線縱向距離的平方和為最小。按這個要求計算回歸方程的方法稱為最小平方法或最小二乘法。所建立的方程是一個二元一次方程式，其標準形式是：
=a+bX(9.5)
式(9.4)為由X推算得來的Y值，即Y的估計值：a稱為截距，它是當X=0時的值，即回歸直線與縱軸的交點：b稱為 回歸系數，它是回歸直線的斜率，其含意是當X每增加一個單位時，相應增(或減)b個單位。當a與b求得后，直線回歸方程就確定了。
二、直線回歸方程的計算法
仍以表9.1資料為例，根據前面的相關分析以及醫學上有關凝血的機理，可知凝血時間依凝血酶濃度而異，且有密切的關系。因此可進一步作由凝血酶濃度(X)推算凝血時間(Y)的回歸方程。求直線回歸方程的步驟如下：
1.列回歸計算表(見表9.1)，計算∑X、∑Y、∑X2、∑Y2、∑XY。
2.計算X、Y、∑(X-X)2、∑(X-X)(Y-Y)
X=∑X/n=15.1/15=1.01
Y=∑Y/n=222/15=14.80
∑(X-X)2=∑X2-(∑X)2/n=0.2093
∑(X-X)(Y-Y)=∑XY-∑X·∑Y/n=-1.7800
3.計算回歸系數b和截距a。b和a兩值計算公式均是根據最小二乘法的原理推算出來的，其公式如下：
(9.5)
a=Y-bX (9.6)
本例b=-1.7800/0.2093=-8.5045
a=14.80-(-8.5045)(1.01)=23.3895
4.列出回歸方程，繪制回歸直線，將求得的b和a的值代入到式(9.4)，即得所求的回歸方程：
=23.3895-8.504X
在凝血酶濃度的實測范圍內，即X=0.8到X=1.2之間，任選兩個X值(一般選相距較遠且直角坐標系上容易讀出者)，代入此回歸方程，即得相應的兩個值。例如：
取 X1=0.8，則1=23.3895-8.5045×0.8=16.59,
X2=1.2 則2=23.3895-8.5045×1.2=13.18。
連接(0.8、16.59)和(1.2 、13.18)兩點所得直線，即為由凝血酶濃度推算凝血時間的回歸直線( 見圖9.9)。須注意回歸直線必通過(χ，y )點，并穿過觀察點群，直線上下各有一些點散布著，否則計算有誤。
三、直線回歸方程的假設檢驗
(一)樣本回歸系數的假設檢驗
根據例9.1資料求得的是樣本回歸系數b，有抽樣誤差的，需作假設檢驗，檢驗其是否是從回歸系數為0的假設總體(即β=0)中隨機抽得的，也就是檢驗b與0的差別有無顯著性。如果差別有顯著性，可認為X與Y間有直線回歸存在。
樣本回歸系數的假設檢驗亦用t檢驗。
H0：β=0 即Y的變化與X無關；
H1：β≠0。
計算公式為：
(9.7)
分母Sb是樣本回歸系數b的標準誤，計算公式為：
(9.8)
分子Sy.x為各觀察值Y距回歸線的標準差，即當X的影響被扣去以后Y方面的變異，可按下式計算：
(9.9)
