基礎(chǔ)醫(yī)學異方差

基礎(chǔ)醫(yī)學異方差第二篇

放寬經(jīng)典模型的假定Dr.Ouyang2學習思路概念成因 后果檢驗消除方法Dr.Ouyang3第五章

異方差Dr.Ouyang4§1什么是異方差Dr.Ouyang5(A)(B)密度儲蓄Y收入X密度儲蓄Y收入X異方差的圖形表示Dr.Ouyang6（A)與(B)的比較：相同點：收入增加，儲蓄平均來說也增加。不同點（A)儲蓄的方差在所有的收入水平上保持不變。（B)儲蓄的方差隨收入的增加而增加。解釋：隨收入增長，人們有更多的備用收入，從而如何支配他們的收入有更大的選擇范圍。Dr.Ouyang7例2：邊錯邊改學習模型人們在學習的過程中，其行為誤差隨時間而減少。例如，在給定的一段時間里，大字出錯個數(shù)與用于打字練習的小時數(shù)的關(guān)系。隨著打字練習小時數(shù)的增加，不僅平均打錯字數(shù)，而且打錯個數(shù)的方差都有所下降。（圖略）注：一般橫截面數(shù)據(jù)容易出現(xiàn)異方差。Dr.Ouyang8以某一行業(yè)的企業(yè)為樣本建立企業(yè)生產(chǎn)函數(shù)模型Yi=Ai1

Ki2

Li3ei

被解釋變量：產(chǎn)出量Y解釋變量：資本K、勞動L、技術(shù)A，

那么：每個企業(yè)所處的外部環(huán)境對產(chǎn)出量的影響被包含在隨機誤差項中。每個企業(yè)所處的外部環(huán)境對產(chǎn)出量的影響程度不同，造成了隨機誤差項的異方差性。這時，隨機誤差項的方差并不隨某一個解釋變量觀測值的變化而呈規(guī)律性變化，呈現(xiàn)復雜型。例3:企業(yè)生產(chǎn)函數(shù)Dr.Ouyang9§2出現(xiàn)異方差時的OLS估計Dr.Ouyang10一、參數(shù)OLS估計的方差增大參數(shù)OLS估計仍然是線性無偏的。Dr.Ouyang11為什么β2的估計值的方差可以直接用這個等式求到？Dr.Ouyang12二、t檢驗失效這是因為直接應用OLS所計算出的參數(shù)的不是最優(yōu)的（即不具備方差最小性），從而導致t檢驗值變小。Dr.Ouyang13§3異方差的檢驗檢驗思路：由于異方差性就是相對于不同的解釋變量觀測值，隨機誤差項具有不同的方差。那么：檢驗異方差性，也就是檢驗隨機誤差項的方差與解釋變量觀測值之間的相關(guān)性及其相關(guān)的“形式”。問題在于用什么來表示隨機誤差項的方差思路一般是，用OLS方法來估計模型，求得隨機擾動項的估計量（注意，該估計量是不嚴格的），即“近似估計量”，再用其平方來表示隨機誤差項的方差。Dr.Ouyang14一、非正式方法1、實際問題：如截面數(shù)據(jù)2、圖解法：......................Dr.Ouyang15.......................................圖解法Dr.Ouyang16Dr.Ouyang172、格萊澤（Glejser)檢驗Dr.Ouyang183、戈德菲爾德—匡特（Goldfeld-Quandt)檢驗（1）

（2）（3）檢驗統(tǒng)計量：（4）判別：

注：“大在分子，小在分母”?？偨Y(jié)：至今沒有很好的檢驗方法！Dr.Ouyang19C個IIIX1Xn(n-c)/2個(n-c)/2個Dr.Ouyang20布魯士-帕根檢驗(BP檢驗)

Breusch-Pagantest(BPtest)假設u2與x的線性組合之間存在某種函數(shù)關(guān)系：E(u2|x1,…,xk)=f(x)=f(d0+d1x1+…+dkxk)常見：f(z)=zu2=d0+d1x1+…+dkxk+e此時，異方差檢驗等價于檢驗：H0:d1=…=dk=0問題：u2isunobservable一個自然的想法：用^u2代替u2

^u2=d0+d1x1+…+dkxk+eFtest/LMtest：H0：overallsignificanceDr.Ouyang21布魯士-帕根檢驗(BP檢驗)

Breusch-Pagantest(BPtest)(i).用OLS估計多元回歸模型：y=b0+b1x1+b2x2+...bkxk+u得到^ui，計算^ui2.(ii).用^ui2對所有解釋變量進行回歸：^ui2=d0+d1xi1+…+dkxik+ei得到擬合優(yōu)度R2u?2.(iii).對回歸方程(ii)進行總體顯著性檢驗：H0:d1=…=dk=0BP檢驗[Breusch-Pagantestforheteroskedasticity

(BPtest)]:LMversionDr.Ouyang22布魯士-帕根檢驗(BP檢驗)

