Stata-簡介及基本操作PPT

1、1Stata 簡介及基本操作簡介及基本操作一、為什么使用一、為什么使用Stata？ Stata 是目前在歐美最為流行的計量軟件，具有操作是目前在歐美最為流行的計量軟件，具有操作簡單、功能強大的特點。由于使用簡單、功能強大的特點。由于使用Stata 的用戶很多，對的用戶很多，對于最新的計量方法，常?？梢韵螺d由用戶寫的于最新的計量方法，常?？梢韵螺d由用戶寫的Stata 命令命令程序程序,十分方便。而官方的十分方便。而官方的Stata版本也經(jīng)常更新，以適應(yīng)版本也經(jīng)常更新，以適應(yīng)計量經(jīng)濟學(xué)迅猛發(fā)展的需要。計量經(jīng)濟學(xué)迅猛發(fā)展的需要。 Stata 11 已于已于2009 年出版，但由于在中國普遍使用年出版

2、，但由于在中國普遍使用的仍是的仍是Stata 10，我們主要介紹，我們主要介紹Stata 10。21. 面臨的問題面臨的問題 研究者經(jīng)常地希望利用多種方法從不同角度對數(shù)據(jù)進研究者經(jīng)常地希望利用多種方法從不同角度對數(shù)據(jù)進行分析，但擁有絕對優(yōu)勢的軟件是不存在的，不存在哪個行分析，但擁有絕對優(yōu)勢的軟件是不存在的，不存在哪個軟件能解決所有的問題。軟件能解決所有的問題。 計量經(jīng)濟學(xué)發(fā)展迅速，新的方法層出不窮，而軟件公計量經(jīng)濟學(xué)發(fā)展迅速，新的方法層出不窮，而軟件公司將這些新方法嵌入到軟件中是存在滯后的。司將這些新方法嵌入到軟件中是存在滯后的。2. 解決方案解決方案 Stata以其操作簡單、運行快速而精確，

3、全面的數(shù)據(jù)管以其操作簡單、運行快速而精確，全面的數(shù)據(jù)管理功能，精致的作圖，強大的統(tǒng)計與計量分析功能，簡練理功能，精致的作圖，強大的統(tǒng)計與計量分析功能，簡練標準的程序與矩陣運算功能、豐富的網(wǎng)絡(luò)資源，在統(tǒng)計學(xué)、標準的程序與矩陣運算功能、豐富的網(wǎng)絡(luò)資源，在統(tǒng)計學(xué)、經(jīng)濟學(xué)、金融學(xué)、心理學(xué)、計算機科學(xué)、物理、化學(xué)等多經(jīng)濟學(xué)、金融學(xué)、心理學(xué)、計算機科學(xué)、物理、化學(xué)等多個學(xué)科領(lǐng)域得到廣泛使用。個學(xué)科領(lǐng)域得到廣泛使用。請瀏覽：請瀏覽：http:/ STATA簡介：發(fā)展歷史簡介：發(fā)展歷史 （1）Stata是一個用于數(shù)據(jù)管理和分析的統(tǒng)計軟件，由是一個用于數(shù)據(jù)管理和分析的統(tǒng)計軟件，由Stata公司（公司（Stata

4、Corp）研制開發(fā)。從）研制開發(fā)。從1985開始，先后開開始，先后開發(fā)了發(fā)了1.1, 1.2,等多個版本。目前等多個版本。目前Stata的最新版本為的最新版本為2009年年7月份推出的月份推出的Stata 11.0。 （2）Stata 11包括四種版本：包括四種版本：Small（小型版）、（小型版）、IC（標（標準版）、準版）、SE（特別版）和（特別版）和MP（多處理器版）。其中屬（多處理器版）。其中屬MP版本最為強大。版本最為強大。MP版與版與SE版的功能完全相同，但版的功能完全相同，但MP版的版的運算速度比運算速度比SE版的要快很多。用戶可以在命令欄輸入版的要快很多。用戶可以在命令欄輸入a

