Name: 竞争风险模型的基于 R 的景观验证 (视频) | JoVE Journal
Uploaded: 2022-09-16T13:10:00.000+00:00
Duration: 5 min 37 s
Description: 学习使用R语言评估竞争风险模型的判别与校准能力，掌握C指数计算、Bootstrap内部验证及外部验证技术。本方案提供可复现的代码示例，帮助研究人员提升生存分析模型的准确性。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本研究由医学科学基金资助 &amp;浙江省科技计划项目（项目编号：2013KYA212）、浙江省自然科学基金一般项目（项目编号：Y19H160126）以及金华市科学技术局重点项目 &amp;#38; 科技局（项目编号：2016-3-005、2018-3-001d 和 2019-3-013）。

监测、流行病学和最终结果（SEER）数据库是一个开放获取的癌症数据库，仅包含去标识化的患者数据（SEER ID：12296-Nov2018）。因此，本研究免除了浙江大学医学院附属金华医院伦理审查委员会的审批要求。
1. 数据准备和 R 软件包准备
<ol>
	<li>准备并导入数据。 
	&#62; 数据集 &#60;- read.csv("…/Breast cancer Data.xlsx&#34;) #导入数据。 
	注意：数据已上传至补充文件1。</li>
	<li>安装并加载R软件包。 
	&#62; 软件包 &#60;- c(&#34;rms&#34;,&#34;cmprsk&#34;,&#34;mstate&#34;,&#34;survival&#34;,&#34;riskRegression&#34;,&#34; 
	prodlim&#34;) 
	&#62; req.pcg &#60;- 函数(pcg){ 
	新 &#60;- pcg[!(pcg %in% installed.packages()[, &#34;Package&#34;])] 
	如果 (length(新)) install.packages(新, dependencies = T) 
	sapply(pcg, require, ch = T) 
	&#160;} 
	&#62; req.pcg(软件包) 
	&#8203;注意：以下操作基于R软件（版本3.6.2）进行，使用rms、cmprsk、mstate、survival、riskRegression和prodlim软件包（http://www.r-projectrg/）。</li>
</ol>
2. 采用两种不同方法建立竞争风险列线图
<ol>
	<li>采用直接法建立竞争风险列线图。 
	&#62; mod_cph &#60;- cph(Surv(Survivalmonths, status) ~ factor1+ factor2+…, 
	&#160;x=T, y=T, surv=T, data=Dataset) 
	&#62; nom &#60;- nomogram(mod_cph, fun=list(function(x) 1-surv_cph(36, x)…), 
	&#160;funlabel=c(&#34;3-year event1 Prob."…), lp=F) #以第36个月为例。 
	&#62; mod_crr &#60;- crr(Survivalmonths, fstatus, failcode=1, cov1=cov) 
	&#62; score &#60;- log(log((1-real.3y),(1-cif.min36)))/(maxbeta/100) 
	&#62; plot(nom)</li>
	<li>采用加权法建立竞争风险列线图。 
	&#62; df.w &#60;- crprep(&#34;Survivalmonths&#34;,&#34; fstatus&#34;, 
	&#160;data=Dataset, trans=c(1,2), cens=0, 
	&#160;keep=c(&#34;factor1&#34;,&#34; factor2"…)) 
	&#62; mod.w &#60;- cph(Surv(Tstart, Tstop, status==1) ~ factor1+factor2+…, 
	&#160;data=df.w, weight=weight.cens, subset=failcode==1, surv=T) 
	​> nom.w &#60;- nomogram(mod.w…)</li>
</ol>
3. 竞争风险列线图的区分能力
<ol>
	<li>区分度C指数
	<ol>
		<li>拟合矩阵 cov 进入竞争风险模型 mod_crr 并获得一个预测矩阵 suv. 
		&#62; suv &#60;- predict.crr(mod_crr, cov)</li>
		<li>获取特定月份的累积发病率 suv 并使用该函数计算C指数 rcorr.cens. 
		&#62; cif36 &#60;- suv[which(suv[,1]==36),][-1] 
		&#62; rcorr &#60;- rcorr.cens(1-cif36, Surv(Dataset$Survivalmonths, Dataset$tumordeath)) 
		&#62; cindex &#60;- rcorr[1]</li>
	</ol>
	</li>
	<li>区分的AUC
	<ol>
		<li>使用该函数评估竞争风险模型的预测性能 评分 （riskRegression 包） 
		&#62; fgr.w &#60;- FGR(Hist(生存月数, fstatus) ~ 因子1 + 因子2 + …, data=数据集, cause=1) 
		&#62; 评分 &#60;- 评分（列表（&#34;Fine-Gray&#34; = fgr.w),</li>
		<li>从……中提取AUC &#34;评分&#34;. 
		> score$AUC</li>
	</ol>
	</li>
