sas计算rr用啥

Posted 2023-05-08

客观缓解率(Objective Response Rate ,ORR)（有的称总体缓解率）：指肿瘤体积缩小达到预先规定值并能维持最低时限要求的患者比例。一般定义客观缓解率为完全缓解加上部分缓解之和。

总体缓解率（ORR）：接受治疗后肿瘤评估至少一次达到完全缓解（CR）或者部分缓解（PR）的患者的比例。

疾病控制率（DCR）：在接受治疗后肿瘤评估至少一次达到完全缓解（CR）、部分缓解（PR）或病变稳定（SD）的患者的比例。

3.疗效指标计算
3.1时间-事件类型
3.1.1 SAS计算：
计算不同组别疾病进展时间（TTP）和无进展生存期（PFS），采用Kaplan-Meier法进行描述性统计，采用Greenwood方法计算标准误差，并提供95%置信区间，绘制生存曲线图。采用log-rank检验比较不同治疗组之间生存曲线是否有差异。

部分原始数据：

原始数据
可私信我发送原版Excel数据，也使用R语言去下载原lung数据。

3.1.2 SAS code：
*
输入：
Lung为数据集（包含subjid受试者、time生存时间、status状态）
每个生存时间生存率计算=本时间点死亡数量/本时间点开始时存活数量。
本时间点删失的不算死亡，但是计算下一个时间点时算死亡，不计入下个时间点计算的分母中。

输出：
1.采用ci表画生存曲线图。
2.lpgrank检验的p值在等价与组检验表单。
3.中位生存时间和CI在lung_3。;

proc lifetest data=lung outsurv=ci method=km atrisk conftype=log;
time time*status(1);
strata /group=sex test=(logrank);
ods output productLimitEstimates=lung_2 quartiles=lung_3;
run;

*code解释：
proc lifetest:调用生存分析过程步。
data=lung：输入lung数据集。
outsurv=ci：输出ci数据表（见如下截图）。
method=km：绘制Kaplan-Meier曲线。
atrisk ：输出lung_2的NumberAtRisk有风险人数（分母），ObservedEvents达到终点事件人数。
time*status(1)：生存时间*状态（删失标识）。
test=(logrank)：采用logrank检验法。
ods output：输出两个数据集lung_2，lung_3（见下文。）
;
CI表
Lung_2表
Lung_3表
生存曲线图
3.1.2 计算原理如下：

设t_1< t_2<t_3<t_4<....<t_k,：表示事件发生的k个时间点；

n_j：表示时刻tj前的生存人数；

d_j：表示t_j时刻发生事件（到达终点）的人数，并记s_j=n_j-d_j;

则生存分布函数（SDF）如下：

hatSleft ( t_j right )=prod_i=1^j(1-fracd_kn_k)
标准误为（Greenwood法）：

hatsigma( hatSleft ( t_j right ))=hatSleft ( t_j right )sqrtsum_k=1^jfracd_kn_ks_k
95%置信区间为，第一种方式（Kalbfleisch&Prentice 1980.,excel暂时未计算出来）：

hatSleft ( t_j right )^exp(pm 1.96s)\ s=hatsigma( hatSleft ( t_j right ))/(hatSleft ( t_j right )*ln(hatSleft ( t_j right )))
95%置信区间为，第二种方式（SAS计算方式）：

hatSleft ( t_j right ) pm 1.96 hatsigma (hatSleft ( t_j right ))
可私信我获取excel计算文件。

两条生存曲线的检验（logrank法）:

H0：两组个体在任意时间点，生存率无差异。

H1：两组个体在任意时间点，生存率有差异。

组别

到达终点事件数

生存个体数

风险个体数(合计)

组A

daj

saj=naj-daj

naj

组B

dbj

sbj=nbj-dbj

nbj

合计

dj

sj

nj

将所有的时间点都罗列出来，计算每个时间点的实际频数和理论频数，构建x^2 统计量

sum_i=1^k(Q_A-E_A)/E_A+sum_i=1^k(Q_B-E_B)/E_B

即可计算出p值。可参考https://zhuanlan.zhihu.com/p/392104512?ivk_sa=1024320u。

风险比率（hazard ratio,HR）的计算：

HR=fracQ_A/E_AQ_B/E_B

3.2率指标sas计算
计算各组ORR、DCR及95%置信区间，采用Fisher精确检验进行比较。（若真实世界研究，则同时采用倾向性评分匹配的方法，对基线进行矫正后采用Fisher精确检验进行比较。）

sas代码如下:

*
输入：
rr为数据集（包含SUBJID受试者编号、ARM组别、ORR是否达到ORR）

