2015-03-02 203 views
-1

我想將12個面板數據集的橫截面相關性測試結果導出到表中,以便將它們與使用不同軟件進行的類似測試進行比較。以下是xtcsd幫助頁面的迴歸和測試說明示例(不幸的是,示例數據集不可用,但xttest2頁面上的類似示例數據集tbl15-1.dta可用)。下面的說明將幫助您瞭解我想要實現:如何導出Stata xtcsd測試結果?

use "http://fmwww.bc.edu/ec-p/data/Greene2000/TBL15-1.dta" 
xtset firm year 
xtreg i f c,fe 
xtcsd, pesaran 

顯示測試統計,我可以用

return list 

如何接取該統計p值?

我發現如何使用命令esttab導出估算結果。 如何將測試結果導出到Stata中的文件?

按照下面@Maarten比伊斯的評論的p值,這裏是我如何使用低級別文件訪問導出的測試結果爲CSV文件:

file open xtcsdfile using xtcsd.csv, write replace 
file write xtcsdfile "pesaran,pvalue" _n 
file write xtcsdfile (r(pesaran)) "," (2*(normal(-abs(r(pesaran))))) _n 
file close xtcsdfile 
+0

我投票結束這個問題作爲離題,因爲它是關於將測試結果導出到文件。 – 2015-03-02 10:57:07

+1

'esttab'是一個命令,而不是一個函數。 – 2015-03-02 13:21:51

回答

2

的Pesaran統計將(漸近)跟隨標準正態分佈如果零假設爲真,那麼:p值是2*(normal(-abs(r(pesaran))))

+0

謝謝,我通過查看'viewsource xtcsd.ado'(它使用來自STATA版本6的'norm()',它已被重命名爲'normal()')來解決這個問題。我的STATA提取Pesaran測試p值的命令變成了'display 2 *(1-normal(abs(r(pesaran))))' – 2015-03-02 11:14:01

+0

我將使用'normal(-something)'而不是'1-normal ),因爲前者在數值上比後者更穩定。欲瞭解更多信息,請訪問:http://www.stata-journal.com/sjpdf.html?articlenum=pr0025 – 2015-03-02 11:55:30