stata建模中的各種小問題筆記

1、保存估計結(jié)果的命令：eststore名稱使用保存結(jié)果的命令：，estimates（名稱）如果你把那個顯示你用過的命令的窗口：窗口操作：WindowsReview如果你把那個顯示變量的窗口：窗口操作：WindowsVariables時間序列填充和擴展時間區(qū)間:命令：tsappend,add（n）增加n個觀測值窗口操作：在上面找dataedit即像一個表格一樣的圖標點開即可編輯數(shù)據(jù)時間序列存在間斷點問題，需要補齊處理：命令：tsfill信息準則赤池信息準則（AIC）判斷判斷模型的最大滯后階數(shù)STATA命令：1 先回歸2 estatic如何看AIC統(tǒng)計量：Breusch-Pagan,Cook-Wei

2、sberg異方差檢驗STATA命令：1 先回歸2 .estathettestvarlist或者在StatisticsPostestimation（倒數(shù)第二個）ReportsandStatistics（倒數(shù)第二個）在里面選擇（hettest）如何看統(tǒng)計量：White異方差檢驗:STATA命令：3 先回歸4 estatimtest,whitevarlist或者在StatisticsPostestimation（倒數(shù)第二個）ReportsandStatistics（倒數(shù)第二個）在里面選擇（imtest）如何看統(tǒng)計量：Ramsey回歸設(shè)定誤差檢驗：STATA命令：1先回歸2.estatovtest或者

3、在StatisticsPostestimation（倒數(shù)第二個）ReportsandStatistics（倒數(shù)第二個）在里面選擇（ovtest）如何看統(tǒng)計量：多重共線性方差膨脹因子檢驗：1 .先回歸2 .estatvif,uncentered或者在StatisticsPostestimation（倒數(shù)第二個）ReportsandStatistics（倒數(shù)第二個）一一在里面選擇（vif）如何看統(tǒng)計量：一般的當最大的方差膨脹因子超過10（相對保守的臨界值定位30）后者平均方差膨脹因子超過1表示模型存在多重共線性的問題。Uncentered用于當模型沒有常數(shù)項時的未中心化的方差膨脹因子。多重共線性的

4、其他偵查方法：R2值高而顯著的t比率?。憾嘀毓簿€性的經(jīng)典”征兆克里安經(jīng)驗法則：僅當來自一個輔助回歸的R2大于得自Y對全部回歸元中的總R2時，多重共線性才算是一個麻煩的問題。做擬合圖（前提是先回歸）STATA命令：1 .解釋變量對成分殘差圖一一用于考察模型形式是否設(shè)定準確。cprplot被解釋變量acprplot被解釋變量2 .增加變量圖一一用于考察數(shù)據(jù)是否存在異常值avplotd被解釋變量3 .擬合值對殘差圖的散點圖一一用于考察殘差是否滿足經(jīng)典的假設(shè)條件rvfplotstdp表示樣本內(nèi)預(yù)測的標 準差stdr表示樣本外預(yù)測的標 準差4 .解釋變量對殘差的散點圖rvpplot被解釋變量STATA寸

5、于數(shù)據(jù)的儲存與重現(xiàn)est命令的用法：（1）儲存回歸結(jié)果：regyx1x2x3（不限于reg,也可儲存ivreg、mvreg、reg3）eststoreA（2）重現(xiàn)回歸結(jié)果：estreplayA（3）對回歸結(jié)果進行進一步分析estforA:sum（對A回歸結(jié)果中的各個變量運行sum命令）在非時間序列的數(shù)據(jù)的情況下，異方差的修正一一用GLS具體的方法如下：1 .quietlyregressyx做回歸2 .predictu,residual取殘差3 .predictyf,xb(xb表示擬合值)將擬合值取出放到y(tǒng)f里4 .genlnu2=ln(uA2)將殘差做平方且取對數(shù)的處理5 .genyf2=yf

6、A2(將yf這個擬合值同上面的殘差做相同的處理)6 .quietlyregresslnu2yfyf2對處理過的殘差對擬合值以及處理過的擬合值做回歸7 .predictnlu2f=exp(xb()再將回歸后的擬合值取出并作對數(shù)處理放到u2f里8 .gensd=sqrt(u2f)將u2f做平方處理predictnl表示模型估計后的非線性預(yù)測，比如指數(shù)預(yù)測xb表示線性預(yù)測exp表示指數(shù)預(yù)測pr表示概率預(yù)測se表示線性預(yù)測(prediction)的標準stdf表示線性預(yù)測(forecast)的標準差取對數(shù)用ln(var)函數(shù)取平方用sqrt(var)函數(shù)然后，利用vwls進行加權(quán)估計vwlsyx,sd

7、(sd)GLS也可以通過regress命令中的weight選項來實現(xiàn)。存在自相關(guān)的修正一一用廣義差分el(mat,i,j)矩陣的第i行第j列自相關(guān)的修正一一用廣義差分具體的方法如下：1. 一階自相關(guān)的修正praisyx,rhotype(regress)praisyx,corerhotype(regress)2.高階自相關(guān)的修正一一以二階自回歸為例quietlyregressD.yx對被解釋變量取差分并且做回歸”。卬”.predictu,resid取出殘差quietlyregul(2).u,noconstant令u對其二階滯后期做自回歸(無截距)matrixmat=e(b)生成矩陣mat將回歸的

8、結(jié)果放到矩陣里(如果是更高階可能有多個自相關(guān)系數(shù))genm=d.y-el(mat,1,1)*l(2)d.y對y的一階差分與y的一階差分的滯后期與權(quán)重的乘積做差分，這個權(quán)重就是mat矩陣里的第一行第一列的系數(shù)，剛好使我們剛剛回歸出來的自相關(guān)系數(shù)(如果是高階可能不只做一個差分，會更為復(fù)雜)genn=d.x-el(mat,1,1)*l(2)d.x對x進行和y一樣的處理方法regmn然后令m對n做回歸關(guān)于參數(shù)約束的模型估計問題(P178張曉炯)STATA命令：cnsreg被解釋變量解釋變量條件ifinweight,constraints(constraints)options其中constraints(constraints)表示線性約束例如：約束規(guī)模報酬不變的估計模型命令如下：1. constraintdefine1lnk+lnl=1將等于1的線性約束(不可以為不等號)條件定義為12. cnsreglnyInklnl,constraints(1)約束為lnk=lnl=0的估計模型：命令如下：1. constraintdefine2InkInl等于0的條件可以不寫出來2. cnsreglnyInkInl,constraints(2)約束條件為l