5、bout 查看所安裝的查看所安裝的Stata所屬的版本。所屬的版本。4 （3）不同版本對樣本容量、變量個數(shù)、矩陣階數(shù)、宏）不同版本對樣本容量、變量個數(shù)、矩陣階數(shù)、宏的字符長度等有著不同的限制。以的字符長度等有著不同的限制。以SE版為例，其最大變版為例，其最大變量個數(shù)為量個數(shù)為32767，最大字符長度為，最大字符長度為244字節(jié)，最大矩陣階數(shù)字節(jié)，最大矩陣階數(shù)為為11000（即（即1100011000）。）。 Stata默認值為：最大變量個數(shù)為默認值為：最大變量個數(shù)為5000，最大矩陣階數(shù)，最大矩陣階數(shù)為為400，最大內(nèi)存為，最大內(nèi)存為50兆兆(在在Stata10中為中為10m)。如果用戶。如果

6、用戶需要更多的內(nèi)存或者更多的變量，可以在命令欄輸入如下需要更多的內(nèi)存或者更多的變量，可以在命令欄輸入如下命令進行擴展。命令進行擴展。. set maxvar # (.set maxvar 5000 ). set maxsize # (.set maxsize 1000 ). set memory #m (.set memory 50m )54. STATA簡介：主要特點與功能簡介：主要特點與功能統(tǒng)計分析與計量分析的結(jié)合：統(tǒng)計分析與計量分析的結(jié)合： （1）單元統(tǒng)計：描述統(tǒng)計、假設(shè)檢驗（參數(shù)、非參數(shù)）、）單元統(tǒng)計：描述統(tǒng)計、假設(shè)檢驗（參數(shù)、非參數(shù)）、ANOVA、質(zhì)量控制、統(tǒng)計作圖。、質(zhì)量控制、統(tǒng)計

7、作圖。 （2）多元統(tǒng)計：）多元統(tǒng)計：MANOVA、主成分、因子分析、典型相、主成分、因子分析、典型相關(guān)、聚類、判別分析、對應(yīng)分析、多維標度。關(guān)、聚類、判別分析、對應(yīng)分析、多維標度。 （3）線性回歸、非線性回歸、工具變量回歸、廣義線性回）線性回歸、非線性回歸、工具變量回歸、廣義線性回歸、分位數(shù)回歸（穩(wěn)健回歸）、系統(tǒng)方程模型（歸、分位數(shù)回歸（穩(wěn)健回歸）、系統(tǒng)方程模型（SUR、聯(lián)立、聯(lián)立方程）、離散選擇模型（二項選擇、排序選擇、多項選擇、方程）、離散選擇模型（二項選擇、排序選擇、多項選擇、條件條件Logit、嵌套、嵌套Logit模型、二元選擇模型等）、計數(shù)模型模型、二元選擇模型等）、計數(shù)模型（泊松回

8、歸、負二項回歸）、截斷與歸并模型、海克曼選擇（泊松回歸、負二項回歸）、截斷與歸并模型、?？寺x擇模型、逐步回歸模型、逐步回歸(stepwise)等。等。6 （4）時間序列分析：時間序列的平滑、相關(guān)圖、）時間序列分析：時間序列的平滑、相關(guān)圖、ARIMAX、GARCH、單位根檢驗、單位根檢驗、Johansen協(xié)整檢驗、協(xié)整檢驗、VAR、VEC、滾動回歸等。、滾動回歸等。 （5）面板數(shù)據(jù)（線性模型、工具變量回歸、動態(tài)面板、分）面板數(shù)據(jù)（線性模型、工具變量回歸、動態(tài)面板、分層混合效應(yīng)、廣義估計方程（層混合效應(yīng)、廣義估計方程（GMM）、隨機邊界模型等）。）、隨機邊界模型等）。 （6）調(diào)查數(shù)據(jù)分析）調(diào)查數(shù)

9、據(jù)分析 （7）生存分析）生存分析 （8）蒙特卡洛模擬、自舉法）蒙特卡洛模擬、自舉法(bootstrap)、刀切法、刀切法（jackknife）7豐富的網(wǎng)絡(luò)資源：豐富的網(wǎng)絡(luò)資源： http:/ http:/www.stata- http:/www.stata-嚴謹、簡練、靈活的程序語言：嚴謹、簡練、靈活的程序語言： Stata提供了嚴謹、簡練而靈活的程序語句，用戶可以編提供了嚴謹、簡練而靈活的程序語句，用戶可以編寫自己的命令和函數(shù)，也可以制作自己的對話框和窗口菜寫自己的命令和函數(shù)，也可以制作自己的對話框和窗口菜單。單。強大的矩陣運算功能：強大的矩陣運算功能：Mata Mata是是Stata的矩陣

