嵌套logit模型在Stata軟件中的擬合

1、現(xiàn)代険防厭學(xué)年第 卷第期現(xiàn)代険防厭學(xué)年第 卷第期中圖分類(lèi)號(hào)文獻(xiàn)標(biāo)識(shí)碼哦】現(xiàn)代険防厭學(xué)年第 卷第期現(xiàn)代険防厭學(xué)年第 卷第期嵌套logi膜型在Stat漱件中的擬合王靜葉冬青摘要 日的 掌握嵌套 模世的原理及其住 軟件中的擬合過(guò)程。方法 運(yùn)用實(shí)例來(lái)介紹該模型的擬 介過(guò)程.結(jié)果 在 軟件中的擬介過(guò)程中需注點(diǎn)數(shù)據(jù)結(jié)構(gòu)的止呦性、構(gòu)建模型中的*層次及解秫變覽、檢驗(yàn) 假設(shè)扯后采用 模塊進(jìn)行嵌套 模型的擬合.結(jié)論嵌套 模別可以解決反應(yīng)變址的存?zhèn)€類(lèi)別之間 的相關(guān)性問(wèn)題正確分析因素対反應(yīng)變空的影響效應(yīng)。關(guān)鍵詞嵌奩 模烈軟件擬合現(xiàn)代険防厭學(xué)年第 卷第期現(xiàn)代険防厭學(xué)年第 卷第期個(gè)類(lèi)別例如反應(yīng)變就關(guān)性范用在 間相關(guān)性越小

2、。若* 均為楔型等價(jià)于條件嵌套 檢型又稱(chēng)廣義極值 模型 肖次由年 提岀 并由年 構(gòu)建 是姜項(xiàng) 模型的擴(kuò)展。原理嵌倉(cāng)模翌主?？梢越鉀Q反應(yīng)變就的外個(gè)類(lèi)別之間存在州關(guān)件的樹(shù)題 而這 點(diǎn)是多項(xiàng) ?？账荒苋萑痰?因?yàn)榻?xiàng)模型耍求并類(lèi)別Z間圧獨(dú)立不相關(guān)的我們稱(chēng)為假設(shè)即任何M類(lèi)別的槪率Zrt不隨其他類(lèi)別的存在與會(huì)而改變。該假設(shè)可以采用方法年來(lái)檢臉可先在條件模刑中擬合包括所有類(lèi)別的楔翌和去除某類(lèi)別后的模型再用方法比較乞個(gè)佔(zhàn)計(jì)系數(shù)之間的矩片 W為它們服從聯(lián)合x(chóng)分布 即進(jìn)行x檢驗(yàn) 若 值 則 不冊(cè)從 假設(shè) 提示不能采川婁項(xiàng)模型 耍用嵌套模熨.在嵌套 模型中將反應(yīng)變址的所右類(lèi)別分成若干少 基金項(xiàng)日年友安做省教（

3、fil A然科研及I!年度安徽省教件廳門(mén)然科研及目年度安嚴(yán)各島校行年教帥科研資助歩日作者簡(jiǎn)介 上靜女 偵士 副教授 研究方向iftlrtt學(xué)分為個(gè)類(lèi)別分為個(gè)亞組$包含包含、楔型第一層高層次為數(shù)亞細(xì)并亞紐內(nèi)包介至少?gòu)膶?zhuān)業(yè)1二將其作為新變fit包含、這樣就構(gòu)成了兩層嵌套第二層低層次為、分別對(duì)應(yīng)于低層次和高層次的無(wú)稈變蚩 模朋的槪章可以衣達(dá)為兩個(gè)簡(jiǎn)甲 模熨槪率的乘枳其中的條件概率3 E -3概率 的計(jì)并麗需定義高層次中毎組的包含值w 分為個(gè)亞俎即對(duì)應(yīng)個(gè) itL. 定義為E 0 u X* £ g A*其中入測(cè)粗了高層次中每個(gè)亞組內(nèi)未觀察因rzM的相之間入越k 說(shuō)明該亞細(xì)內(nèi)并類(lèi)別之則兩層嵌奩模

