使用蒙特卡罗模拟计算功率，*2部分:使用power运行模拟

在上一篇文章中，向您展示了如何使用蒙特卡罗模拟计算t检验的功率。在本文中，将向您展示如何将模拟集成到Stata的power命令中，以便您可以轻松地为一系列参数值创建自定义表格和图形。

统计学家在计划一项科学研究时，很少为一组假设计算功率。我们通常计算一系列参数值的功率，并选择一组经济上和逻辑上可行的现实假设。例如，下面我使用power onemean以10为增量计算50到100的样本大小。该表显示了假设的参数值，包括alpha水平、零值和替代假设下的平均值，平均值(delta)之间的标准化差异、标准偏差和每个样本大小的功率。

还使用了下面的graph选项绘制了样本大小范围内的功率。然后，使用表格和图表来选择一个满足我的研究的功率要求的样本大小。

除了样本大小，power命令还允许您为其他参数输入一系列值，例如标准偏差，平均值或alpha级别。并且power将创建结果的表格和图表。

您还可以将自己的方法添加到power的强大功能套件中。让我们将上一篇文章中的t测试模拟程序添加到power中，看看它是如何工作的。

回想上一篇文章，我们创建了一个名为simttest的程序来计算t检验的功效。该程序接受五个输入参数，创建一个假随机假设的伪随机数据集，对虚假设进行检验，并返回假设检验的结果。

我们使用simulate多次运行程序并将结果保存到名为reject的变量中。

然后，根据零假设被拒绝的次数比例来计算功率。

您可以通过创建一个名为power_cmd_mymethod的程序将此模拟方法添加到power中，其中mymethod是power命令的名称。让我们调用程序power_cmd_simttest。

下面的代码块定义了power_cmd_simttest。 注意，它与我们的simttest程序非常相似。它从capture program drop开始，然后是program和version 15.1。 接下来，像在Simttest中那样使用syntax定义输入参数。这里我添加了一个名为reps（）的新输入参数，它是模拟的重复次数。

程序的中间部分运行模拟并总结结果。在simulate和summarize之前都要使用quietly，这会抑制其输出的显示。 这里simulate运行程序simttest，并像以**样将结果保存到变量reject中。请注意，simulate和simttest中的所有输入参数都是使用syntax定义本地宏。summarize计算reject的平均值并将其存储在标量r(mean)中。

代码块的底部部分返回功率和其他参数。 标量power返回变量reject的平均值，其他参数通过syntax传递的本地宏。

现在，可以通过输入power simttest来运行模拟。

非常有效！甚至可以为一系列样本大小制作表格和图表。

如果您只想考虑一系列样品尺寸范围，可以在此处停止。 但是，如果希望要为其他参数（如m0，ma和sd）输入一系列值，则需要再编写一个小程序。 该程序必须命名为power_cmd_mymethod_init，因此我们将命名我们的程序power_cmd_simttest_init。

下面的代码块定义了power_cmd_simttest_init，并像其他程序一样从从capture program drop和program开始。请注意，程序定义以sclass选项开头。sreturn local pss_colnames行初始化输出表中双列引号中列出的参数的列。sreturn local pss_numopts行允许为双引号中的参数*numlists。

现在，可以使用power simttest来计算假设不同替代假设的一系列均值的功率。甚至可以针对不同的样本大小执行此操作。

可以通过在graph选项中*xdimension绘制功率分析的结果。

到目前为止，我已经展示了如何使用蒙特卡罗模拟计算功率以及如何将这些模拟集成到power中。 从一个简单的t测试示例开始，以便我们可以专注于编程并使用power onemean检查我们的工作。

更多stata软件的相关信息请登录科学软件网