10、運算語言，實現(xiàn)各種運算。的矩陣運算語言，實現(xiàn)各種運算。兼容性：兼容性： 在在Stata中可以直接運行其它程序；其他語言編寫的程序中可以直接運行其它程序；其他語言編寫的程序可以作為插件（可以作為插件（Plugin）直接嵌入到）直接嵌入到Stata中。中。8二、二、Stata 的窗口的窗口 安裝好安裝好 Stata 后，點擊電腦桌面上的后，點擊電腦桌面上的Stata 圖標，即可圖標，即可打開打開Stata。此時可以看到，在最上方有一排菜單，即。此時可以看到，在最上方有一排菜單，即“File Edit Data Graphics Statistics User Window Help”。在此之下，有

11、四個窗口，分別為：在此之下，有四個窗口，分別為： 左上左上“Review”（歷史窗口）：此窗口記錄著自啟（歷史窗口）：此窗口記錄著自啟動動Stata以來執(zhí)行過的命令。以來執(zhí)行過的命令。 左下左下“Variables”（變量窗口）：此窗口記錄著目前（變量窗口）：此窗口記錄著目前Stata 內(nèi)存中的所有變量。內(nèi)存中的所有變量。 右上右上“Results”（結(jié)果窗口）：此窗口顯示執(zhí)行（結(jié)果窗口）：此窗口顯示執(zhí)行Stata 命令后的輸出結(jié)果。命令后的輸出結(jié)果。 右下右下“Command”（命令窗口）：在此窗口輸入想（命令窗口）：在此窗口輸入想要執(zhí)行的要執(zhí)行的Stata 命令。命令。9 為了使屏幕分割更

12、美觀實用，可以用鼠標將以上窗口為了使屏幕分割更美觀實用，可以用鼠標將以上窗口拉到任意大小與位置。然后點擊菜單拉到任意大小與位置。然后點擊菜單Edit”“Preferences” “General Preferences” “Windowing” “Locksplitter”，就可以鎖定當前畫面，而在以后重啟，就可以鎖定當前畫面，而在以后重啟Stata 時時自動顯示這個畫面設(shè)置。自動顯示這個畫面設(shè)置。10Stata 的主界面的主界面11三、三、Stata 操作操作1將數(shù)據(jù)導(dǎo)入將數(shù)據(jù)導(dǎo)入Stata 打開打開Stata 軟件后，點擊軟件后，點擊Data Data Editor 圖標圖標,即可打即可打

13、開一個類似開一個類似Excel的空白表格。然后，用的空白表格。然后，用Excel 打開文件打開文件“實驗數(shù)據(jù)實驗數(shù)據(jù).xls”，復(fù)制文件中的相關(guān)數(shù)據(jù)，并粘貼到，復(fù)制文件中的相關(guān)數(shù)據(jù)，并粘貼到Data Editor 中。中。 導(dǎo)入數(shù)據(jù)的另一方法是，點擊菜單導(dǎo)入數(shù)據(jù)的另一方法是，點擊菜單“File”“Import”，然后導(dǎo)入各種格式的數(shù)據(jù)。但這種方法有時不如直接從然后導(dǎo)入各種格式的數(shù)據(jù)。但這種方法有時不如直接從Excel 表中粘貼數(shù)據(jù)來得方便直觀。表中粘貼數(shù)據(jù)來得方便直觀。12 關(guān)閉關(guān)閉 Data Editor 后，即會看到左下方的后，即會看到左下方的“Variables”窗口出現(xiàn)了窗口出現(xiàn)了8

14、個變量，分別為個變量，分別為year、gdp、index、index1978、gov、pop、deflator、gcons。13 此時，可以點擊此時，可以點擊 Save 圖標（也可以點擊菜單圖標（也可以點擊菜單“File”“Save”），將數(shù)據(jù)存為），將數(shù)據(jù)存為Stata 格式的文件（擴展名為格式的文件（擴展名為dta），），比如比如wanger_law.dta。 這樣，以后就可以用這樣，以后就可以用Stata 直接打開這個數(shù)據(jù)集了（不直接打開這個數(shù)據(jù)集了（不需要再從需要再從Excel 表中粘貼過來）。表中粘貼過來）。 打開的方式有兩種。可以點擊打開的方式有兩種。可以點擊Open 圖標（也可以

