作者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)