4、熨變成尬層嵌套模型?，F(xiàn)代険防厭學(xué)年第 卷第期在軟件中的擬合采用數(shù)抑;與衛(wèi)生統(tǒng)計(jì)作昌單位 安做悅科人學(xué)渝行病與衛(wèi)生統(tǒng)計(jì)學(xué)系 合肥現(xiàn)代険防區(qū)學(xué)年第卷第期文件來(lái)擬合。該資料中足 個(gè)家庭對(duì)為地類(lèi)飯店的選項(xiàng)還有一些變呈 家庭的收入、孩子數(shù) 、飯店的等級(jí)、平均每人毎頓飯的花費(fèi)、家到飯店的距肉前兩個(gè)變就在每個(gè)家庭內(nèi)是不變的而后三個(gè)變雄對(duì)于不同的 飯店込變化的。歆分析哪些因素影響家庭對(duì)飯店的選擇并計(jì) 舁出每個(gè)家庭迭擇每類(lèi)飯店的槪率人小從而作出預(yù)測(cè).數(shù)據(jù)結(jié)構(gòu)在擬介之詢(xún)注總數(shù)據(jù)結(jié)構(gòu)的正確性 例如 第一個(gè)家庭 選擇 中的項(xiàng)則數(shù)據(jù)躱入為為配比組此處為家庭序號(hào)為不選擇 為選擇人樣木蚩資料 這樣的數(shù)據(jù)準(zhǔn)備匸作比較繁瑣。構(gòu)

5、建嵌套懊刊中的并層次及其解秤變毘擬合嵌套模型時(shí)先將個(gè)飯店分為 個(gè)亞組毎 個(gè)亞組包含若干個(gè)選項(xiàng) 見(jiàn) 中的閘述 所以耍采用命令來(lái)產(chǎn)生一個(gè)能識(shí)別第一層選項(xiàng)的新分類(lèi)變盤(pán) 見(jiàn) 數(shù)據(jù)結(jié)構(gòu)再分別確定第一層次和第二層次的解釋變址 、和為低層次的解釋變st而島層次的解釋變新產(chǎn)生 以為參照組 產(chǎn)生分別與、 之間的交互項(xiàng)作為島丿/次的解禪變呈現(xiàn)代険防區(qū)學(xué)年第卷第期假設(shè)的檢驗(yàn)在條件 模童中 采用 模塊 以 為反應(yīng)變 就指定配比組變址設(shè)定解科變覽個(gè)運(yùn)行結(jié)果見(jiàn)表.將哄保仔到變1S中。再在模取菇礎(chǔ)上加上運(yùn)行結(jié)果見(jiàn)表址后將它打 進(jìn)行檢臉結(jié)果見(jiàn)農(nóng)A若有 配比組則數(shù)據(jù)共仃檢臉麻ffi接近不能認(rèn)為滿(mǎn)足 假設(shè)可以采用現(xiàn)代険防區(qū)學(xué)年第

