第一篇:個人總結-面板數據 stata
一、面板數據如何從混合最小二乘、固定效應、隨即效應中選擇
1、混合最小二乘、固定效應
.xtreg vol1 FI share1 share2 share3 , fe Fixed-effects(within)regression Number of obs = 289Group variable: code Number of groups = 36R-sq: within = 0.0210 Obs per group: min = 7 between = 0.0259 avg = 8.0 overall = 0.0095 max = 11 F(4,249)= 1.34corr(u_i, Xb)=-0.7580 Prob > F = 0.2569 vol1 Coef.Std.Err.t P>|t| [95% Conf.Interval] FI.0274564.0192466 1.43 0.155-.0104505.0653634 share1.0043386.1112655 0.04 0.969-.2148029.2234801 share2.0075372.111155 0.07 0.946-.2113867.2264611 share3.0058343.1113069 0.05 0.958-.2133888.2250574 _cons-.5947285 11.1241-0.05 0.957-22.50405 21.31459 sigma_u.02625677 sigma_e.05223138 rho.20172957(fraction of variance due to u_i)F test that all u_i=0: F(35, 249)= 0.87 Prob > F = 0.6879
F檢驗,原假設為個固定效應都相同,拒絕則選擇固定效應,反之為混合最小二乘?!藭r選混合最小二乘
2、隨即效應還是混合最小二乘
.xtreg roa_a stateshr size tl , reRandom-effects GLS regression Number of obs = 77Group variable: id Number of groups = 14R-sq: within = 0.4859 Obs per group: min = 3 between = 0.8053 avg = 5.5 overall = 0.6927 max = 11Random effects u_i ~ Gaussian Wald chi2(3)= 106.77corr(u_i, X)= 0(assumed)Prob > chi2 = 0.0000 roa_a Coef.Std.Err.z P>|z| [95% Conf.Interval] stateshr-.0023371.0011081-2.11 0.035-.0045089-.0001653 size.0008558.000203 4.22 0.000.0004578.0012537 tl-.1132165.0112583-10.06 0.000-.1352823-.0911507 _cons.0922799.0109761 8.41 0.000.0707672.1137926 sigma_u.00103144 sigma_e.00170742 rho.2673593(fraction of variance due to u_i)..xttest0Breusch and Pagan Lagrangian multiplier test for random effects roa_a[id,t] = Xb + u[id] + e[id,t] Estimated results: Var sd = sqrt(Var)roa_a.0000114.0033722 e 2.92e-06.0017074 u 1.06e-06.0010314 Test: Var(u)= 0 chi2(1)= 2.87 Prob > chi2 = 0.0903.原假設為var(u)=0,拒絕則采用隨即效應——此時選隨機
3、固定還是隨即.hausman fe re Coefficients(b)(B)(b-B)sqrt(diag(V_b-V_B))fe re Difference S.E.stateshr-.0208393-.0416243.0207849.022253 size.0287271.0156759.0130511.0068713 tl.0134735.2663396-.252866.2696783 b = consistent under Ho and Ha;obtained from xtreg B = inconsistent under Ha, efficient under Ho;obtained from xtreg Test: Ho: difference in coefficients not systematic chi2(3)=(b-B)'[(V_b-V_B)^(-1)](b-B)= 3.61 Prob>chi2 = 0.3070
原假設為隨機,拒絕則選固定——此時選隨機
二、縮尾會使顯著性降低
第二篇:面板數據和混合數據分析相關總結
面板數據和混合數據分析相關總結
這是我在查閱各種資料后得出的關于面板數據的總結,最近在做面板的實證論文,所以需要這個,歡迎大家繼續擴充,只要是關于面板的都行,關于具體如何在Eviews6中實現的更好,不甚感激。
*橫截面的異方差與序列的自相關性是運用面板數據模型時可能遇到的最為常見的問題,此時運用OLS可能會產生結果失真,因此為了消除影響,對我國東、中、西部地區的分析將采用不相關回歸方法(SeeminglyUnrelated Regression, SUR)來估計方程。而對于全國范圍內的估計來說,由于橫截面個數大于時序個數,所以采用截面加權估計法(Cross SectionWeights, CSW)。
*一般而言,面板數據可用固定效應(fixed effect)和隨機效應(random effect)估計方法,即如果選擇固定效應模型,則利用虛擬變量最小二乘法(LSDV)進行估計;如果選擇隨機效應模型,則利用可行的廣義最小二乘法(FGLS)進行估計(Greene ,2000)。它可以極大限度地利用面板數據的優點,盡量減少估計誤差。至于究竟是采用固定效應還是隨機效應,則要看Hausman 檢驗的結果。
*單位根檢驗:在進行時間序列的分析時,研究者為了避免偽回歸問題,會通過單位根檢驗對數據平穩性進行判斷。但對于面板數據則較少關注。隨著面板數據在經濟領域應用,對面板數據單位根的檢驗也逐漸引起重視。面板數據單位根的檢驗主要有Levin、Lin 和Chu 方法(LLC 檢驗)(1992 ,1993 ,2002)、Im、Pesaran 和Shin 方法(IPS 檢驗)(1995 ,1997)、Maddala 和Wu 方法(MW檢驗)(1999)等。
*協整檢驗:協整檢驗是考察變量間長期均衡關系的方法。在進行了各變量的單位根檢驗后,如果各變量間都是同階單整,那么就可以進行協整檢驗了。面板協整檢驗理論目前還不成熟,仍然在不斷的發展過程中,目前的方法主要有:
(1)Kao(1999)、Kao and Chiang(2000)利用推廣的DF和ADF檢驗提出了檢驗面板協整的方法,這種方法零假設是沒有協整關系,并且利用靜態面板回歸的殘差來構建統計量。
(2)Pedron(i1999)在零假設是在動態多元面板回歸中沒有協整關系的條件下給出了七種基于殘差的面板協整檢驗方法。和Kao的方法不同的是,Pedroni的檢驗方法允許異質面板的存在。
(3)Larsson et a(l2001)發展了基于Johansen(1995)向量自回歸的似然檢驗的面板協整檢驗方法。這種檢驗的方法是檢驗變量存在共同的協整的秩。
*一般的順序是:先檢驗變量的平穩性,當變量均為同階單整變量時,再采用協整檢驗以判別變量間是否存在長期均衡關系。如果變量間存在長期均衡的關系,我們可以通過誤差修正模型(ECM)來檢驗變量間的長期因果關系;如變量間不存在協整關系,我們將對變量進行差分,然后通過向量自回歸模型(VAR),檢驗變量間的短期因果關系。
關于平穩性檢驗和協整檢驗、因果檢驗流程圖
↗ 同階單整→協整檢驗→協整?(YES:EG兩步法 for 長期因果關系;NO:誤差修正模型ECM/VEC for 短期因果關系)
平穩?(單位根檢驗)
↘非同階單整→差分使平穩→VAR→Granger因果檢驗 for 短期因果關系
關于面板數據模型選擇回歸與檢驗流程圖
混合固定(main:個體固定)隨機(main:個體隨機)▏▁▁▁▁▁▁▁▁▁▁▁▁▏▁▁▁▁▁▁▁▁▁▁▁▁▁ ▏
▏先回歸估計▏先回歸估計
↓Cross-section:fixed↓Cross-section:random
F檢驗Hausman檢驗
▏▏
H0:混合H1:個體固定HO:個體隨機H1:個體固定
--
Output:▏▏
If:If:
F=(Cross-section F Stat.)>Fa(df1,df2)H=(Cross-section Random Stat.)>χ2a(df1)
or Prob. Then:reject H0,accept H1Then:reject H0,accept H 1是先做F檢驗還是先做Hausman檢驗?。蛔鯢檢驗的時候,Fixed and Random、comm和Cross-section specific選項應該怎么設置??;另外我看高鐵梅上面對面板的分類有些不同,能說說有啥區別么? 以Eviews6為例,來說明一下面板模型的選擇問題: F檢驗是用來在混合模型和固定效應模型中做出選擇,而Hausman檢驗是用來在固定效應模型和隨機效應模型中做出選擇,所以不存在孰先孰后的問題; 由于我們通常估計的個體效應而不是時刻效應,所以我們進行回歸和檢驗的時候,Period選擇None。 回歸的時候,具體操作設置如下,Depedent Variable里填因變量,Common Coefficients里填自變量(包括截距項c),Cross-Section視回歸需要選擇None、Fixed、Random,Period選擇None,可以依次實現混合回歸、個體固定回歸、個體隨機回歸。然后在個體固定回歸之后,進行F檢驗進行模型選擇決策1;在個體隨機回歸之后,進行Hausman檢驗進行模型選擇決策2,從而最終得出最佳回歸。 表2-1: 回歸分析相關命令一覽 命令用途 anova 方差和協方差分析 heckman Heckman 篩選模型 intreg 離散型變量模型,包括Tobit、cnreg 和intreg ivreg 工具變量法(IV 或2SLS) newey Newey-West 標準差設定下的回歸 prais 針對序列相關的Prais-Winsten, Cochrane-Orcutt, or Hildreth-Lu 回歸 qreg 分量回歸 reg OLS 回歸 sw 逐步回歸法 reg3 三階段最小二乘回歸 rreg 穩健回歸(不同于方差穩健型回歸,即White 方法)sureg 似無相關估計 svyheckman 調查數據的Heckman 篩選模型 svyintreg 調查數據的間斷變量回歸 svyregress 調查數據的線性回歸 tobit Tobit 回歸 treatreg treatment 效應模型 truncreg 截斷回歸 表2-2: 時間序列命令一覽 命令用途 clemao1 允許結構突變的單位根檢驗 zandrews dfuller dfgls pperron coin 單方程協整檢驗 dwstat 參考dwstat2 , durbina2 durbinh 表2-3: Panel Data 模型相關命令一覽I 命令模型 統計描述相關命令: xtdes 變量類型,數據類型描述 xtsum 基本統計量 xttab 按表格形式列示 xtpattern 面板數據的模式 估計相關命令: xtreg 面板數據模型(固定效應、隨機效應) xtregar 含有AR(1)干擾項的固定效應和隨機效應面板數據模型 xtgls 截面-時序混合模型,可處理異方差、組內序列相關和組間相關性 xtpcse OLS or Prais-Winsten models with panel-corrected standard errors xtrchh Hildreth-Houck random coefficients models xtivreg 面板模型的工具變量或兩階段最小二乘法估計 xtabond Arellano-Bond(1991)線性動態面板數據模型估計 xtabond2 Arellano-Bover(1995)系統GMM 動態面板數據模型估計 xttobit Tobit 隨機效應面板模型 xtintreg Random-effects interval data regression models xtlogit Fe, Re, Pa logit models xtprobit Re, Pa probit models xtcloglog Re, Pa cloglog models xtpoisson Fe, Re, Pa Poisson models xtnbreg Fe, Re, Pa negative binomial models xtfrontier 面板隨機前沿模型 xthtylor Hausman-Taylor estimator for error-components models 表2-4: Panel Data 模型相關命令一覽II 命令模型 假設檢驗相關: test Wald 檢驗,如時間效應聯合顯著性檢驗 xttest0 隨機效應檢驗 xttest1 面板序列相關檢驗 xttest2 ads xtserial Wooldridge 一階序列相關檢驗 xtab Arellano 面板一階序列相關檢驗 hausman Hausman 檢驗 面板單位根和協整相關: xtunit stata提供的檢驗方法 ipshin IPS(2003)面板單位根檢驗 levilin Levin,Lin和Chu(LLC, 2002)面板單位根檢驗 madfuller Sarno-Taylor(1998)面板單位根檢驗 xtfisher Maddala和Wu(1999),基于P 值的面板單位根檢驗 表2-5: Post-estimation Commands 命令名稱用途 adjust 列示預測結果的均質,適于多種回歸分析,可分組列示 estimates 估計結果的存儲、再顯示、列表比較等 hausman Hausman 模型識別檢驗 lincom 獲得參數的線性組合,在Logit 模型中可以獲得系數線性組合的OR 值 linktest 但方程link 識別檢驗,用y 對Oy 和Oy2 回歸 lrtest 似然比(LR)檢驗 mfx 計算邊際效應和彈性系數 nlcom 系數的非線性組合 predict 獲得擬合值、殘差等 predictnl 獲得非線性估計的擬合值、殘差等 test 線性約束的假設檢驗,Wald 檢驗 testnl 非線性約束的假設檢驗 vce 列示參數估計值的方差-協方差矩陣 表2-6: 二維圖種類一覽 圖形種類簡單描述 scatter scatterplot line line plot connected connected-line plot scatteri scatter with immediate arguments area line plot with shading bar bar plot spike spike plot dropline dropline plot dot dot plot rarea range plot with area shading rbar range plot with bars rspike range plot with spikes rcap range plot with capped spikes rcapsym range plot with spikes capped with symbols rscatter range plot with markers rline range plot with lines rconnected range plot with lines and markers tsline time-series plot tsrline time-series range plot mband median-band line plot mspline spline line plot lowess LOWESS line plot lfit linear prediction plot qfit quadratic prediction plot fpfit fractional polynomial plot lfitci linear prediction plot with CIs qfitci quadratic prediction plot with CIs fpfitci fractional polynomial plot with CIs function line plot of function histogram histogram plot kdensity kernel density plot 表2-7: 二維圖選項一覽 選項類別簡單描述 added line options draw lines at specified y or x values added text option display text at specified(y,x)value axis options labels, ticks, grids, log scales title options titles, subtitles, notes, captions legend option legend explaining what means what scale(#)resize text, markers, and line widths region options outlining, shading, aspect ratio, size aspect option constrain aspect ratio of plot region scheme(schemename)overall look by(varlist,...)repeat for subgroups nodraw suppress display of graph name(name,...)specify name for graph saving(filename,...)save graph in file advanced options difficult to explain 表2-9: 模擬分析相關命令一覽 命令用途備注 抽樣相關: corr2data 產生具有指定相關性的數據僅適用于模擬相關分析 drawnorm invnorm(uniform())產生服從標準正態分布的隨機數函數,可調節均值和方差 matuniform(r,c)產生均勻分布函數 sample 從現有數據中進行非重復隨機抽樣參考bsample sim arma 產生服從ARIMA 過程的隨機變量需要下載 Bootstrap 相關: bootstrap bs bstat bsample MC 相關: simulate MC simulation jknife 類似于MC permute postfile 存儲MC 的結果 statsby exp list stata11常用命令 注:JB統計量對應的p大于0.05,則表明非正態,這點跟sktest和swilk檢驗剛好相反; dta為數據文件; gph為圖文件; do為程序文件; 注意stata要區別大小寫; 不得用作用戶變量名: _all _n _N _skip _b _coef _cons _pi _pred _rc _weight double float long int in if using with 命令: 讀入數據一種方式 input x y 1 4 2 5.5 3 6.2 4 7.7 5 8.5 end su/summarise/sum x 或 su/summarise/sum x,d 對分組的描述: sort group by group:su x %%%%% tabstat economy,stats(max)%返回變量economy的最大值 %%stats括號里可以是:mean,count(非缺失觀測值個數),sum(總和),max,min,range,%% sd,var,cv(變易系數=標準差/均值),skewness,kurtosis,median,p1(1%分位 %% 數,類似地有p10, p25, p50, p75, p95, p99),iqr(interquantile range = p75 – p25)_all %描述全部 _N 數據庫中觀察值的總個數。_n 當前觀察值的位置。_pi 圓周率π的數值。list gen/generate %產生數列 egen wagemax=max(wage)clear use by(分組變量)set more 1/0 count %計數 gsort +x(升序)gsort-x(降序)sort x 升序;并且其它變量順序會跟著改變 label var y “消費” %添加標簽 describe %描述數據文件的整體,包括觀測總數,變量總數,生成日期,每個變量的存儲類型(storage type),標簽(label)replace x5=2*y if x!=3 %替換變量值 replace age = 25 in 107 %令第107個觀測中age為25 rename y2 u %改變變量名 drop in 2 %刪除全部變量的第2行 drop if x==.刪去x為缺失值的所有記錄 keep if x<2 %保留小于2的數據,其余變量跟隨x改變 keep in 2/10 %保留第2-10個數 keep x1-x5 %保留數據庫中介于x1和x5間的所有變量(包括x1和x5),其余變量刪除 ci x1 x2,by(group)%算出置信區間,不過先前對group要先排序,即sort group; %by的意思逐個進行 cii 12 3.816667 0.2710343, level(90)%已知均值,方差,計算90%的置信區間 cii 10 2 %obs=10,mean=2,以二項分布形式,計算置信區間 centile x,centile(2.5 25 50 75 97.5)%取分位數 correlate/corr x y z %相關系數 pwcorr x y,sig %給出原假設r=0的命令 %如果變量非服從正態分布,則spearman x y regress/reg mean year %回歸方程建立 reg y x,noconstant %無常數項 predict meanhat %預測擬合值 predict e,residual %得到殘差 estat hettest % 異方差檢驗 dwstat % Durbin-Watson自相關檢驗 vif % 方差膨脹因子 logit y x1 x2 x3(y取0或1,是被解釋變量,x1-x3是被解釋變量)%logit回歸 probit y x1 x2 x3(y取0或1,是被解釋變量,x1-x3是被解釋變量)%probit回歸 tobit y x1 x2 x3(y取值在0和1之間,是被解釋變量,x1-x3是被解釋變量)%tobit回歸 sktest e %殘差正態性檢驗 p>0.05則接受原假設,即服從正態分布; %% sktest是基于變量的偏度和斜度(正態分布的偏度為0,斜度為3)swilk x %基于Shapiro-Wilk檢驗 %%p值越小,越傾向于拒絕零假設,也就是變量越有可能不服從正態分布 xi %生成虛擬變量 tabulat gender,summ(math)%用gender指標對math進行分類,返回兩類math的mean、std、freq tabulate=tab %gen f=int((shengao-164)/3)*3+164 組距為3 tabulate 變量名 [, generate(新變量)missing nofreq nolabel plot ] %%%%% generate(新變量)// 按分組變量產生啞變量 nofreq // 不顯示頻數 nolabel // 不顯示數值標記 plot // 顯示各組頻數圖示 missing // 包含缺失值 cell // 顯示各小組的構成比(小組之和為 1)column // 按欄顯示各組之構成(各欄總計為 1)row // 按行顯示各組之構成(各行總計為 1)%%%%% 求和,求最??? mod(x,y)%求余數 means %返回三種平均值 di normprob(1.96)di invnorm(0.05)di binomial(20,5,0.5)di invbinomial(20,5,0.5)di tprob(10,2)di invt(10.0.05)di fprob(3,27,1)di invfprob(3,27,0.05)di chi2(3,5)di invchi2(3,0.05)stack x y z,into(e)%把三列合成一列 xpose,clear %矩陣轉置 append using d:
主站蜘蛛池模板:
亚洲精品国产成人精品|
极品少妇被猛的白浆直喷白浆|
两口子交换真实刺激高潮|
亚洲日本va午夜在线电影|
精品日产一卡2卡三卡4卡在线|
亚洲国产五月综合网|
国产女人好紧好爽|
日日摸夜夜添夜夜添无码免费视频|
免费无码av片在线观看播放|
精品人妻系列无码人妻不卡|
国产精品9999久久久久仙踪林|
国产女人高潮叫床免费视频|
国产精品xxx大片免费观看|
成人无遮挡18禁免费视频|
国产福利视频一区二区|
日本丰满熟妇videossex8k|
国产亚洲精品一区二区三区|
伊人久久精品无码av一区|
日韩亚洲国产激情一区浪潮av|
国内精品视频一区二区三区八戒|
日韩欧美一区二区三区永久免费|
8x国产精品视频|
中文字幕乱码亚洲无线码|
国产美女精品视频线免费播放|
国产精品进线69影院|
亚洲日产精品一二三四区|
女人爽到高潮潮喷18禁网站|
呻吟国产av久久一区二区|
99re6热在线精品视频播放|
国产精品福利视频导航|
亚洲 欧美 日韩 国产综合 在线|
久久国产精品99精品国产|
永久免费观看国产裸体美女|
国产午夜高潮熟女精品av|
狂野欧美性猛xxxx乱大交|
亚洲av无码久久精品蜜桃|
亚洲人成绝网站色www|
精品卡一卡二卡3卡高清乱码|
久久99热人妻偷产国产|
2020国产精品精品国产|
国产精品区一区二区三在线播放|
第三篇:stata命令總結
第四篇:stata命令總結