发帖
雷达卡
目前常用的样本量计算软件PASS虽然能够计算COX回归的样本量，但是从下图截图中可以看到其中一个重要参数 R-Squared of X1 with Other X’s是比较难以获得的。
查询帮助，关于该参数的定义为：
This is the R-Squared that is obtained when X1 is regressed on the other X’s (covariates) in the model.
也即COX模型中其他协变量对X1的决定系数。对于连续性变量而言，该系数还是比较容易计算得到的，但是如果X1为分类变量，X1的决定系数似乎很难理解。查阅该方法的原始文献,分别是以下两篇文章：
Hsieh, F.Y. and Lavori, P.W. 2000. 'Sample-Size Calculations for the Cox Proportional Hazards Regression Model with Nonbinary Covariates', Controlled Clinical Trials, Volume 21, pages 552-560.（下载网址：https://yunpan.cn/cY8GWH5AVvPyb 访问密码 d6f5）
Schoenfeld, David A. 1983. 'Sample-Size Formula for the Proportional-Hazards Regression Model', Biometrics, Volume 39, pages 499-503.（下载网址：https://yunpan.cn/cY8G5WrtjVmnz 访问密码 46e2）
从第一篇文章可以看到，PASS计算的COX回归样本量很可能只适用于协变量为非二分类的连续性变量，相应的方法也见于R（来源：http://artax.karlin.mff.cuni.cz/ ... piCont.default.html）。
但是从R中其他类似的COX回归样本量计算方法中可以看到，似乎该决定系数（R中命名为rho2） 也存在相应的计算方法：
来源：http://artax.karlin.mff.cuni.cz/ ... /html/ssizeEpi.html
ρ=corr(X_1, X_2)=(p_1-p_0)\times√{\frac{q(1-q)}{p(1-p)}},
其中 p=Pr(X_1=1), q=Pr(X_2=1), p_0=Pr(X_1=1|X_2=0), and p_1=Pr(X_1=1 | X_2=1).
上述公式显示似乎并不清楚，查阅原始文献Sample size formula for proportional hazards modelling of competing risks（下载网址：https://yunpan.cn/cYWEXIEPfy23z 访问密码 09c4），公式及说明如下：
但是该公式似乎只对COX模型中存在两个2分类的协变量时有效，对于存在3个及以上的2分类协变量，以及如果同时存在2分类的协变量以及连续性变量时，似乎目前还没有很好的方法。所以大家在采用软件计算COX回归的样本量时应该注意该方法的适用范围，然而软件本身并没有就这点详细说明。
如果大家知道有针对上述局限的样本量计算方法，欢迎交流！
京公网安备 11010802022788号
