使用蒙特卡羅模擬計算功率，*2部分:使用power運行模擬

時間：2019-12-25作者：北京天演融智軟件有限公司瀏覽：163

在上一篇文章中，向您展示了如何使用蒙特卡羅模擬計算t檢驗的功率。在本文中，將向您展示如何將模擬集成到Stata的power命令中，以便您可以輕松地為一系列參數(shù)值創(chuàng)建自定義表格和圖形。

統(tǒng)計學家在計劃一項科學研究時，很少為一組假設計算功率。我們通常計算一系列參數(shù)值的功率，并選擇一組經(jīng)濟上和邏輯上可行的現(xiàn)實假設。例如，下面我使用power onemean以10為增量計算50到100的樣本大小。該表顯示了假設的參數(shù)值，包括alpha水平、零值和替代假設下的平均值，平均值(delta)之間的標準化差異、標準偏差和每個樣本大小的功率。

還使用了下面的graph選項繪制了樣本大小范圍內(nèi)的功率。然后，使用表格和圖表來選擇一個滿足我的研究的功率要求的樣本大小。




除了樣本大小，power命令還允許您為其他參數(shù)輸入一系列值，例如標準偏差，平均值或alpha級別。并且power將創(chuàng)建結果的表格和圖表。

您還可以將自己的方法添加到power的強大功能套件中。讓我們將上一篇文章中的t測試模擬程序添加到power中，看看它是如何工作的。

回想上一篇文章，我們創(chuàng)建了一個名為simttest的程序來計算t檢驗的功效。該程序接受五個輸入?yún)?shù)，創(chuàng)建一個假隨機假設的偽隨機數(shù)據(jù)集，對虛假設進行檢驗，并返回假設檢驗的結果。

我們使用simulate多次運行程序并將結果保存到名為reject的變量中。

然后，根據(jù)零假設被拒絕的次數(shù)比例來計算功率。

您可以通過創(chuàng)建一個名為power_cmd_mymethod的程序?qū)⒋四M方法添加到power中，其中mymethod是power命令的名稱。讓我們調(diào)用程序power_cmd_simttest。




下面的代碼塊定義了power_cmd_simttest。 注意，它與我們的simttest程序非常相似。它從capture program drop開始，然后是program和version 15.1。 接下來，像在Simttest中那樣使用syntax定義輸入?yún)?shù)。這里我添加了一個名為reps（）的新輸入?yún)?shù)，它是模擬的重復次數(shù)。

程序的中間部分運行模擬并總結結果。在simulate和summarize之前都要使用quietly，這會抑制其輸出的顯示。 這里simulate運行程序simttest，并像以**樣將結果保存到變量reject中。請注意，simulate和simttest中的所有輸入?yún)?shù)都是使用syntax定義本地宏。summarize計算reject的平均值并將其存儲在標量r(mean)中。




代碼塊的底部部分返回功率和其他參數(shù)。 標量power返回變量reject的平均值，其他參數(shù)通過syntax傳遞的本地宏。




現(xiàn)在，可以通過輸入power simttest來運行模擬。

非常有效！甚至可以為一系列樣本大小制作表格和圖表。

如果您只想考慮一系列樣品尺寸范圍，可以在此處停止。 但是，如果希望要為其他參數(shù)（如m0，ma和sd）輸入一系列值，則需要再編寫一個小程序。 該程序必須命名為power_cmd_mymethod_init，因此我們將命名我們的程序power_cmd_simttest_init。




下面的代碼塊定義了power_cmd_simttest_init，并像其他程序一樣從從capture program drop和program開始。請注意，程序定義以sclass選項開頭。sreturn local pss_colnames行初始化輸出表中雙列引號中列出的參數(shù)的列。sreturn local pss_numopts行允許為雙引號中的參數(shù)*numlists。




現(xiàn)在，可以使用power simttest來計算假設不同替代假設的一系列均值的功率。甚至可以針對不同的樣本大小執(zhí)行此操作。




可以通過在graph選項中*xdimension繪制功率分析的結果。







到目前為止，我已經(jīng)展示了如何使用蒙特卡羅模擬計算功率以及如何將這些模擬集成到power中。 從一個簡單的t測試示例開始，以便我們可以專注于編程并使用power onemean檢查我們的工作。

更多stata軟件的相關信息請登錄科學軟件網(wǎng)








北京天演融智軟件有限公司專注于matlab,stata,spsspro,pscad,gams,nvivo,cyme,gms,mathematica等

• 詞條

  詞條說明

• 疫情期間，可申請遠程使用Stata

  由于新型冠狀病毒（COVID-19）在**各地蔓延，StataCorp LLC于2020年3月12日對于Stata軟件的高校用戶做出了相應的調(diào)整措施。新型冠狀病毒（COVID-19）的情況每天都在變化，各單位和高校都在考慮如何較好地**其員工和學生的健康和安全。許多大學正在將課程陸續(xù)過渡到在線平臺，很多單位也要求員工進行遠程工作。對于購買過Stata網(wǎng)絡版，學生實驗室版本或校園版許可的客戶，我們希

• Intel Fortran Compiler 語言編譯器

  Intel Fortran Compiler for Windows針對Windows的一個完整的Fortran開發(fā)環(huán)境? 兼容微軟Visual Studio 2013及以后版本。? 沒有Microsoft Visual Studio？也沒有沒問題：基于Microsoft Visual Studio 2015 shell的Fortran開發(fā)環(huán)境不需要購買其他任何產(chǎn)品。? 從熟悉的Microsoft

• EVS軟件CONNECTIVITY ASSESSMENT模塊介紹

  連接性評估應用連通性評估可有效應用于:? 分析(如化學數(shù)據(jù))?√選擇性地移除錯誤區(qū)域，例如投影到數(shù)據(jù)范圍之外的羽流碎片。??3D巖性模型√明確評估材料的不連續(xù)區(qū)域?√鎖定有足夠儲量的礦體進行開采??? 具有導致不連續(xù)區(qū)域（塊）的夾持層的地層模型。? ?√逐層厚度細分后，去除微小的層碎片。?利用分析數(shù)據(jù)進行連通性評估當一

• 在 MODFLOW 6 中使用交換功能

  在邊界條件不能充分描述彼此之間關系的情況下，如何讓模型在?MODFLOW?6?中平穩(wěn)運行？?在?MODFLOW?6?中，您可以選擇通過交換功能將模型相互關聯(lián)，而不僅僅是使用邊界條件。這就允許?MODFLOW?6?計算每個模型之間的流量，就像它是一個大型非結構化網(wǎng)格的一部分。通過交換連接模型，允許在具