資源描述:
《離散選擇模型舉例122》由會員上傳分享,免費在線閱讀,更多相關內容在行業(yè)資料-天天文庫。
1、一.二元離散選擇模型1.二元響應模型(Binaryresponsemodel)我們往往關心響應概率,其中x表示各種影響因素(各種解釋變量,包括虛擬變量)。根據(jù)不同的函數(shù)形式可以分為下面三類模型:線性概率模型(Linearprobabilitymodel,LPM)、對數(shù)單位模型(logit)、概率單位模型(probit):三種模型估計的系數(shù)大約有以下的關系:2.偏效應(1)如果解釋變量是一個連續(xù)型變量,那么他對p(x)=p(y=1
2、x)的偏效應可以通過求下面的偏導數(shù)得出來:,偏效應的符號和該解釋變量對應的系數(shù)的符號一致;兩個解釋變量偏效應之比等于它們各自的估計系數(shù)之比。
3、(2)如果解釋變量是一個離散性變量,則從變化到+1時對概率的影響大小為:上面的其他解釋變量的取值往往取其平均值。3.估計方法與約束檢驗極大似然估計;三種常見的大樣本檢驗:拉格朗日乘數(shù)檢驗、wald檢驗、似然比檢驗。4.Stata程序語法(以Probit為例)probitdepvar[indepvars][weight][ifexp][inrange][,level(#)nocoefnoconstantrobustcluster(varname)score(newvar)asisoffset(varname)maximize_options]predict[type]n
4、ewvarname[ifexp][inrange][,statisticrulesasifnooffset]wherestatisticisppredictedprobabilityofapositiveoutcome;thedefaultxblinearpredictionstdpstandarderroroftheprediction二.具體的例子1.數(shù)據(jù):美國1988年的CPS數(shù)據(jù)2.模型:估計成為工會成員的可能性,模型形式如下:參加工會的概率=F(潛在經(jīng)驗potexp、經(jīng)驗的平方項potexp2、受教育年限grade、婚否married、工會化程度high);
5、解釋變量:Potexp=年齡-受教育年限-5;grade=完成的受教育年限;married:1表示婚,0未婚;high:1表示高度工會化的行業(yè),否則為0。3.估計的結果3.1probitunionpotexppotexp2grademarriedhighunionCoef.Std.Err.zP>z[95%Conf.Interval]potexp.0835091.01560875.350.000.0529166.1141016potexp2-.0015308.0003179-4.820.000-.0021538-.0009078grade-.042078.0189089
6、-2.230.026-.0791388-.0050171married.0622516.11258360.550.580-.1584083.2829115high.5612953.0996625.630.000.3659613.7566292_cons-1.468412.2958112-4.960.000-2.048192-.88863323.2dprobitunionpotexppotexp2grademarriedhigh給出了,如果要求偏效應還需要對其乘以估計的系數(shù)beta;uniondF/dxStd.Err.zP>zx-bar[95%C.I.]potexp.0
7、226964.00415295.350.00018.884.014557.030836potexp2-.000416.000085-4.820.000519.882-.000583-.00025grade-.0114361.0051379-2.230.02613.014-.021506-.001366married*.0167881.03011370.550.580.641-.042234.07581high*.1470987.02470055.630.000.568.098687.195511obs.P.216pred.P.1904762(atx-bar)(*)dF
8、/dxisfordiscretechangeofdummyvariablefrom0to1,對離散變量。此外,如果想針對某些解釋變量的特定取值進行計算,可以用下面的語句:matrixmyx=(8,64,10,1,1)dprobit,at(myx)uniondF/dxStd.Err.zP>zx[95%C.I.]potexp.0261573.00443085.350.0008.017473.034841potexp2-.0004795.0000978-4.820.00064-.000671-.000288grade-.0131799.0065759-2.