作者: StatGuest (統計過客) 看板: P_yhliu
標題: Re: 如何比較兩個迴歸模型的係數是否有差
時間: Sun May 13 08:39:56 2007
※ 引述《yhliu (老怪物)》之銘言：
> ※ 引述《Asterid.bbs@ptt.cc (asterid)》之銘言：
> > 想請問各位，倘若我用STATA跑兩個模型(logistic regression)
^^^^^^^^^^^^^^^^^^^
老眼昏花沒看清楚!
> > 第一個模型放的是main effect,第二個模型裡則加了x4進去，
> > 那該如何比較這兩個模型裡，
> > 或者將同一個model拆成兩組不同的樣本去跑(因此樣本數會不同)...
> > MODEL3: male: y= x1+x2+x3
> > MODEL4: female: y= x1+x2+x3
> > ..比較x1對y的影響是否會隨著性別的差異而有差....
> 令兩樣本之係數估計結果分別是 b, b', 則
> t = (b-b')/√[(se(b))^2+(se(b'))^2]
^^^^^^^^^^^^^^^^^^^
(ase(b))^2+(ase(b'))^2
> 可用於檢定兩樣本同模型對應係數之差異是否顯著. t 之
> 自由度可取兩樣本迴歸模型殘差自由度中較小的.
> 一般就相同目的, 在不同假設下, 可能考慮合併樣本模型
> Y = α + γG + β1*x1 + β2*x2 + β3*x3
> + γ1*G*x1 + ε
logistic reg. 模型應為
logit(p) = α + γG + β1*x1 + β2*x2 + β3*x3
+ γ1*G*x1
> 其中 ε 是隨機誤差項. 以一般統計軟體之模型表示法,
> 即
> y = G + x1 + x2 + x3 + G*x1
> 檢定 γ1, 即 G*x1 之係數是否為 0.
> 兩種方法結果不一定相同, 因後者假設沒有 G*x2, G*x3
> 等項, 且假設兩樣本具有相同誤差項變異數.
就標準 logistic 模型而言, 此與原方法(分開跑模型)的
差異只在於沒有 G*x2, G*x3 兩項.
> > 這問題我想蠻久了，一時找不到答案，有點急..煩請各位高手協助我，感謝^^
--
yhliu 在 07/05/13 8:46:16 從 163.15.188.87 修改這篇文章