</ol>
4. 竞争风险模型的校准能力
<ol>
	<li>具有95%置信区间的竞争风险模型校准曲线
	<ol>
		<li>获取一个数据框，其中包含每个个体在特定失效时间的累积发病率。 
		&#62; cif36 &#60;- data.frame(cif36) 以第36个月为例。 
		&#62; cif36.36_o 的列名&#60;-c(&#34;36m&#34;)</li>
		<li>根据估计的累积发病率将队列分为五个亚组，并计算每个亚组的平均预测累积发病率。 
		&#62; group36 &#60;- cut(cif36$`36m`, 
		quantile(cif36$`36m`, seq(0, 1, 0.2)) 
		include.lowest = TRUE, labels = 1:5) 
		&#62; 平均值36 &#60;- as.vector(by(cif36$`36m`, group36, mean))</li>
		<li>使用该函数计算观察到的累积发病率，即实际累积发病率 姜黄素，然后在特定失效时间点获得带有95%置信区间的观察累积发病率。 
		&#62; cum36 &#60;- cuminc(Dataset$Survivalmonths, Dataset$fstatus, group36) 
		&#62; obs36 &#60;- 时间点(cum36, Dataset$Survivalmonths)$est[c(1:5),36] 
		&#62; obs36var &#60;- timepoints(cum36, Dataset$Survivalmonths)$var[c(1:5), 36] 
		&#62; df &#60;- 数据框（mean36，obs36，obs36var）</li>
		<li>使用该函数以预测累积发病率作为 x 轴，观察累积发病率作为 y 轴绘制校准曲线 ggplot. 
		&#62; ggplot(df) + geom_point(aes(x=mean36, y=obs36), col=&#34;红色&#34;)+ 
		&#160;geom_point(aes(x=mean36, y=obs36), col=&#34;红色&#34;,pch=4)+ 
		&#160;geom_line(col=&#34;红色&#34;，aes(x=均值36，y=观测值36))+ 
		&#160;geom_errorbar(col=&#34;红色&#34;，aes(x=均值36，y=观测值36+1.96 
		*sqrt(obs36var)) 
		&#160;ymin = obs36 - 1.96 × √(obs36var)，ymax = obs36 + 1.96 
		*sqrt(obs36var) 
		geom_abline(lty=3,lwd=2,col=c(rgb(0,118,192, 
		maxColorValue=255)))</li>
	</ol>
	</li>
	<li>竞争风险模型风险评分的校准曲线
	<ol>
		<li>评估所有变量的每个水平，并获得总RS值。 
		&#62; Dataset$factor1[Dataset$factor1==1] &#60;- factor1.scale[&#34;Factor1_水平1&#34;] 
		&#62; … #例如，Dataset$histology[Dataset$histology==1]&#60;-组织学.比例尺[&#34;组织学1&#34;] 
		&#62; Dataset$rs &#60;- Dataset$factor1 + Dataset$factor2 + Dataset$factor3 + … 
		注意：通过将每位患者的各变量得分相加，计算其总RS得分。</li>
		<li>统计不同总风险评分的频次，并计算其观察到的累积发病率。 
		&#62; rs.freq &#60;- as.data.frame(table(数据集$rs)) 
		&#62; obs.36 &#60;- 向量（模式=&#34;数值的&#34;，长度 = nrow(rs.freq) 
		&#62; for (i in 1:nrow(rs.freq)) { 
		&#160;数据集 &#60;- subset(数据集, 数据集$rs == rs.freq[i,1]) 
		&#160;cif.dataset &#60;- cuminc(数据集$生存月数, 数据集$死亡3) 
		&#160;cif36.dataset &#60;- 时间点(cif.dataset,36) 
		&#160;obs.36[i] &#60;- cif36.dataset$est[1]}</li>
		<li>设置x轴的范围，并计算总风险评分的预测累积发病率。 
		&#62; RS &#60;- range(nom$total.points) 
		&#62; x.36 &#60;- seq(min(RS), max(RS), 0.01) 
		&#62; pre.36 &#60;- 1-(1-cif.min36)^exp(x.36*maxbeta/100)</li>
		<li>以风险评分绘制校准曲线。 
		&#62; plot(x.36, pre.36, type='l'…) 
		&#62; par(new=TRUE) 
		> plot(as.vector(rs.freq[,1]), obs.36…)</li>
	</ol>
	</li>
