公衛醫師醫學統計學輔導:直線回歸

一、直線回歸方程的意義 計算出相關係數後,如果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)