姜啟源之統(tǒng)計(jì)回歸模型-資料下載頁(yè)
2025-04-29 05:35本頁(yè)面
【正文】 年齡的關(guān)系 ,并進(jìn)行統(tǒng)計(jì)預(yù)測(cè) . 序號(hào) 年齡 冠心病 序號(hào) 年齡 冠心病 序號(hào) 年齡 冠心病 序號(hào) 年齡 冠心病 1 20 0 26 35 0 51 44 1 76 55 1 … … … … … … … … … … … … 25 34 0 50 44 0 75 55 1 100 69 1 分析與假設(shè) 被觀察者獨(dú)立選取 x~被觀察者 年齡 , Y~患病情況 (Y=1~患病 , Y=0~不患病 ) 20 25 30 35 40 45 50 55 60 65 7000 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 91A g eCHD20 25 30 35 40 45 50 55 60 650 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 8A g eProportion of CHD年齡段 段中點(diǎn) 人數(shù) 患病人數(shù) 患病比例 2029 10 1 3034 32 15 2 … … … … … 6069 10 8 合計(jì) 100 43 按年齡段 分組 統(tǒng)計(jì)患病人數(shù)及 比例 患病比例隨年齡增大而遞增 ,是介于 0與 1之間的 S型曲線 . 無(wú)法建立前面那樣的回歸模型 ,需要對(duì)數(shù)據(jù)進(jìn)行預(yù)處理 . 患病比例 y是 年齡段中點(diǎn) x時(shí) Y的平均值 分析與假設(shè) Y取值 0, 1 。 y 取值 [0, 1] ????? ????? 332210 xxxy用普通方法建立回歸方程 20 25 30 35 40 45 50 55 60 65 7000 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 91A g eProportion of CHD? y取值不一定在 [0,1]中 . ? 誤差項(xiàng) ε只能取值 0,1, 不具有正態(tài)性 , 且具有異方差性 . 違反普通回歸分析的前提條件 ! 當(dāng)因變量 Y為一個(gè) 二分類 (或多分類 )變量時(shí)，需要用到 新的回歸模型 . )|( xYEy ?Y的 條件期望 Logit 模型 )|1()( xYPx ???)|( xYEy ?Y的 (條件 )期望 π(x)~年齡 x的患病概率 (患病比例 y) ))(1)(()|( xxxYD ?? ??方差 π(x) ~ S型曲線 , 取值 [0,1] xxeex10101)( ????? ???? xxx10))(1)(ln ( ???? ???Logit模型 (Logistic回歸模型 ) ))(1 )(ln ())(( xxxL o g i t ??? ??π(x)的變換 ),( ????取值 連接函數(shù)， Logistic模型 反函數(shù) Logit 模型 數(shù)據(jù)預(yù)處理 : 將年齡分成 k(=8)組 . xi~第 i組年齡 , ni~被觀察人數(shù) , mi~患病人數(shù) , i=1,…, k iii nm /??患病概率 β0,β1~回歸系數(shù) ),( iinB ?設(shè) mi服從二項(xiàng)分布 回歸系數(shù)可用極大似然法估計(jì)得到 . iiii xL og i t 10)1l n ()( ????? ????Logit 模型 模型求解 Logit模型可用 MATLAB命令 glmfit求解 b = glmfit(x, y, ’distr’, ’link’) [b,dev,stats] = glmfit(x, y, ’distr’, ’link’) x~自變量數(shù)據(jù)矩陣 (第 1列自動(dòng)添加列向量 1). y~因變量數(shù)據(jù)向量 (對(duì) distr =binomial, y可取矩陣 : 第 1列為 “成功”次數(shù) , 第 2列為觀察次數(shù) ). ’distr’ ~估計(jì)系數(shù)所用分布 (’binomial’,’poisson’ 等 ),缺省時(shí)為 ’ normal’ . ’link’ ~’logit’,’probit’ 等 (缺省時(shí)為’ logit’). b~回歸系數(shù)的估計(jì)值 , dev~擬合偏差 , stats~統(tǒng)計(jì)指標(biāo) 模型求解 參數(shù) 估計(jì)值 標(biāo)準(zhǔn)差 β0 β1 [yhat, dylo, dyhi] = glmval(b, x, 39。logit39。) 年齡段 年齡 x 患病比例 ( 實(shí)際值 ) 患病比例 (預(yù)測(cè)值 y) 置信區(qū)間 2029 [, ] … … … … … 6069 [, ] 編程計(jì)算 20 25 30 35 40 45 50 55 60 6500 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 9A g eProportion of CHD回歸曲線 自變量為 x時(shí) y的預(yù)測(cè)值 yhat及置信度為 95%的置信區(qū)間 擬合偏差 模型評(píng)價(jià)與結(jié)果分析 xxxxL o g i t 10))(1 )(l n ())(( ????? ????2210 xx ??? ??? Logit模型是否需要引入 x2項(xiàng)？ 用似然比統(tǒng)計(jì)量計(jì)算： pval = 1 chi2cdf(devdev2,1) = ? 選用 Probit模型 (另一種廣義線性模型 )結(jié)果如何？ )()( 10 xx ??? ??? xxxo b i t 101 ))(())((Pr ???? ???? ? Ф是正態(tài)概率分布函數(shù) (S型曲線 ) 模型中 引入 x2項(xiàng)不能顯著提高擬合程度 . glmfit中需將 logit改為 probit 20 25 30 35 40 45 50 55 60 6500 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 9A g eProportion of CHD D a t aL o g i t m o d e lP r o b i t m o d e l年齡段 年齡 x 患病比例( 實(shí)際值 ) 預(yù)測(cè)值 1（ Logit） 預(yù)測(cè)值 2（ Probit） 2029 … … … … … 6069 兩個(gè)模型的擬合程度不相上下 . 模型評(píng)價(jià)與結(jié)果分析 Probit模型 求解 參數(shù) 估計(jì)值 標(biāo)準(zhǔn)差 β0 β1 擬合偏差 xxxxxL o g i t 1 0 5 3 8 ))(?1 )(?l n ())(?( 10 ??????? ?????β1的直觀解釋 Odds~事件發(fā)生 (患病 )概率與不發(fā)生 (不患病 )概率之比 . 模型評(píng)價(jià)與結(jié)果分析 )( 10)(1)()( xexxxO d d s ???? ????年齡 x的人患與不患病概率之比 年齡增加 1歲的 Odds比 (發(fā)生比率 ) 11010 )1()()1( ?????eeexO d d sxO d d s xx??? ???年齡增加 1歲 Odds比的對(duì)數(shù) ))( )1(ln (1 xO d d sxO d d s ???年齡增加 k歲后的 Odds )()( 1 xO d d sekxO d d s k ???模型評(píng)價(jià)與結(jié)果分析 0 5 0 )20(? ??xxxxxL o g i t 1 0 5 3 8 ))(?1 )(?l n ())(?( 10 ??????? ?????20歲 的青年人患冠心病的概率 發(fā)生比 (患與不患冠心病的概率之比） 0 5 9 )20( ?O d d s10年后 30歲 人的發(fā)生比 1 6 9 5 9 )30( 110 ??? ? ?eO d d s60歲 時(shí) )60( ?O d d s)(? * ?x? 48* ?x48歲時(shí)患冠心病的概率會(huì)大于不患冠心病的概率 . 年齡增加 1歲患病概率的變化很小 . Logit回歸模型 ?)( ??? xx?是 20歲的 倍 140 ?? ?e? 因變量是定性變量的回歸分析作為一種有效的數(shù)據(jù)處理方法已被廣泛應(yīng)用，尤其在醫(yī)學(xué)、社會(huì)調(diào)查、生物信息處理等領(lǐng)域 . ??????miii xxxxL og i t10))(1)(l n ())(( ?????mxx ,1 ?自變量 可以是定量變量或定性變量 ? 可以用 逐步回歸 方法建立多元 Logit模型和 Probit模型 , 逐個(gè)地加入自變量 (包括某個(gè)自變量的高次項(xiàng)及某些自變量的交叉變量 ), 并且實(shí)時(shí)地進(jìn)行模型比較檢驗(yàn) , 選擇與數(shù)據(jù)擬合較好的模型 . 模型評(píng)述 多元 Logit模型
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1