</ol>
5. 竞争风险模型的决策曲线分析
<ol>
	<li>载入 stdca 函数以执行决策曲线分析。 
	&#62; source(&#34;stdca.R&#34;)</li>
	<li>从列线图中提取多项式方程，用于计算生存概率。 
	&#62; nomogramEx(nomo = nom) 
	&#62; Dataset$predictors &#60;- A * (Dataset$rs ^3) + B * (Dataset$rs ^2) + C * Dataset$rs + D #predictors 为通过已建立的列线图计算得出的癌症特异性死亡预测概率</li>
	<li>执行决策曲线分析。 
	&#62; stdca(data = Dataset, outcome = &#34;status&#34;, ttoutcome = &#34;Survivalmonths&#34;, timepoint = 36, 
	&#160;predictors = &#34;predictors&#34;, cmprsk = TRUE, smooth = FALSE, probability = FALSE) 
	注：在存在竞争风险的情况下评估结局时，应将 cmprsk 设为 TRUE。</li>
</ol>
6. 使用自助法进行内部验证
<ol>
	<li>使用自助法获取平均预测累积发病率。
	<ol>
		<li>重采样原始数据集（数据集）与 替换 用于生成自助数据集（数据集_in）。建立竞争风险模型（mod.in_crr) 与自举数据集一起使用该函数 predict.crr 预测 mod.in_crr 并循环 b 次以生成 suvall.in. 
		B=b 
		suvall.in &#60;- 列表() 
		for(j in 1:B){ 
		&#160;数据集_in &#60;- 数据集[样本(c(1:nrow(数据集)), nrow(数据集), 
		replace = TRUE),] 
		连接（数据集_输入） 
		&#160;覆盖范围 &#60;- model.matrix(~因子1 + 因子2 + …)[,-1] 
		&#160;_crr 中的修饰 &#60;- crr(生存月数, fstatus, failcode=1, cov1=cov.in) 
		分离（Dataset. inner） 
		suv. in &#60;- predict.crr(mod. in _crr, cov) 
		suvall.in[[j]] &#60;- suv.in}</li>
		<li>获取某个月份的平均预测累积发病率。 
		cif36all. 内层 &#60;- 向量（模式=&#34;数值的&#34;，长度=nrow(数据集) 
		对于（k 从 1 到 B）{ 
		&#160;cif36all. 内部&#60;- cif36all. inner + suvall. inner[[k]][which(suvall. inner[[k]][,1] == 36),][-1] 
		} 
		cif36.in &#60;- cif36all.in/B</li>
	</ol>
	</li>
	<li>使用函数通过内部交叉验证计算C指数 rcorr.cens. 
	rcorr. 内部 &#60;- rcorr.cens(1-cif36.in, Surv(Dataset$Survivalmonths, Dataset$tumordeath)) 
	cindex. 内部 &#60;- rcorr. inner[1]</li>
	<li>使用交叉内部验证进行校准。 
	&#8203;注意：具有内部验证的竞争风险模型校准曲线的代码与第4节中的代码类似，而 suv 被替换为 suv.in.</li>
</ol>
7. 竞争风险模型的外部验证
<ol>
	<li>利用外部数据获取预测的累积发病率。使用外部数据变量矩阵（cov.ex）计算预测的累积发病率。 
	suv.ex &#60;- predict.crr(mod_crr,cov.ex) 
	cif36.ex &#60;- suv.ex [which(suv.ex $time==&#34;36&#34;),][-1]</li>
	<li>使用外部验证计算C指数。 
	rcorr.ex &#60;- rcorr.cens(1-cif36.ex,Surv(Dataset.ex$Survivalmonths,Dataset.ex$tumordeath)) 
	cindex.ex &#60;- rcorr.ex[1]</li>
	<li>使用外部验证进行校准。 
	注意：具有内部验证的竞争风险模型的校准曲线代码与第4节中的代码类似，只是将suv替换为suv.ex。</li>
</ol>