15、點擊圖標（也可以點擊菜單菜單“File”“Open”），然后尋找要打開的），然后尋找要打開的dta 文件的位文件的位置。置。 另外一種方法是在命令窗口輸入以下命令（假設(shè)文件在另外一種方法是在命令窗口輸入以下命令（假設(shè)文件在E 盤的根目錄）并回車（按盤的根目錄）并回車（按Enter 鍵），鍵），use E:wagner_law.dta,clear其中，選擇項其中，選擇項“clear”表示可以替代內(nèi)存中的已有數(shù)據(jù)。表示可以替代內(nèi)存中的已有數(shù)據(jù)。14 如果要關(guān)閉一個數(shù)據(jù)集（如果對數(shù)據(jù)集進行了改動，別如果要關(guān)閉一個數(shù)據(jù)集（如果對數(shù)據(jù)集進行了改動，別忘了先存盤忘了先存盤Save），以便使用另外一個數(shù)據(jù)集

16、，可以在命），以便使用另外一個數(shù)據(jù)集，可以在命令窗口輸入：令窗口輸入：. clear 這樣，內(nèi)存中所有的當前數(shù)據(jù)都被清空，然后可以再打這樣，內(nèi)存中所有的當前數(shù)據(jù)都被清空，然后可以再打開另外一個數(shù)據(jù)集。開另外一個數(shù)據(jù)集。152變量的標簽變量的標簽 在變量窗口，每個變量的在變量窗口，每個變量的“名字名字”（Name）旁邊顯示了）旁邊顯示了其其“標簽標簽”（label）。但目前的標簽過于簡略，缺乏變量的）。但目前的標簽過于簡略，缺乏變量的解解釋信息。釋信息。 如果想將變量如果想將變量“gov”的標簽改為的標簽改為“government expenditures”，可進行如下操作。以鼠標的右鍵點擊變量

17、名，可進行如下操作。以鼠標的右鍵點擊變量名“gov”，然后選擇，然后選擇“Edit variable label”，輸入，輸入“government expenditures”即可。此時，再去看變量即可。此時，再去看變量“gov”的標簽，就已的標簽，就已經(jīng)改為經(jīng)改為“government expenditures”。另外，在右上角的結(jié)果。另外，在右上角的結(jié)果窗口出現(xiàn)了以下一行命令：窗口出現(xiàn)了以下一行命令：. label variable tc government expenditures 這表明，直接在命令窗口輸入上面這行命令，也能達到一這表明，直接在命令窗口輸入上面這行命令，也能達到一樣的

18、效果，但不如直接在變量上點擊右鍵操作方便。樣的效果，但不如直接在變量上點擊右鍵操作方便。注意：注意：Stata 中字母的大小寫是嚴格區(qū)分的，因此中字母的大小寫是嚴格區(qū)分的，因此Stata 建建議對于變量名一律使用小寫字母。議對于變量名一律使用小寫字母。163審視數(shù)據(jù)審視數(shù)據(jù) 一個數(shù)據(jù)集可能很大，而我們常希望看到數(shù)據(jù)的概貌。一個數(shù)據(jù)集可能很大，而我們常希望看到數(shù)據(jù)的概貌。想看數(shù)據(jù)集中的變量名單、標簽等，可以在命令窗口輸入：想看數(shù)據(jù)集中的變量名單、標簽等，可以在命令窗口輸入：. describe 如果想給整個數(shù)據(jù)集加上一個標簽，以說明這個數(shù)據(jù)如果想給整個數(shù)據(jù)集加上一個標簽，以說明這個數(shù)據(jù)集來自集來

