【正文】 gretl 可以處理的數(shù)據(jù)類型分為 3 種：(1) 橫斷面數(shù)據(jù) (crosssectional) (2) 時間序列 (time series) (3) 追蹤數(shù)據(jù) (Panel，亦有人譯為「縱橫資料」。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 楊奕農 CYCU gretl 入門手冊 47 圖 用 gretl 匯入 Excel 數(shù)據(jù)的警告訊息。 圖 點選下載資料。在此頁的左上方，[Indices] 項下，有 [US] 和 [World] 兩個選擇，因為我們示范下載的是臺股大般指數(shù)，所以點選 [World]，進入下一網頁后，選擇其中的 [Asia/Pacific]，會出現(xiàn)如圖 的畫面。開啟之后畫面如下。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 楊奕農 CYCU gretl 入門手冊 39 從 Excel 匯入數(shù)據(jù) 假設我們要匯入的 Excel 數(shù)據(jù)文件為 的工作表 sheet1，其為 19591994 年，樣本數(shù)共 36 筆的數(shù)據(jù)，變量有 3 個，名稱分別是 obs、CT、YT (如下圖所示)。例如，我們的虛無假設是 H0: b2 = 0 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 36 gretl 入門手冊 楊奕農 CYCU 這樣的虛無假設檢定已經不需要再任何特別的步驟了，因為之前跑完回歸的結果，就已經有足夠的信息讓我們來進行判斷。顯示出來的檢定結果如下，其中 JB 值 (JarqueBera 統(tǒng)計量的簡稱) = ，其對應之 pvalue = ，所以結果和 DoornikHansen 常態(tài)性檢定相同，推論都是『無法拒絕「殘差為常態(tài)分配」的虛無假設』。 殘差無自我相關之檢定 要進行殘差無自我相關之檢定，要回到回歸結果窗口，在選單上點選 [檢定/殘差自我相關檢定]，然后在隨后出現(xiàn)的對話窗口中的「欲檢定之落后期數(shù)」 項中填入你想檢定的落后期數(shù) (此例仍為 8)， 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 32 gretl 入門手冊 楊奕農 CYCU 按確定之后會出現(xiàn)以下的檢定結果窗口。初步來看，我們會認為落后期填入 8 時，LM = ，其 p 值，P(ChiSquare(8) ) = ，據(jù)此得到無法拒絕落后 18 期都沒有自我相關異質變異之結果。而因變數(shù)與殘差無相關之檢定，在文獻中比較少看到，所以以下分別介紹如何檢定殘差是否符合假設 (2)(3)，以及 (5) 的檢定。以下先介紹常見的的回歸殘差的檢定，后節(jié)再介紹回歸系數(shù)的檢定。你若將鼠標移到十字符號上，它還會顯示該點的資料標簽，當我們想要看看哪些點是可能的「極端值」(outlier) 時，這個功能就會很有幫助。如果你希下一次執(zhí)行「一般最小平方法」時，Ct 自動被選為應變量，則可以勾選 [設定為預設變量] 這個選項。 敘述統(tǒng)計 若要看變量的敘述統(tǒng)計，只要同時按 [shift 鍵] 可復選 Ct 和 Yt 這兩個變數(shù)數(shù)據(jù)，然后鼠標指針停在陰影區(qū)之上，再按鼠標右鍵，選 [敘述統(tǒng)計] 如下圖，就可以查看所選的這兩個變量的敘述統(tǒng)計。 (5) 快速工具圖示列 在主窗口最下方，是快速工具圖示列，可以更方便地開始我們所需要的各種功能 (如下圖所示)。 . 變數(shù) 大部份對變量進行檢定、或繪圖的功能都在此項下，包含常性檢定、單根檢定、計算 Gini 系數(shù)等 . 模型 此選單下之功能提供了各種常見的計量模型，包含一般最小平方法 (Ordinary Least Squares ) 、加權最小平方法 (Weighted Least Squares)、異質變異修正模型、兩步驟最小平方法 (TwoStage Least Squares，2SLS)、高精準度之一般最小平方法、Logit、Probit、Heckit、Logistic、Possion 模型；在時間序列模型方面有自我回歸移動平均 (ARIMA)、一般化自我相關條件異質變異 (ARCH/GARCH)、向量自我回歸 (Vector Autoregression)、向量誤差修正模型 (VECM)、EngleGranger EngleGranger 共整合檢定、Johansen Johansen 共整合檢定；另外 gretl 也可以估計 Panel、非線性模型 (Nonlinear models)、最小誤差絕對誤差回歸 (LAD)、分量回歸 (quantile regression)、最大概似法、GMM 一般動差法、聯(lián)立方程式(Simultaneous equations) 等。 安裝好 gretl 之后，用鼠標在桌面的圖標上點兩下，或者是由 Windows XP 的 [開啟/所有程序/gretl] 群組中的 gretl 選項來開啟皆可。如果你想要下載的是 Windows 版本，在左邊的選單中點 gretl for Windows，然后找文中第一行的 3 2009年4月時的正式釋出的版本為 ，這個版本是我翻譯的第一個中文版本，只有選單是中文的。譯的, 在中文環(huán)境下安裝完, 自動就是中文版), 馬上就會上手操作 (加上我即將要在網站上釋出的 gretl 入門手冊來說明一些需要特別說明的功能, 在教學、學習上就更不會有問題了)。 其實我原本想要直接翻譯英文手冊，這樣的中文手冊就會自動隨附在安裝檔中。 自由軟件在使用授權上賦予軟件使用者以下的四種自由： 〔自由 0〕使用的自由： 可以不受任何限制地來使用軟件。對于經費有限的教師和學生來講，都不需要再擔心是否侵犯別人的版權。 gretl 雖然源自于 Linux 操作系統(tǒng), 但它的 Windows 版本 也十分好用, 特別是初學者幾乎不需學習任何語言指令, 只要用點選選單、填入數(shù)據(jù)的方式，就可以估計大部份的統(tǒng)計檢定和計量模型，功能十分完整，繪圖部份則由著名的 gnuplot 來支持顯示，其多樣性不在話下，相較于其它有名的商業(yè)收費統(tǒng)計軟件而言，一點也不遜色。星期的奮斗，勉強搶時間寫出這本書的 beta 版 (就是開始讓大家試用PDF 版)，以分享推廣的精神，先在網絡上公開下載，提供免費的教育目非商業(yè)性使用)。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 關于本書 這本書并非 gretl 官方英文版手冊的翻譯，它的定位只是「簡易操作手冊」，所以才命名為「gretl 入門手冊」，因此本書并非為 GNU General Public License 授權的 PDF 電子文件，但仍然歡迎教學及個人免費下載使用 (見注 1 之說明)。 所以我先行翻譯了主要的軟件操作接口之后 [注2]，再來就是決定動手寫出此書，當做公益推廣。gretl 的原始版本是在 linux 操作系統(tǒng)下所開發(fā)的，但是它的 Windows 版本也同步地在更新成長，目前 3(2010年11月) 的正式釋出的版本為 。但請注意，CVS 版是下一正式版本的測試版，內含新功能和舊 bug 的修正，但由于新功能還在測試中，所以也可能比較容易有一些新的 bug。. 資料 此項下可以增加、編輯數(shù)據(jù)變量、增加或移除樣本、進行數(shù)據(jù)類型轉換 (將轉換成時間序列 (time series)、橫剖面 (crosssectional)、或追踨資料 (panel d. 檢視 在此項下可以執(zhí)行變量繪圖、進行主成份分析 (Principal ponents)、顯示的敘述統(tǒng)計、相關系數(shù)矩陣、交叉自我相關圖等。注意，gretl 的變量只能為英文字符組成，而且英文的大小寫被視為不相同的變數(shù)。 請點選 [Ramanathan] 的數(shù)據(jù)庫，并雙點擊選擇其中的 [data36] Disposable incoand consumption 這個選項 (如下圖) 就可以開啟該數(shù)據(jù)文件了。我們可以在選單上點 [模型/一般最小平方法]，如下圖。下一章再介紹如何在 gretl 進行這些檢定。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 第 3 章 回歸模型的檢定 一般初學者在跑完回歸模型之后，通常會急著去解釋回歸估計的結果，或者馬上對回歸系數(shù)進行檢定，但其實這樣的步驟并不是正確的。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 楊奕農 CYCU gretl 入門手冊 25 若想要讓所得到之估計式也將具有常態(tài)分配的性質，則應該要加上第 5 個假設： (5) 殘差為常態(tài)性 (normality) 不過，一般來說，只要樣本數(shù)夠大，一般最小平方法所得到估計式之分配，將可漸進為常態(tài)分配 (Hill, Griffiths, and Judge, 2001)。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 楊奕農 CYCU gretl 入門手冊 27 殘差同質變異之檢定: ARCH LM 檢定 如果你的樣本數(shù)據(jù)是時間序列數(shù)據(jù)的話，通常會再進行 ARCH 類型的異質變異檢定。我們先回到回歸結果窗口，點選 [儲存/殘差平方]，然后 gretl 會要求你輸入要將殘差平方儲存成什么變量名稱，通常我們就用 gretl 給的預設名稱即可。 中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 34 gretl 入門手冊 楊奕農 CYCU 除了 DoornikHansen 常態(tài)性檢定之外，gretl 也提供了其它常態(tài)性檢定，例如很常見的 JarqueBera 常態(tài)性檢定。但由于如何修正回歸模型的計量方法，已超過本手冊的介紹范圍，無法在此深入說明，大家可自行參考計量的書籍。以下僅介紹匯入 Excel 和 Eviews 的資料文件做為示范。我們開啟 Eviews 所附的示范數(shù)據(jù) ，該檔一般是安裝在 C:/Program Files/Eviews6/Example Files/Sample Programs/logl 的子目錄下。在連上 首頁之后，點選 [INVESTING]→[Market Stats] (如圖 所示)。 然后計算機就會問你是否要儲存檔案，選 [儲存檔案] 后按 [確定]，就可以讓你中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 44 gretl 入門手冊 楊奕農 CYCU 決定數(shù)據(jù)要存到那一個子目錄。 圖 用 Excel 開啟剛剛下載的數(shù)據(jù) csv 文件。若你還中原大學國際貿易學系 Chung Yuan Christian University, Taiwan. 50 gretl 入門手冊 楊奕農 CYCU 需要人工加注一些說明，可以在主選單上，點選 [數(shù)據(jù)/編輯數(shù)據(jù)信息]，在其中輸入數(shù)據(jù)，存盤以利日后所用 (例如下圖所示)。 如果你要輸入的是橫斷面數(shù)據(jù)，就在前一個畫面的數(shù)據(jù)類型對話窗口中選擇「橫斷面數(shù)據(jù)」即可，按 [下一頁] 鈕后，看你是否要馬上接著手動輸入數(shù)據(jù)；若是，則勾選「開始輸入數(shù)據(jù)」，再按「套用」鈕，之后 gretl 會問你變量的名稱，請以英文替你的變量取一個名稱 (此例填入 y)，再按確定就會看到出現(xiàn)一個「編輯數(shù)據(jù)」窗口，我們就可以開始手動輸入資料 (如下圖所示)。 結果可以顯示成以下的樣本，最后再儲存起來即可。 我的經驗是，先用 Excel 開啟你剛剛下載的 csv 檔 (以此例而言，)，如圖 。 點此項 [^TWII] 就會進入臺股大盤指數(shù)信息主頁，要下載歷史數(shù)據(jù)的話，接著點選左方的 [Historical Prices] (見圖 )，即可開始填入你要下載的數(shù)據(jù)開始日期 (Start Date)，和結束日期 (End Date) (如圖 ；還可以選擇你要下載的是日資料 (Daily)，或者是周 (Weekly)、月 (Monthly) 資料。即使 Excel 中的數(shù)據(jù)時間順序是「遞減」而不是「遞增」的，也不用去 (重新排序) 處理它，gretl 在讀 Excel 檔時，會自動判斷來修正。另外，要再次提醒，gretl 對變量的命名，其英文大小寫是有區(qū)別的。但有的時候，可能是由于某種經濟或財務理論的關系，我們必需建立其它的虛無假設，例如想檢定看看 b2 是否等于 1，則虛無假設就要寫成 H0: b2 = 1 在 gretl 要做這種檢定，在模型結果窗口上，點選 [檢定/系