<ol><li>Andersen, P. K., Gill, R. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cox%27s+regression+model+for+counting+processes%3A+A+large+sample+study%5BTitle%5D)+AND+%22The+Annals+of+Statistics%22%5BJournal%5D)">Cox's regression model for counting processes: A large sample study.</a> The Annals of Statistics. 10 (4), 1100-1120 (1982).</li>
<li>Lubsen, J., Pool, J., vander Does, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+practical+device+for+the+application+of+a+diagnostic+or+prognostic+function%5BTitle%5D)+AND+%22Methods+of+Information+in+Medicine%22%5BJournal%5D)">A practical device for the application of a diagnostic or prognostic function.</a> Methods of Information in Medicine. 17 (2), 127-129 (1978).</li>
<li>Harrell, F. E., Lee, K. L., Mark, D. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multivariable+prognostic+models%3A+Issues+in+developing+models%2C+evaluating+assumptions+and+adequacy%2C+and+measuring+and+reducing+errors%5BTitle%5D)+AND+%22Statistics+In+Medicine%22%5BJournal%5D)">Multivariable prognostic models: Issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors.</a> Statistics In Medicine. 15 (4), 361-387 (1996).</li>
<li>Hung, H., Chiang, C. -T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Estimation+methods+for+time-dependent+AUC+models+with+survival+data%5BTitle%5D)+AND+%22The+Canadian+Journal+of+Statistics+%2F+La+Revue+Canadienne+de+Statistique%22%5BJournal%5D)">Estimation methods for time-dependent AUC models with survival data.</a> The Canadian Journal of Statistics / La Revue Canadienne de Statistique. 38 (1), 8-26 (2010).</li>
<li>Moons, K. G. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Risk+prediction+models%3A+I.+Development%2C+internal+validation%2C+and+assessing+the+incremental+value+of+a+new+%28bio%29marker%5BTitle%5D)+AND+%22Heart%22%5BJournal%5D)">Risk prediction models: I. Development, internal validation, and assessing the incremental value of a new (bio)marker.</a> Heart. 98 (9), 683-690 (2012).</li>
<li>Fu, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Real-world+impact+of+non-breast+cancer-specific+death+on+overall+survival+in+resectable+breast+cancer%5BTitle%5D)+AND+%22Cancer%22%5BJournal%5D)">Real-world impact of non-breast cancer-specific death on overall survival in resectable breast cancer.</a> Cancer. 123 (13), 2432-2443 (2017).</li>
<li>Fine, J. P., Gray, R. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+proportional+hazards+model+for+the+subdistribution+of+a+competing+risk%5BTitle%5D)+AND+%22Journal+of+the+American+Statistical+Association%22%5BJournal%5D)">A proportional hazards model for the subdistribution of a competing risk.</a> Journal of the American Statistical Association. 94 (446), 496-509 (1999).</li>
<li>Wu, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Establishing+a+competing+risk+regression+nomogram+model+for+survival+data%5BTitle%5D)+AND+%22Journal+of+Visualized+Experiments%22%5BJournal%5D)">Establishing a competing risk regression nomogram model for survival data.</a> Journal of Visualized Experiments. (164), e60684(2020).</li>
<li>Zhang, Z., Geskus, R. B., Kattan, M. W., Zhang, H., Liu, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nomogram+for+survival+analysis+in+the+presence+of+competing+risks%5BTitle%5D)+AND+%22Annals+of+Translational+Medicine%22%5BJournal%5D)">Nomogram for survival analysis in the presence of competing risks.</a> Annals of Translational Medicine. 5 (20), 403(2017).</li>
<li>Zhang, Z. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Overview+of+model+validation+for+survival+regression+model+with+competing+risks+using+melanoma+study+data%5BTitle%5D)+AND+%22Annals+Of+Translational+Medicine%22%5BJournal%5D)">Overview of model validation for survival regression model with competing risks using melanoma study data.</a> Annals Of Translational Medicine. 6 (16), 325(2018).</li>
<li>Newson, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Confidence+intervals+for+rank+statistics%3A+Somers%27+D+and+extensions%5BTitle%5D)+AND+%22Stata+Journal%22%5BJournal%5D)">Confidence intervals for rank statistics: Somers' D and extensions.</a> Stata Journal. 6 (3), 309-334 (2006).</li>
<li>Davison, A. C., Hinkley, D. V., Schechtman, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Efficient+bootstrap+simulation%5BTitle%5D)+AND+%22Biometrika%22%5BJournal%5D)">Efficient bootstrap simulation.</a> Biometrika. 73 (3), 555-566 (1986).</li>
<li>Roecker, E. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Prediction+error+and+its+estimation+for+subset-selected+models%5BTitle%5D)+AND+%22Technometrics%22%5BJournal%5D)">Prediction error and its estimation for subset-selected models.</a> Technometrics. 33 (4), 459-468 (1991).</li>
<li>Steyerberg, E. W., Harrell, F. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Prediction+models+need+appropriate+internal%2C+internal-external%2C+and+external+validation%5BTitle%5D)+AND+%22Journal+of+Clinical+Epidemiology%22%5BJournal%5D)">Prediction models need appropriate internal, internal-external, and external validation.</a> Journal of Clinical Epidemiology. 69, 245-247 (2016).</li>
<li>Zhang, Z., Chen, L., Xu, P., Hong, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Predictive+analytics+with+ensemble+modeling+in+laparoscopic+surgery%3A+A+technical+note%5BTitle%5D)+AND+%22Laparoscopic%2C+Endoscopic+and+Robotic+Surgery%22%5BJournal%5D)">Predictive analytics with ensemble modeling in laparoscopic surgery: A technical note.</a> Laparoscopic, Endoscopic and Robotic Surgery. 5 (1), 25-34 (2022).</li>
</ol>

