我們很高興地向您介紹Sata16的新功能。下面，我們列出了該版本的所有亮點。接下來會介紹更多關于前13個功能。我們介紹每種功能所使用的詞語，您也可以在向現有和潛在的Stata用戶介紹它們時使用。
大多數這些功能都將會使所有學科領域的研究員們感到興奮。在適當的情況下，我們將強調哪些學科領域是人們最感興趣的或者提供不同領域的用戶如何與該功能相關的建議。我們將發送（單獨）列表，按照學科對這些功能進行排序。
1. Lasso
2. Reporting
3. Meta-analysis
4. Choice models
5. Python integration
6. New in Bayesian analysis—Multiple chains, predictions, and more
7. Panel-data ERMs
8. Import data from SAS and SPSS
9. Nonparametric series regression
10. Multiple datasets in memory
11. Sample-size analysis for confidence intervals
12. Nonlinear DSGE models
13. Multiple-group IRT models
14. xtheckman
15. Multiple-dose pharmacokinetic modeling
16. Heteroskedastic ordered probit models
17. Graph sizes in printer points, centimeters, and inches
18. Numerical integration
19. Linear programming
20. Stata in Korean
21. Mac interface now supports Dark Mode and native tabbed windows
22. Do-file Editor—Autocompletion and more syntax highlighting
1. Lasso
Lasso是一種用于模型選擇、預測和推斷的機器學習技術。
新的lasso命令使用線性、泊松、logit或probit回歸模型的偏差為連續、計數和二進制結果選擇“最佳”預測器。
例如，如果輸入
. lasso linear y x1-x500
lasso將選擇指定協變量的子集-例如x2，x10，x11和x21。然后，您可以使用標準的predict命令獲得y的預測值。
如果您有一個二進制或計數的結果，您可以用同樣的方法使用lasso logit、lasso probit或lasso poisson。如果您更喜歡使用彈性網絡或平方根lasso方法來選擇變量，那么可以使用elasticnet或sqrtlasso命令。
有時，變量選擇或預測是您使用lasso的最終目標。而其他時候，您可能對估算和檢驗系數感興趣。Stata 16提供了11個命令，允許您參數估計、標準誤差和置信區間，并在使用lasso方法從潛在控制變量中選擇的同時，對感興趣的變量執行檢驗。這些命令是：dsregress、dslogit、dspoisson、poregress、pologit、popoisson、poivpoisson、xporegress、xpologit、xpopoisson以及xpoivregress.
ds系列命令執行double-selection lasso，po系列命令執行partialing-out lasso，以及xpo系列命令執行cross-fit partialing-out lasso。他們對具有連續、二進制和計數結果的模型執行此操作。它們甚至可以在模型中處理內生協變量，以獲得連續的結果。目前文獻討論了許多基于lasso的推斷方法。我們提供了一些這樣的方法，以便研究人員能夠選擇他們最喜歡的。事實上，文獻中還有更多基于lasso的推斷方法，研究人員通?？梢允褂?strong>lasso、sqrtlasso和elasticnet中提供的工具來實現其他方法。
lasso和elasticnet命令是用于變量選擇和預測的標準lasso工具。lasso推斷工具實現了主要由計量經濟學家開發的新方法。然而，這些推斷方法將在所有學科中流行，因為它們提供了一種檢驗和解釋目標變量系數的方法。
用戶可以在新的lasso參考手冊中輕松了解lasso的所有功能。
2.報告功能
Stata的報告功能允許您創建Word、PDF、Excel和HTML文檔，這些文檔將Stata結果和圖形與格式化的文本和表格結合在一起。無論您創建的文檔類型如何，都可以依靠Stata的集成版本控制功能來確保您的報告是可復制的。
想要在數據更改時同步更新動態報告嗎？Stata的報告功能使這一點變得簡單。使用更新的數據集重新運行創建報表的命令或do文件，并自動更新報告中所有Stata結果。
當然，Stata16具有新的和改進的報告功能，但同樣重要的是，Stata的所有報告功能現在都記錄在新的報告參考手冊中。本手冊包括許多新的示例，用于演示了工作流程，并提供了有關使用Stata自定義創建Word、PDF、Excel和HTML文檔的指導。
Stata 16中的新報告功能：
l dyndoc和markdown命令現在除了以前創建的HTML文檔之外，還創建了Word文檔?，F在，您可以輕松地將所有的Stata輸出和帶有Markdown格式文本的圖形合并，以創建定制的Word文檔
l Do文件編輯器現在為Markdown語言元素提供語法突出顯示。
l putdocx命令現在允許您包含頁眉、頁腳和頁碼。它使大文本塊的編寫更加容易。
l html2docx命令將HTML文檔（包括css）轉換為Word文檔。
l Docx2pdf命令將word文檔轉換為PDF格式。
3. Meta分析
Stata16有一套執行meta分析的新命令。這一套命令可以讓您探索和組合不同研究的結果。例如，如果您收集了20項關于特定藥物對血壓影響的研究的結果，您可以總結這些研究，并使用Meta分析估計總體效果。
新的Meta分析組件很廣泛，但它的獨特之處在于它的簡潔性。
例如，您可以輸入
. meta set effectsize stderr
聲明預先計算的效果大小，或使用meta esize從匯總數據來計算效果。使用這種方法，您可以執行隨機效應、固定效應或共同效應meta分析。
為了估算整體效果大小及其置信區間，獲取異質性統計數據等等，只需輸入：
. meta summarize
且可視化結果就像打字一樣簡單：
. meta forestplot
但是meta套件提供了更多功能。
Meta回歸分析和亞組分析可以讓您評估研究的異質性。這些可以通過meta regress和meta forestplot, subgroup()或者meta summarize, subgroup()獲得。
您可以調查潛在的發布偏差。使用meta funnelplot從外觀上檢查漏斗圖的不對稱性；使用meta bias正式地檢驗漏斗圖的不對稱性；使用meta trimfill的剪補法評估發布的偏差。
您甚至可以使用meta summary，cumulative（）執行累積meta分析。
所有的meta分析功能都記錄在新的meta分析參考手冊中。
4. 選擇模型
在Stata16中，我們引入了一套新的統一的命令套件，用于建模選擇數據。添加了用于匯總選擇數據的新命令。重命名并改進了用于擬合選擇模型的現有命令。我們甚至添加了一個新的命令，用于面板數據擬合混合Logit模型。我們在新的“選擇模型參考手冊”中將它們一起記錄下來。
這里是最好的部分:在選擇合適的模型后，margins起了作用。這意味著您現在可以輕松地解釋所選擇模型的結果。雖然選擇模型中估計的系數通常幾乎無法解釋，但是margins允許您根據結果詢問和回答非常具體的問題。假設您在模擬交通的選擇，您可以回答以下問題：
?預計有多少比例的旅客選擇航空旅行？
?每增加1萬美元的收入，乘汽車旅行的概率如何變化？
?如果機場的等待時間增加30分鐘，這對每種交通方式的選擇有何影響？
還有什么新的部分呢？現在，在擬合選擇模型之前，先cmset數據。例如：
.cmset personid transportmethod
然后，使用cmsummary、cmchoiceset、cmtab和cmsample來探索、匯總和查找您數據中潛在的問題。
您可以使用cm系列的估計命令來選擇以下選擇模型之一：
?cmclogit 條件logit（McFadden選擇）模型
?cmmixlogit 混合logit模型
?cmxtmixlogit 面板數據混合logit模型
?cmmprobit 多項式Probit模型
?cmroprobit 排序Probit模型
?cmrologit 排序logit模型
與其他命令不同，cmxtmixlogit不僅被重命名和改進。它在Stata16中也是全新的，并且適用于面板數據的混合logit模型。
5. Python集成
在stata 16中，可以從Stata中嵌入并執行python代碼。Stata新的python命令允許您從Stata輕松調用python，并在Stata中輸出python結果。
您可以交互地或在do文件和ado文件中調用python，以便利用python的廣泛語言特性。還可以直接通過Stata執行python腳本文件（.py）。
此外，我們還介紹了Stata函數接口（sfi）python模塊，它提供了Stata和python之間的雙向連接。此模塊允許您訪問Stata的當前數據集、框架、宏、標量、矩陣、值標簽、特征、全局Mata矩陣等。
所有這些都意味著您現在可以直接在Stata中使用任何python包。例如，可以使用matplotlib繪制三維圖形?？梢允褂胣umpy進行數值計算。您可以使用scrappy從網頁上抓取數據?？梢酝ㄟ^TensorFlow和Scikit Learn訪問其他機器學習技術，如神經網絡和支持向量機；以及更多的技術。
最后，Stata的do文件編輯器現在包含了Python語言的語法突出顯示。
雖然高級用戶和程序員最有可能使用Python集成，但Stata中Python的可用性將激發所有學科領域中的更多的用戶。
6.貝葉斯分析中的新功能—多鏈，預測等等）
多鏈條：
基于MCMC (Markov chain Monte Carlo)樣本的貝葉斯推斷只有在Markov鏈收斂時才有效。評估這種收斂性的一種方法是模擬和比較多個鏈。
新的nchans（）選項可以與bayes:prefix和bayesmh命令一起使用。例如，您可以輸入：
. bayes, nchains(4): regress y x1 x2
將生產四條鏈。這些鏈將結合在一起產生更精確的最終結果。然而，在解釋結果之前，可以通過圖形比較鏈來評估收斂。還可以使用現在由Bayes報告的Gelman-Rubin收斂診斷來評估收斂性：當模擬多個鏈時，使用的regress和其他貝葉斯估計命令。當您關心非收斂性時，可以使用bayesstats grubin命令進一步研究，以獲得模型中每個參數的單個gelman–rubin診斷。
貝葉斯預測：
貝葉斯預測是來自后驗預測分布的模擬值。這些預測對于檢驗模型擬合和預測樣本外觀測很有用。使用bayesmh擬合模型后，可以使用bayespredict來計算這些模擬值或函數，并將它們保存在新的Stata數據集中。例如，您可以輸入：
. bayespredict (ymin:@min({_ysim})) (ymax:@max({_ysim})), saving(yminmax)
來計算模擬值的最小值和最大值。然后，可以使用其他的后置估計命令（如bayesgraph）來獲得預測的摘要信息。bayespredict創建的數據集可能包含數據集中每個觀測的數千個模擬值。有時，您不需要所有這些單獨的值。為了獲得后驗，如后驗均值或中位數，您可以使用bayespredict、pmean或bayespredict、pmedian?；蛘?，您可能對模擬值的隨機樣本感興趣。例如，您可以使用bayesreps、nreps（100）來獲得100個副本。
最后，您可能希望使用后驗預測p值，也稱為PPPs或貝葉斯預測p值來評估模型的擬合優度。PPPs測量觀測數據和復制數據之間的協議，可以使用最新的bayesstats ppvalues命令計算。例如，使用我們前面的示例：
. bayesstats ppvalues {ymin} {ymax} using yminmax
7.面板數據的ERMs
擴展回歸模型（ERMs）是上次版本的一個新功能。ERM命令適用于解釋觀測數據中三個常見問題的模型：內生協變量、樣本選擇和單獨或聯合處理。
在Stata16中，我們介紹了用于擬合面板數據ERM模型的xtererss、xteintreg、xteprobit和xteoprobit命令。這意味著ERM現在可以解釋上面提到的三個問題以及面板內相關性。這些新命令適用于隨機效應線性、區間、概率和有序概率回歸模型。它們允許一個或所有方程中的隨機效應，并且允許隨機效應在方程間相互關聯。
使用觀測（非實驗）數據的所有學科的研究人員都對ERM模型感興趣，并對這些命令的新面板數據版本感到興奮。然而，不同的學科領域對這些模型的討論也不同。
上面，我們提到ERMs模型解決的問題稱為內生協變量、樣本選擇、處理以及面板內相關。雖然這個術語在經濟學等一些學科中很常見，但其他學科可能會使用其他術語。
?研究人員可能會要求多層次（兩層）數據模型來解釋組內相關性，而不是面板數據和面板內相關性。
?研究人員可能會要求處理未觀察到的或未測量到的干擾因子，而不是內生協變量。
?與樣本選擇不同，研究人員可能關注的是信息丟失、不可忽視的無反應或隨機丟失的結果（MNAR）的試驗。
?研究人員可能會詢問因果推斷或估計平均處理效果（ATEs）的方法，而不是治療方法。
值得注意的是，盡管所有的學科領域都對ERM模型感興趣，但它們通常使用不同的語言。
8.從SAS和SPSS中導入數據
通過使用Stata16新的導入SAS和SPSS的命令，您現在可以導入存儲成SAS（.sas7bdat）和SPSS（.sav）格式的數據。通過對話框可以在導入數據之前輕松瀏覽數據，并根據需要選擇要加載到Stata中的變量和觀察的子集。例如，
此外，使用新的import saxport8和export saxport8命令，可以將SAS XPORT版本8的傳輸文件導入和導出到Stata中?，F有的import saxport和export saxport命令與SAS XPORT版本5傳輸文件一起使用，并已重命名為import saxport5和export saxport5。
9.非參數序列回歸
Stata16中新的npregress系列命令適用于使用多項式、B樣條或協變量樣條近似因變量平均值的非參數序列回歸模型。這意味著您不需要指定任何預定的功能形式。只需指定您希望在模型中包含哪些協變量。例如，鍵入
. npregress series wineoutput rainfall temperature i.irrigation
npregress系列報告的不是系數，而是報告結果，連續變量的平均邊際效應和分類變量的對比。結果表明，降雨的平均邊際效應為1，灌溉的對比度為2。這種對比可以解釋為灌溉的平均處理效果。
作為非參數回歸，未知均值近似為協變量的一系列函數。然而，我們仍然可以從參數模型中得到推論。我們只使用margins。我們可以輸入：
. margins irrigation, at(temperature=(40(5)90))
并得到在40、50、…、90度的溫度下進行灌溉的預期效果表。我們可以用marginsplot來繪制結果。
更重要的是，npregress系列可以適用于部分參數（半參數）模型。
10.內存中的多個數據集
現在可以將多個數據集加載到內存中。輸入：
. use people
然后people.dta被載入內存。接下來，輸入：
. frame create counties
. frame counties: use counties
內存中有兩個數據集。people.dta位于名為default的框架中，countries.dta位于名為countries的框架中。當前的框架仍為default。大多數Stata命令使用當前框架中的數據。例如，如果輸入：
. list
會列出people.dta，如果輸入：
.frame counties: list
然后會列出counties.dta?；蛘吣梢酝ㄟ^鍵入使counties成為當前框架：
.frame change counties
現在list就會列出counties數據。
導航框架很容易，連接它們也很容易。假設兩個數據集都有一個名為countycode的變量，該變量以相同的方式標識縣。鍵入：
. frlink m:1 countycode, frame(counties)
默認框架中的每個人都會連接到counties框架中的一個縣。這意味著您現在可以使用frget命令將變量從countries框架復制到當前框架?；蛘呖梢允褂?strong>frval（）函數直接訪問countries框架中變量的值。例如，如果我們將每個人的收入設置到default框架中，在counties框架中有縣收入的中位數，那么我們可以通過鍵入來生成一個包含相對收入的新變量：
. generate rel_income = income / frval(counties, median_income)
這只是開始。雖然本例中僅使用兩個框架，但一次最多可以在內存中有100個框架，并且這些框架之間可以有許多鏈接。
11.置信區間的樣本量分析
新的ciwidth命令執行精度和樣本大?。≒rSS）分析，即置信區間（CIs）的樣本大小分析。當計劃研究時使用此方法，并且希望在使用CIs進行推斷時以最佳方式分配資源。換言之，當您想要估計在計劃研究中獲得所需的CI精度所需的樣本量時，可以使用此方法。
ciwidth生成的樣本大小、精度以及：
?一個平均值的CI
?一個方差的CI
?兩個獨立均值的CI
?兩個配對均的CI
控制面板界面允許您選擇分析類型和輸入假設以獲得所需的結果。
ciwidth允許在可自定義的表和圖形中顯示結果。
ciwidth還為您提供了添加自己方法的工具。
12.非線性DSGE模型
在Stata15中，我們引入了用于擬合線性DSGE模型的dsge命令，這是模型用于經濟學和金融學的時間序列模型。這些模型是傳統預測模型的替代方案。這兩個模型都試圖解釋總體經濟現象，但DSGE模型是基于微觀經濟理論推導出的模型。
在Stata 16的更新中，dsgenl命令適用于非線性DSGE 模型。大多數DSGE模型都是非線性的，這意味著您不再需要手動行線性化它們。當您將方程輸入dsgenl時，它會為您將它們線性化。
用dsgenl估算模型的參數后，可以得到過渡矩陣和政策矩陣；確定模型的穩態；估計方程組所隱含的變量的方差、協方差和自方差；創建并繪制脈沖響應函數。
這可能是宏觀經濟學家和任何在央行工作的人最喜歡的特征。
13.多組IRT模型
IRT模型探討了潛在（未被觀察到的）特征與測量該特征各方面的項目之間的關系。這通常發生在標準化檢驗中，其中感興趣的特征是能力，例如數學能力。IRT模型設計了一組項目（測試問題），并對這些未被觀察到的特征進行了評估。教育、心理學和健康領域的研究人員經常采用IRT模型。
Stata的irt命令適合1、2和3參數邏輯模型。它們還適合分級響應、名義響應、部分信用和評級尺度模型及它們的任何組合。在擬合模型后，irtgraph繪制項目特征曲線、測試特征曲線、項目信息函數和測試信息函數。
在Stata 16的更新中，irt命令允許跨組比較。采用任何現有的irt命令，添加一個group（varname）選項，并適應相應的多組模型。例如，鍵入
.irt 2pl item1-item10, group(female)
便匹配了兩組2PL模型
IRT模型對潛在性狀的群體特異性均值和方差進行估計。還可以估計一個或多個項目項目組難度和識別參數。通過約束，您可以精確地指定哪些參數允許更改，哪些參數在組之間被約束為相等。
您甚至可以使用似然比測試來比較有約束和沒有約束的模型，以執行一個基于IRT模型的差異項功能測試。
還有更多！
