2013-06-13 08:45閱讀：

大家在計算嶺回歸的時候是不是都遇到了這樣的問題：最后的k值都定了，標準線性回歸系數(shù)和非標準線性回歸系數(shù)都算出來了，但是就是沒有t檢驗值和顯著性概率？我也很困擾，而且困擾了我好幾天，這幾天把整個論壇都找遍了也沒找到想要的信息，二是看到一位腦殘的同志狂發(fā)SPSS 18.0可以直接計算嶺回歸，但是我用何曉群的《應用回歸分析》里面的數(shù)值代進去是不對，這對于我們SPSS的菜鳥來說是個悲劇�。。�！
好在狗屎運，在百度找了好久，終于功夫不負有心人，找到了，現(xiàn)在我把原文章和文章出處貼上來，至于大家會問，這段代碼可行嗎？放心，我試過了！感謝這位好老師！淚奔！我知道板上有很多更我一樣的菜鳥，每次遇到這類問題都是束手無策，有的時候還是要靠自己啊，誒！
具體做法就是在你的SPSS的安裝文件中找到ridge regression.sps這個文件，雙擊打開，看到里面很多代碼，是的。然后找到和下面代碼一樣的那段代碼，將褐色的那部分加進去。保存（ms不能直接保存），那就另存到非安裝文件夾里面，比如桌面，然后拷回來。這就是可以計算t和p值的ridge regression.sps插件了。你可以用下面的進行計算了
INCLUDE'ridge regression文件的路徑\ridge regression.sps'.
ridgereg enter=自變量
/dep=應變量
/start=
/stop=
/inc=
我自己安裝的SPSS不能給出嶺回歸的顯著性檢驗，也不知道是什么原因。代06級統(tǒng)計學專業(yè)的實用回歸分析時就遇到了這個問題。SPSS的教程關于這個部分并未作詳細描述，在網上搜索也沒找到解決方案。所以當時就把這個問題忽略了。今年又代07級統(tǒng)計學，再次遭遇這個問題，經過多方嘗試與摸索，終于找到一個簡單可行的解決方案！先說發(fā)現(xiàn)這個方法的過程。
最開始想的是通過syntax窗口編程，在原有嶺回歸語句上添加上給出顯著性檢驗的語句。對SPSS程序語法不是太熟，但是SPSS有個特點是，對于窗口化操作的命令，最后不選“OK”而選“Paste”的話，就能夠進入syntax窗口，看到具體的程序�；谶@點，可以仿造相應程序編一個顯著性檢驗的語句。但是，這個想法失敗了。

看來要對嶺回歸的程序包添加相應語句才行。打開Ridge regression.sps，里面是密密麻麻關于嶺回歸的程序。還好，每個主要步驟都作了說明。查找一下，找到了嶺回歸系數(shù)估計的部分，果然沒有顯著性檢驗的語句。
還好，通過原有語句不難得到檢驗統(tǒng)計量的值。但是僅僅給出這個值，，對于做檢驗而言，很不方便。跟一般回歸一樣，要給出顯著性P值才好。這個關系不難根據(jù)P值意義得到，關鍵是相應的程序語句怎么寫。對照系數(shù)估計上面方差分析的部分（方差分析部分給出了F檢驗顯著性P值），嘗試著寫程序，終于成功！保存之后，再做嶺回歸就能給出顯著性檢驗的P值了！
沒用過別的SPSS版本，有的版本好像沒這個問題。提出這個方法，希望可以幫到有同樣問題的人。
附：在Ridge regression.sps中添加的語句(紅色傾斜的部分，就是添加的語句，當然，為了給出相應結果，原有語句作了少許改動。)
*---------------------------------------------------------------------------.
* Calculate raw coefficients from standardized ones, compute standard errors
* of coefficients, and an intercept term with standard error. Then print
* out similar to REGRESSION output.
*---------------------------------------------------------------------------（從這里開始是給出系數(shù)估計）
. compute beta={b;0}.
. compute b= ( b &/ std ) * sy.
. compute intercpt=ybar-t(b)*t(xmean).
. compute b={b;intercpt}.
. compute xpx=(sse/(sst*(n-nv-1)))*inv(xpx+(k &* ident(nv,nv)))*xpx*
inv(xpx+(k &* ident(nv,nv))).
. compute xpx=(sy*sy)*(mdiag(1 &/ std)*xpx*mdiag(1 &/ std)).
. compute seb=sqrt(diag(xpx)).
. compute seb0=sqrt( (sse)/(n*(n-nv-1)) + xmean*xpx*t(xmean)).
. compute seb={seb;seb0}.
. compute rnms={varname,'Constant'}.
. compute ratio=b &/ seb.
. compute ppp=2*(1-tcdf(abs(ratio),n-nv-1)).
. compute bvec={b,seb,beta,ratio,ppp}.
. print bvec/title='--------------Variables in the Equation----------------'
/rnames=rnms /clabels='B' 'SE(B)' 'Beta' 'T' 'sig'.
. print /space=newpage.
end if.

分享

N