19、自“Wagner Law 1978-2009”，可輸入命令：，可輸入命令：. label data Wagner Law 1978-2009如果再次使用命令如果再次使用命令“describe”，就會看到數(shù)據(jù)集的標簽，就會看到數(shù)據(jù)集的標簽“Wagner Law 1978-2009”。如果想看變量如果想看變量 gov、gcons與與gdp的具體數(shù)據(jù)，可使用命令的具體數(shù)據(jù)，可使用命令:. list gov gcons gdp注意：注意：如果數(shù)據(jù)較長，可用如果數(shù)據(jù)較長，可用Ctrl+Break終止顯示。終止顯示。17 有時我們想對數(shù)據(jù)集的一部分執(zhí)行命令，比如只想有時我們想對數(shù)據(jù)集的一部分執(zhí)行命令，比如

20、只想看變量看變量gov 與與gdp 的前的前5 個數(shù)據(jù)，則可輸入命令個數(shù)據(jù)，則可輸入命令:. list gov gdp in 1/5 也可以通過邏輯關(guān)系來定義數(shù)據(jù)集的子集。如果要也可以通過邏輯關(guān)系來定義數(shù)據(jù)集的子集。如果要列出所有滿足條件列出所有滿足條件“year 1995”的變量的變量gov與與gdp 的數(shù)的數(shù)據(jù)，則可以使用以下命令據(jù)，則可以使用以下命令:. list year gov gdp if year=1995其中，其中，“=”表示表示“大于等于大于等于”。其他表示關(guān)系的邏輯。其他表示關(guān)系的邏輯符號為符號為“= =”（等于），（等于），“”（大于），（大于），“”（小于），（小于），

21、“=2001反之，如果只想保留滿足反之，如果只想保留滿足“year2001”條件的觀測值，條件的觀測值，而刪去所有其他觀測值而刪去所有其他觀測值:. keep if year=2001195考察變量的統(tǒng)計特征考察變量的統(tǒng)計特征如果想看變量如果想看變量 gov、gcons和和gdp的統(tǒng)計特征，可輸入命令的統(tǒng)計特征，可輸入命令:. summarize gov gcons gdpVariable Obs Mean Std. Dev. Min Maxgov 3214454.13 19217.2 1122.09 76299.93gcons 3211441.76 12730.41 480 44396.9g

22、dp 3280670.91 94125.313645.218 340506.9滿足條件的統(tǒng)計：滿足條件的統(tǒng)計：. summarize gov gcons gdp if year2000如果不指明變量，則將顯示數(shù)據(jù)集中所有變量的統(tǒng)計指標。如果不指明變量，則將顯示數(shù)據(jù)集中所有變量的統(tǒng)計指標。summarize如果要顯示內(nèi)存中某些變量之間的相關(guān)系數(shù)，可輸入命令如果要顯示內(nèi)存中某些變量之間的相關(guān)系數(shù)，可輸入命令:. correlate gov gcons pop gdp205畫圖畫圖 Stata 具有很強的畫圖功能。如果想看變量具有很強的畫圖功能。如果想看變量q 的直方圖的直方圖（假定組寬為（假定組寬

23、為1000），可輸入以下命令），可輸入以下命令:. histogram q, width(8000) frequency05101520Frequency020000400006000080000government expenditures210.00001.00002.00003.00004Density020000400006000080000government expenditureskernel = epanechnikov, bandwidth = 6128.97Kernel density estimate 由于直方圖不連續(xù)，如果想看連續(xù)的經(jīng)驗分布圖，可使由于直方圖不連續(xù)，如果

24、想看連續(xù)的經(jīng)驗分布圖，可使用以下命令：用以下命令：. kdensity gov22如果要畫如果要畫 gov 與與gdp 之間的散點圖，則可輸入以下命令之間的散點圖，則可輸入以下命令:scatter gov gdp020000400006000080000government expenditures0100000200000300000400000gdp23如果想在散點圖上同時畫出回歸直線，可以使用如下命令如果想在散點圖上同時畫出回歸直線，可以使用如下命令:. twoway (scatter gov gdp)(lfit gov gdp)02000040000600008000001000002