输出：
OneWayFreqs:单因子频数
Binomial:二项式比例
BinomialTest:二项式比例检验
单组率和95CI在二项式比例表中。

;

ods output Binomial= Binomial;
proc freq data=rr;
by arm;
table ORR/binomial alpha=0.05;
run;
4.关于删失的一些考虑
事件/删失

事件/删失原因

分类

事件

进展

进展

事件

死亡

死亡

事件

没有基线影像学检查

没有基线影像学检查

删失

发生事件前开始了新的抗肿瘤治疗

开始新的抗肿瘤治疗

删失

发生事件的日期距上一次影像学检查日期大于两个评估周期

连续两次影像学检查缺失

删失

服药后无疗效评价或疗效评价的结果均为NE

服药后无疗效评价

删失

没有发生事件，且受试者仍在继续服药

继续治疗

删失

没有发生事件，且受试者已出组

失访

在计算生存率时对删失的处理：本时间点删失的不算死亡，但是计算下一个时间点时算死亡，不计入下个时间点计算的分母中。具体可看excel表格中对生存率的计算。

5.参考文献
1.CDE.抗肿瘤药物临床试验终点技术指导原则.

2.CDE.抗肿瘤药物上市申请临床数据收集技术指导原则.

3.CDE.抗肿瘤药物临床试验技术指导原则.

4.渣渣东的微信公众号.https://mp.weixin.qq.com/s/U8lxK9EmitQmb3wr_FNVRA.

5.陈峰等.临床试验统计学.

6.知乎王不留行.https://zhuanlan.zhihu.com/p/392104512?ivk_sa=1024320u. 参考技术A SAS计算RR的方法是通过计算受试者组（treatment group）和对照组（control group）之间的比值来计算的。RR即为受试者组和对照组之间的比值，又称为风险比（risk ratio）或相对风险（relative risk）。RR的计算公式为：RR = P(受试者组患病) / P(对照组患病) ，即计算受试者组患病的概率除以对照组患病的概率，得到受试者组和对照组之间患病的风险比。RR大于1，表明受试者组患病的风险比对照组要高。RR小于1，表明受试者组患病的风险比对照组要低。RR等于1，表明受试者组和对照组患病的风险比一样。SAS计算RR的方法也可以使用统计模型，如回归分析，比如Cox比例危险模型，以探究影响结果的因素可能性。 参考技术B 在SAS中，计算风险比（RR）通常需要使用PROC FREQ和PROC MEANS两个过程来完成。

首先，使用PROC FREQ过程来计算出各组别的发生次数，然后计算出相应的事件概率或者危险概率。接下来，使用PROC MEANS过程计算出每组别的均值以及标准误差，进而计算RR。

具体而言，以下是计算RR的步骤：

1. 使用PROC FREQ过程计算出各组别的发生次数
2. 根据发生次数计算出每组别的事件概率或危险概率
3. 使用PROC MEANS过程计算出每组别的均值和标准误差
4. 计算出每组别之间的RR，公式为RR = (事件概率1 / 事件概率2)

总之，通过使用SAS中的PROC FREQ和PROC MEANS过程，可以比较容易地计算出RR，从而评估不同组别之间的风险差异

时间片轮转算法（RR）能不能用在作业调度上?还是只能用在进程调度上？

导师给出的答案是不能，但是吧友似乎有一些别的看法。
先来看作业调度的定义，高级调度又称长程调度或作业调度，他的调度对象是作业。主要功能是根据算法，决定将外存处于后备队列的哪几个作业调入内存，为他们创建进程，分配必要的资源，并将它们放入就绪队列。分时和实时系统不设置高级调度。
其次简单的叙述一下低级调度（进程调度）：对象为进程，功能是决定就绪队列中的哪个进程应获得处理机，并分配。
最后来看一下RR的核心，进程切换：在RR调度算法中，应在合适进行进程的切换，可分为两种情况：1，若一个时间片尚未用完，正在运行的进程便已经完成，立刻激活调度程序，将他从就绪队列中删除，再调度就绪队列中对手的进程运行，并启动一个新的时间片。2.在一个时间片用完时，计时调度器中断处理程序被激活。如果进程尚未运行完毕，调度程序将把它送往就绪队列的末尾。
可见，RR算法在切换进程时只涉及了就绪队列的c、a、o作，也就是只有低级调度参与，并没有从外存将作业调入内存，另外，RR算法常用于分时系统，分时系统也是不设置作业调度的。 参考技术A 我也想知道这个问题 参考技术B 不能用在作业调度上
我就知道这个，并且我也不知道为啥。