大家好
我们正在计划一项事件发生时间分析研究。我们的目标是比较两组并评估事件发生时间是否存在显著差异。根据我们的初步估计,我们预计两组之间可能没有显著差异。然而,为了确保研究的有效性,我们希望为每个组选择适当的样本量,这样如果确实存在差异,我们就可以自信地拒绝零假设。
我们打算使用 SAS JMP 软件来计算所需的样本量,但我们不确定此过程的具体细节。有人可以提供指导或向我们提供资源,说明如何使用 SAS JMP 计算事件发生时间研究的样本量吗?
谢谢你!
感谢马克的回复:)
我已经尝试过这个,但我不确定我从这个计算中得到的样本量是否可以用于比较研究?!
谢谢!
嘿@萨巴,
不幸的是,正如您所怀疑的那样,可靠性寿命探索器不太适合在比较两组时确定样本量。目前还没有一个平台可以直接做到这一点,所以您需要查看我们的愿望清单以提出建议!
您可以采取的一种解决方法(假设您有 JMP Pro)是利用许多平台中提供的一键模拟功能。首先创建一个数据表,其中包含一个组列、一个事件发生时间列、一个审查列,然后是一个用于模拟事件发生时间的列。以下是该数据表的示例:

对于事件发生时间列,我简单地粘贴了模拟时间列中的值。模拟列中引用了表变量(不用担心;我会在我的回复中分享此表,以便您试用)。
接下来,我们可以在 Analyze>>Reliability and Survival 菜单下运行 Fit Life by X。像这样填写:

单击“运行”,然后向下滚动,直到看到标有“测试结果”的表:

右键单击标有 Prob>ChiSq 的列并选择 Simulate(它位于弹出菜单的底部附近)。您将看到如下窗口:

选择切换事件时间和切换模拟时间。确保审查指标列引用模拟列,否则由此生成的模拟时间将不正确。
将样本数量设置为任意数字,以及随机数(我使用了 12345),然后单击“确定”。模拟完成后,您将获得如下所示的数据表:

为了节省时间,我决定提前停止(大约 1000 次运行)。其中一个表脚本标记为功效分析。如果您从具有 p 值的列运行一键模拟,JMP 会自动假定您可能正在进行自定义功效分析。如果您运行该脚本,它将打开分布平台的一个特殊实例:

在这种情况下,对于 alpha 为 0.05 的情况,我们有大约 30% 的功效(准确地说是 %29.577)来检测我在表中指出的差异(100 个时间单位 vs. 150 个时间单位)。此时,我可以返回到表格,为每个组添加更多行(最简单的方法可能是滚动到该组的最后一行,右键单击,选择“添加行”,添加所需的行数,然后填写组号),然后重复整个过程。
我承认这很粗糙而且耗时,但这是一个开始。至少,你可以得到一个不错的样本大小范围以供考虑。不过一定要看看愿望清单并为此添加建议!
嗨,卡金
感谢您的回复
我认为这会很有帮助。我对你解释的这部分有一个疑问,“确保审查指标列引用模拟列,否则它将不正确,从而生成模拟时间。”我应该如何确定这一点?
谢谢你!