25、00000300000400000gdpgovernment expendituresFitted values其中，其中，“l(fā)fit”表示表示“l(fā)inear fit”（線性擬合）。（線性擬合）。246生成新變量生成新變量 如果要估計以下方程：如果要估計以下方程：tttgdpgov )ln()ln(10為了估計這個方程，需要在為了估計這個方程，需要在 Stata 中對原變量取對數(shù)，中對原變量取對數(shù)，可使用命令可使用命令“generate”。. generate lngov=log(gov). generate lngdp=log(gdp). reg lngov lngdp如果要生成如果要生成“

26、虛擬變量虛擬變量”,定義為：定義為：1995119950ttdumtStata命令為：命令為：generate dum=(if year=1995)25刪除變量：刪除變量：drop 變量名變量名如：如：drop dum修改變量值修改變量值replace 變量名變量名=表達式表達式 條件條件如：如：replace dum=0 if year=199026 在執(zhí)行在執(zhí)行 Stata 命令時，有時需要調(diào)用許多變量，而某些命令時，有時需要調(diào)用許多變量，而某些變量名可能很長。此時，如果在命令窗口一一輸入變量名，變量名可能很長。此時，如果在命令窗口一一輸入變量名，可能較費事?？赡茌^費事。 解決方法之一是：

27、可以直接在左下角的解決方法之一是：可以直接在左下角的“變量窗口變量窗口”單擊需要的變量，則該變量名就會顯現(xiàn)在命令窗口。單擊需要的變量，則該變量名就會顯現(xiàn)在命令窗口。 解決方法之二是，如果有以下變量解決方法之二是，如果有以下變量lnq1, lnq2, , lnq30，而只想使用其中的前而只想使用其中的前15 個變量，則可以用個變量，則可以用lnq1-lnq15 來簡來簡略地表示這略地表示這15 個變量。個變量。 解決方法之三是，用解決方法之三是，用“*”號來節(jié)省變量名的書寫。號來節(jié)省變量名的書寫。 假設(shè)想將內(nèi)存中所有以假設(shè)想將內(nèi)存中所有以“l(fā)n”開頭的變量都去掉，則可開頭的變量都去掉，則可輸入命

28、令輸入命令:. drop ln*277Stata 的計算器功能的計算器功能 Stata 也可以作為計算器來使用。只要輸入命令也可以作為計算器來使用。只要輸入命令“display expression”即可。即可。 比如比如“display log(2)”將計算將計算ln 2。288線性回歸分析線性回歸分析使用使用 OLS 對方程進行估計，可輸入命令：對方程進行估計，可輸入命令：reg lngov lngdp上表中的上表中的“_cons”表示常數(shù)項，表示常數(shù)項，“R-squared”顯示顯示R2 =0.9676，“Adj R-squared”顯示顯示0.9665。檢驗整個方程顯著性的。檢驗整個方

29、程顯著性的F 統(tǒng)計統(tǒng)計量之量之p 值（值（Prob F）為）為0.0000，顯示這個回歸方程是高度顯，顯示這個回歸方程是高度顯著的。殘差平方和為著的。殘差平方和為1.77，方程標準差，方程標準差Root MSE為為0.24。29如果要顯示估計系數(shù)的協(xié)方差矩陣，可輸入命令，如果要顯示估計系數(shù)的協(xié)方差矩陣，可輸入命令，. . vcevce其中，其中，“vce”vce”表示表示“variance covariance matrix variance covariance matrix estimated”estimated”。在進行回歸時，如果不要常數(shù)項，可以加上選擇項在進行回歸時，如果不要常數(shù)項，

30、可以加上選擇項“noconstant”noconstant”reg lngov lngdp,noconstantreg lngov lngdp,noconstant30 如果要計算被解釋變量的擬合值如果要計算被解釋變量的擬合值 ，并將其記為，并將其記為lntchat，可輸入命令：可輸入命令：. predict lntchaty 如果要計算如果要計算“殘差殘差”（residual），并將其記為），并將其記為e1，可輸入命令：可輸入命令：. predict e1,residual 其中，選擇項其中，選擇項“residual”表示預(yù)測殘差。表示預(yù)測殘差。319Stata 的日志的日志 如果希望在每次使用如果希望在每次使用 Stata 時，儲存其運行結(jié)果，時，儲存其運行結(jié)果，可點擊菜單可點擊菜單“File”“Log”“Begin”，然后輸入日志，然后輸入日志（log）的文件名，并存儲在你指定的位置。從此以后，）的文件名，并存儲在你指定的位置。從此以后，你在你在Stata中的所有操作及其輸出結(jié)果，都將被記錄在中的所有操作及其輸出結(jié)果，都將被記錄在此日志中，直至你