作者声明不存在任何竞争利益。

本研究比较了通过两种不同方法建立的竞争风险列线图，并对所建立的列线图进行了评估与验证。具体而言，本研究提供了基于直接法建立列线图的逐步教程，以及计算C指数和绘制校准曲线的方法。
R 软件中的 rms 包被广泛用于构建和评估 Cox 比例风险模型，但不适用于竞争风险模型。对于具有多个结局的模型，Zhang 等人报道了使用 R 软件中的 riskRegression 包对竞争风险模型进行验证的方法，该方法通过计算 AUC 和 Brier 评分来评估模型的区分能力，并绘制校准曲线以验证其校准能力9,10。
然而，riskRegression 软件包存在一些不足之处，例如无法计算C指数。C指数是指在所有可评估的有序对中，一致有序对所占的比例（C指数 = 一致可评估有序对 / 所有可评估有序对）3,...

近年来，随着精准医学的发展，不断有新的预后因素被发现，结合分子与临床病理因素的预后模型在临床实践中正受到越来越多的关注。然而，诸如Cox比例风险模型等非图形化模型，其结果以系数值形式呈现，临床医生难以直观理解1。相比之下，列线图（nomogram）是一种回归模型（包括Cox回归模型、竞争风险模型等）的可视化工具，是一种用于近似图解计算数学函数的二维图示2。它能够量化临床模型中各变量的不同水平，并计算风险评分（RS），从而预测患者的预后。
模型评估是模型构建中的关键环节，通常有两个被广泛接受的评估特征：区分度和校准度。在临床模型中，区分度指的是模型将发生事件的个体与未发生事件的个体区分开的能力，例如死亡患者与存活患者的区分，通常使用一致性指数（C指数）或受试者工作特征曲线下面积（AUC）来表征3,4。校准度是指将模型预测的概率与实际观察到的概率进行比较的过程，校准曲线已被广泛用于展示这一特性。此外，模型验证（包括内部验证和外部验证）是模型构建中的重要步骤，只有经过验证的模型才能进一步推广应用5。
Cox 比例风险模型是一种在医学研究中用于分析预后因素与生存状态之间关联的回归模型。然而，Cox 比例风险模型仅考虑两种结局状态 [Y (0, 1)]，而研究对象常常面临多于两种的状态，此时会出现竞争风险 [Y (0, 1, 2)]1。总生存期（OS）定义为从起始日期（例如治疗开始）到因任何原因死亡的时间，是生存分析中最重要的终点指标。然而，OS 无法区分癌症特异性死亡与非癌症特异性死亡（如心血管事件及其他无关原因），从而忽略了竞争风险的影响6。在这些情况下，应优先采用竞争风险模型来预测考虑了竞争风险的生存状态7。Cox 比例风险模型的构建与验证方法已较为成熟，而关于竞争风险模型的验证报道则较少。
在我们之前的研究中，建立了一个特定的竞争风险列线图，该列线图结合了列线图与竞争风险模型，并基于竞争风险模型构建了风险评分估计方法8。本研究旨在介绍已建立的竞争风险列线图的不同评估与验证方法，该列线图可作为临床医生在考虑竞争风险的前提下预测预后的实用工具。