6、卷第期現(xiàn)代険防區(qū)學(xué)年第卷第期> I5個(gè)解釋變的條件logi（權(quán)釵（織型1 ）結(jié)果口變系數(shù)siZ值尸值95% CI«&A406166210.067M2-0.910.363-0.1946496*0.0713254046590010.2832041Q590.558-039138.07237-412443960.0993056-1460014-O.4W4254 -0.04916H7收人啖申1-G07375060.01774447】600001(*529-).0389721孩子數(shù)類(lèi)或1(14105386021370511.920.055-O.OM3157 - 08293928(t

7、t -«5 955324表2 7個(gè)解釋變J|的條件logits型（樓型2）結(jié)果系數(shù)Z值PiA95%C7收入-0.1367799O.O35M793«20.000-0.2070404-).0665193飯店編級(jí)0.30666220J41829I2600310 O2S6X20 5IU642E*-O.1977SO5a(M71653-4.190.000).2901927-0.1053082收人啖型1-0.0390183aoowiH7150000-410574455 - -0.0205911收人«« J0.0407053(100MMOS5.0600000.02494

8、62 -0.0564644孩子敗黃尿1-0.23987570.1063674-2.260.024-0.44852-0.03H9<M孩子效啖蟲(chóng)3-0389862dll 43797-3.400001-0.6135662-0.1652061注】對(duì)數(shù)佢煞值二-4M.90834現(xiàn)代険防區(qū)學(xué)年第卷第期現(xiàn)代険防區(qū)學(xué)年第卷第期械倉(cāng)模型的擬介采用 棋塊 以 為反應(yīng)變割 折定配比組變啟 分別設(shè)定第厘和第二夷的解釋變空 運(yùn)行結(jié)宋見(jiàn)農(nóng)。嵌套 ?？沾驐l件 模型的似然比檢膾更加說(shuō)期r該數(shù)據(jù)適合用磁套 模型。嵌套 樓型分析后發(fā)現(xiàn)平均毎人毎頓飯的花費(fèi)、 家到飯店的灰離、毎個(gè)家庭的收入、孩了數(shù)妙響家庭對(duì)飯店的 選擇。進(jìn)-

9、步采用 命令計(jì)舜毎個(gè)家庭選擇毎種飯店的預(yù) 測(cè)槪率計(jì)克公式見(jiàn)前賂 從而進(jìn)行冋代預(yù)測(cè) 若某家庭 選擇 類(lèi)飯店的険測(cè)概率中 類(lèi)所對(duì)應(yīng)的碩測(cè)概率垠人 則 該家庭判為選抒 類(lèi)飯店在 個(gè)家庭的 次選擇中現(xiàn)代険防區(qū)學(xué)年第卷第期現(xiàn)代険防國(guó)年第 卷第期預(yù)測(cè)結(jié)果與實(shí)處選擇結(jié)果相比次選擇是止確的UP 1E確率為X -T*3 型1和櫃型2之間基異的hausman檢區(qū)paHset ( b)(B)左值(l»B)壟值的標(biāo)冷觀收人-0.0616621-0“477990.07511780.0576092飯店等級(jí)0.16590010.3066622-0.1407621O245I3OKfem-0 244396-0.197

10、7505-0.04664S60.0R76173收入矣血1-0 0737506-0.0390183-0.03473230.015049出子散啖?H10.41053&6-0.23987570.65041430.1853533注】檢笛后p值接近aos.不能認(rèn)為摘足ila假設(shè).可以采用hcit哀4 ««logit m型的忍合結(jié)果吐ZfflLP也95%a飯店花費(fèi)-0.09443520.OM02-2.78aoo6-0.1611131 -O.O277572飯E尊級(jí)0.17937590.1268951.41ai57-O.O693338 - O.42HO855距宴-0.1745797

11、0.0435352-4.030.000-OJ595152 *-0.0896443類(lèi)取收人類(lèi)型14.02875020.0 16242-2.47aoi3-O.O515332-O.OO59672收人晏5! 3004583730.00891095.14QOOO0.02S3722 - 0.0633024-007(M1640.BM359-0.510614-03437058 0 2028729孩子數(shù)類(lèi)別3-03626W10.1171277-3.100002-0.5922041-0.1330721類(lèi)型次SH5.71575813321711450.0141.143415 1Q2881庚8! 21.7212221.1520021.490135-O.53666OR-3 979105渙宙3I.4665M0.416907S3320.0000.34642-2J837M桂似能比檢段.護(hù)值"0快學(xué)領(lǐng)域中普遍使用的原因。參考文獻(xiàn)討論木文所閘述的嵌奩 模型雖然可以耳成但為了與爹水平模醴區(qū)別開(kāi) 這里的被翻譯成 占r 而且第 層為址島肚 第二層為次島層 以次類(lèi)推 在姜水平楔 乜中水平”也叫 但水平 為滋低水平 水平 為較島水平以次類(lèi)推。筆者對(duì)該數(shù)據(jù)進(jìn)行條件預(yù)測(cè) 發(fā)現(xiàn)其険測(cè)正確率為模世分析忘也進(jìn)行了回代離于嵌套 模型的預(yù)測(cè)正確率是因?yàn)闂l件模型的- 對(duì)數(shù)似然函數(shù)現(xiàn)代険防國(guó)年第 卷第