Breusch-Pagantest(BPtest)如果懷疑異方差僅與某些特定解釋變量xj-xj+q有關(guān)，則在第(ii)回歸時,只需將^ui2對xj-xj+q，并進行F/LMtest即可.Dr.Ouyang23懷特檢驗

TheWhiteTest但是，異方差可能是來源于：x高次項(2次項、交互項)之間可能存在函數(shù)關(guān)系：假定E(u2|x1,…,xk)=f(x1,…,xk)是x1,…,xk的二次函數(shù),即所有解釋變量的一次項、平方項和交互項的線性組合,則(當k=2時)：^u2=d0+d1x1+d2x2+d3x12+d4x22+d5x1x2+eFtest/LMtest:H0:d1=0,d2=0,d3=0,d4=0,d5=0懷特異方差檢驗(TheWhitetestforheteroskedasticity):LMversion問題：當k很大時，上述回歸方程中的解釋變量迅速增加。Dr.Ouyang24懷特檢驗的一種變形

AlternateformoftheWhitetest回顧：?i=^b0+^b1xi1+^b2xi2+...+^bkxik

?i2是x1…xk的平方項和交互項的一個線性組合在一定條件下，可以用

?i

和

?i2的某種線性組合，來代替x1…xk的一次項、平方項和交互項線性組合，故可將^u2對?i,?i2進行回歸:于是：^u2=a0+a1?i+a2

?i2+eFtest/LMtest:H0:a1=0,a2=0Note：onlytestingfor2restrictionsnowDr.Ouyang25§4、異方差的修正Dr.Ouyang26Dr.Ouyang27

Dr.Ouyang28Dr.Ouyang29當未知X.....0.....X0...Dr.Ouyang30Dr.Ouyang31MoreonWLSWLSisgreatifweknowwhatVar(ui|xi)lookslike 如果我們知道Var(ui|xi)的形式，WLS很棒Inmostcases,won’tknowformofheteroskedasticity 在大多數(shù)情況下，我們并不清楚異方差的形式Dr.Ouyang32FeasibleGLS

可行GLSInthiscase,youneedtoestimateh(xi) 此時，你需要估計h(xi)Typically,westartwiththeassumptionofafairlyflexiblemodel,suchas 我們可以從一個非常靈活的方程形式入手 Var(u|x)=s2exp(d0+d1x1+…+dkxk)Sincewedon’tknowthed,mustbeestimated 由于d未知，我們必須對它進行估計。為什么不用前面BP檢驗中的水平線性模型？Dr.Ouyang33FeasibleGLS(continued)

可行GLSOurassumptionimpliesthat 我們的假定意味著 u2=s2exp(d0+d1x1+…+dkxk)v, whereE(v|x)=1.ln(u2)=a0

+d1x1+…+dkxk+eWhereE(e)=1andeisindependentofx 其中E(e)=1且e獨立于xNow,weknowthat?isanestimateofu,sowecanestimatethisbyOLS

現(xiàn)在，我們知道?是u的一個估計，所以我們可以通過OLS對其進行估計。Dr.Ouyang34FeasibleGLS(continued)

可行GLSNow,anestimateofhisobtainedas?=exp(?),andtheinverseofthisisourweight 對h的估計可以通過?=exp(?)得到，其倒數(shù)為我們的權(quán)重So,whatdidwedo? 那么，我們做了什么呢？RuntheoriginalOLSmodel,savetheresiduals,?,squarethemandtakethelog 對原方程做OLS回歸，保存殘差?，平方之，并取自然對數(shù)Regressln(?2)onalloftheindependentvariablesandgetthefittedvalues,?

將ln(?2)對全部解釋變量回歸，得到預測值?DoWLSusing1/exp(?)astheweight 將1/exp(?)作為權(quán)重，做WLSDr.Ouyang35Dr.Ouyang36三、具體例子1988年美國的研究開發(fā)與支出18個工業(yè)行業(yè)，RD～Sale1、OLS:Se=(990.99)(0.0083)t=(0.1948)(3.8434)R2=(0.4783)t值在0.002水平上是統(tǒng)計顯著的Dr.Ouyang372、Park檢驗t=(0.8572)(1.1626)R2=0.0779無法拒絕同方差性Dr.Ouyang383、Glejser檢驗t=(0.8525)(2.0931)R2=0.2150t=(-0.5032)(2.3704)R2=0.2599t=(3.7601)(-1.6175)R2=0.1405①②③