问题描述
我正在尝试在事件发生时的两个采样时间内在SAS中进行采样大小计算。
这里是情况:
- 假设两个样本都遵循指数分布
- 假设在替代假设下给定的恒定危险比,我们称为hr(第2组与第1组)
- 我们将使用logrank测试。
- 给出应计时间a,并跟进时间f
- 还给出了第1组的指数危害,称为exph1
- 假设两组之间的样本数量比为1:1
- 所需的标称功率为p
现在我的代码如下:
proc power;
twosamplesurvival test=logrank
accrualtime = a
followuptime = f
refsurvexphazard= exph1
hazardratio = hr
power = p
/* eventstotal = . /*events total */
/* ntotal= . /*total sample size */
;
run;
您可以取消注释eventstotal =。或ntotal =。取决于您是要计算请求的事件数还是实际的样本总数。
在随访结束时,它们不应被视为相同的对象,如果事件未发生,则主题将被正确审查。
但是,对于事件总数和样本总数,我总是得到相同的数字。我在这里做错了什么?
我实际上知道如何手工计算,并且我对请求的事件编号的手工计算非常接近SAS输出(SAS给出的值稍大但非常接近),但是我的总样本量比事件大得多数字。
由于机密原因,我无法透露上述参数的任何特定初始值。有人可以帮忙吗?真的会很感激。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)