<table><tbody><tr><td>R 软件</td><td>无</td><td>不适用</td><td>版本 3.6.2 或更高&amp;nbsp;</td></tr><tr><td>计算机系统</td><td>Microsoft&amp;nbsp;</td><td>Windows 10</td><td>&amp;nbsp;Windows 10 或更高版本</td></tr></tbody></table>

基于r语言的竞争风险模型景观验证

An R-Based Landscape Validation of a Competing Risk Model

Cox比例风险模型在临床环境中广泛应用于生存分析，但其无法处理多种生存结局。与传统的Cox比例风险模型不同，竞争风险模型考虑了竞争事件的存在，并结合列线图（nomogram）——一种图形化计算工具，为临床医生提供精确的预后预测。本研究报道了一种建立竞争风险列线图的方法，包括评估其区分能力（即一致性指数和曲线下面积）、校准能力（即校准曲线）以及净收益（即决策曲线分析）。此外，还通过重抽样原始数据集进行内部验证，并利用外部数据集对已建立的竞争风险列线图进行外部验证，以展示其外推能力。该竞争风险列线图可作为临床医生在考虑竞争风险因素下预测预后的有力工具。

本研究从SEER数据库中获取了乳腺癌患者的数据，并将其作为示例数据。SEER数据库提供了代表美国约34.6%人口的癌症数据，且已获得访问该数据库的许可（参考编号：12296-Nov2018）。
建立了两个列线图（图1），二者均包含组织学类型、分化程度、T分期和N分期，分别采用直接法和加权法构建。各变量水平的评分及其对应的总评分概率几乎相同，但也可观察到一些细微差异。Zhang等人提出了一种“加权”方法用于构建竞争风险列线图，该方法首先将原始数据转换为加权数据（使用函数crprep），然后利用加权数据构建Cox回归模型（使用函数coxph），最后基于Cox回归模型建立竞争风险列线图（使用函数nomogram）9。相比之下，本研究中的“直接”方法与“加权”方法完全不同。简而言之，本研究使用竞争风险模型（使用函数crr）生成的参数替代了Cox回归模型（使用函数coxph）中的参数，最终利用这些...

本方案介绍了用于评估竞争风险模型区分能力和校准能力的 R 语言代码，以及对该模型进行内部和外部验证的代码。

基于R语言的竞争风险模型景观验证

Cox 比例风险模型在临床环境中广泛应用于生存分析，但它无法处理多重生存……

本方案介绍了多种评估和验证竞争风险列线图的方法，该列线图在生存分析中考虑了竞争事件的存在。本方案可作为 R 语言中风险回归（risk regression）软件包的补充，例如通过自助法（bootstrap method）进行 C 指数计算及内部验证。首先进行 C 指数区分度分析时，需将矩阵 cob 拟合到竞争风险模型 mod CRR 中，并通过执行相应命令获得预测矩阵 SUV。从 SUV 获取特定月份的累积发病率，并使用 R-core 提供的函数计算 C 指数。对于 AUC 区分度，使用 risk regression 包中的 score 函数评估竞争风险模型的预测性能，然后通过执行命令从 score 中提取 AUC。为了获得竞争风险模型在95%置信区间内的校准曲线，首先获取一个数据框，其中包含特定失效时间下每个个体的累积发生率。然后根据估计的累积发生率将队列分为五个亚组，并计算每个亚组的平均预测累积发生率。使用 cuminc 函数计算观察到的累积发生率，即实际累积发生率。然后通过执行命令，获得在特定失效时间下带有95%置信区间的观察累积发病率。以预测累积发病率为X轴，观察累积发病率为Y轴，绘制校准曲线。使用ggplot函数。对于竞争风险模型的风险评分校准曲线，评估所有变量的每个水平，并通过执行命令获得总风险评分（RS）。统计不同总风险评分的频数，计算其观察到的累积发生率。设定X轴的范围，并计算总风险评分的预测累积发生率。然后通过执行命令，以风险评分绘制校准曲线。为了使用自助法（bootstrap method）获得平均预测累积发生率，需对原始数据集进行有放回重采样，生成自助数据集。数据集输入。然后建立一个竞争风险模型：使用Bootstrap数据集拟合mod NCRR，并在循环中B次调用函数predict CRR来对mod NCRR进行预测，以生成SUV all in。接着，计算特定月份的平均预测累积发生率。使用R-core sends函数通过区间交叉验证计算C指数。进行校准，使用外部验证，通过执行命令，利用外部数据及外部数据变量矩阵获得预测累积发生率：Code-X。然后通过执行命令，使用外部验证计算C指数。采用直接法和加权法得到两个列线图，结果显示各变量水平对应的评分以及总评分对应的发生概率几乎相同。尽管观察到一些轻微差异，竞争风险模型的校准曲线接近于一致性直线，且各组中观察频率的95%置信区间均落在一致性直线上，表明该模型具有良好的校准能力。使用内部和外部验证绘制的校准曲线显示，所构建的模型在内部验证中具有良好的校准能力，但在外部验证中校准能力较差。获得了竞争风险列线图的决策曲线分析结果，展示了随着阈值概率的增加，净收益的变化情况。在进行竞争风险列线图的内部验证时，通过重采样原始数据集生成bootstrap数据集至关重要。除了bootstrap方法外，还可采用随机化展示和K4方法来生成用于内部验证的数据集。利用我们基于R语言的竞争风险模型全景验证方法，临床医生在实际应用中能够更便捷地考虑竞争风险，从而进行预后分析。

