作者coldcurrent (under the net)
看板Statistics
標題[程式] STATA 在probit模型下跑Chow test
時間Fri Jul 29 15:19:44 2011
[軟體程式類別]:STATA
[程式問題]:跑完probit後的統計檢定
[軟體熟悉度]:中
[問題敘述]:
各位好,我的模型有將男性和女性樣本分組,但是跑的模型相同。
而希望檢定不同性別的迴歸結果是否有差異。
線性模型有Chow test做類似的檢定:
H0: 兩組資料模型及參數完全一致
Ha: 兩組資料模型相同但係數值不同
在非線性模型的檢定中,目前打算用likelihood ratio test,
也就是用STATA內建的lrtest。
然而LR test比較常用在檢定變數增加或減少時,對模型的解釋力是否改變。
需要設定模型內幾個變數為零,成為受限制模型後,
再用非受限和受限制模型的log-likelihood值進行檢定。
但我的檢定並非如此。
搜尋STATA的說明後,LR test下有Chow-type LR test,
解釋如下:
Chow-type tests are appropriate for hypotheses that specify that all
coefficients of a model do not vary between disjointed subsets of the data.
似乎是說這類型檢定的虛無假設為:
兩個不重疊(disjointed)的子集(subsets)間,所有的參數沒有顯著差異。
想請問做過類似檢定的先進們,我的理解是否有誤?
[程式範例]:
我的模型是這樣:
1.跑總樣本的probit
probit y x1 x2 x3
est sto all
2.跑女性樣本的probit
probit y x1 x2 x3 if gen == 0
est sto female
3.跑男性樣本的probit
probit y x1 x2 x3 if gen == 1
est sto male
4.進行LR test
此處是參考STATA軟體內-help lrtest-的範例
lrtest (all) (female male)
而結果為:
Likelihood-ratio test LR chi2(15) = 13.21
Prob > chi2 = 0.5857
Assumption: (all) nested in (male, female)
在ptt上我簡化了原本模型,所以LR chi2的數值不一樣。
想問,0.5857應該落入了拒絕域,那是拒絕"do not vary"的虛無假設嗎?
(代表兩個群體的迴歸結果有差異?)
這樣的做法是否達成我原本要檢定的目標?
謝謝板上的各位!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.24.15.58
※ 編輯: coldcurrent 來自: 114.24.15.58 (07/29 15:21)