an-r-based-landscape-validation-of-a-competing-risk-model

Research

科研

JoVE Journal

Cancer Research

癌症研究

竞争风险模型的基于 R 的景观验证

面对未来不确定的环境条件，迫切需要能够管理水生系统的工具和方法。本文提供了开展针对性流域评估的方法，使资源管理者能够构建基于景观的累积效应模型，并将其应用于情景分析管理框架中。...

基于定量情景分析框架的流域规划

本文介绍了一种基于Cox比例风险回归模型和竞争风险回归模型构建列线图的方案。当生存分析中存在竞争事件时，竞争风险方法是一种更为合理的分析方法。...

建立用于生存数据的竞争风险回归列线图模型

基于 CogStack 利用电子健康记录实现精神病风险实时检测与预警系统

该研究发现，在多原发结直肠癌（MPCC）患者术后，男性性别、肿瘤分级差以及较晚的肿瘤-淋巴结-转移分期与较差的癌症特异性生存率（CSS）相关。我们建立了一个列线图，用于预测MPCC患者的CSS，以辅助临床治疗决策。...

用于预测术后多原发结直肠癌患者癌症特异性生存的竞争风险列线图

本研究利用机器学习模型和竞争风险分析，评估了结直肠印戒细胞癌患者的预后系统。研究发现，与pN分期相比，阳性淋巴结的对数比（LODDS）是更优的预后预测指标，具有良好的预测性能，可通过可靠的生存预测工具辅助临床决策。...

基于机器学习模型的三种淋巴结分期系统在结直肠印戒细胞癌中的预测性能比较

我们描述了通过将肿瘤细胞或类器官注射到小鼠盲肠中建立原位结直肠肿瘤模型，并从此模型中后续分离循环肿瘤细胞（CTCs）的方法。...

结直肠癌原位小鼠模型中循环肿瘤细胞的分离

一个 在体 内 小鼠坐骨神经周围浸润模型

小鼠原位胰腺导管腺癌同种移植瘤免疫表型分析的实验步骤旨在对肿瘤免疫微环境进行特征分析。通过手术将肿瘤原位植入小鼠体内。收获体积为200–600...

通过流式细胞术对小鼠原位同系移植（同源性）KPC 胰腺导管腺癌进行免疫表型分析

本文提供了一种基于细胞三维形态来鉴别健康与病变细胞的工作流程。我们描述了如何利用基于三维表面生成的二维投影轮廓，训练自组织映射网络，从而实现对所研究细胞群体的客观聚类。...

基于形态学的健康与病变细胞区分方法：利用傅里叶变换与自组织映射

本文介绍的诊断用FL-DNA试剂盒是一种省时且操作简便的方法，可用于确定结直肠癌病变存在的可靠概率。...

通过粪便DNA完整性检测评估结直肠癌风险与患病率

学习使用R语言评估竞争风险模型的判别与校准能力，掌握C指数计算、Bootstrap内部验证及外部验证技术。本方案提供可复现的代码示例，帮助研究人员提升生存分析模型的准确性。

基于R语言的竞争风险模型景观验证

探索更多视频