R語言關(guān)于生存分析知識點(diǎn)總結(jié)
生存分析處理預(yù)測特定事件將要發(fā)生的時間。 它也被稱為故障時間分析或分析死亡時間。 例如，預(yù)測患有癌癥的人將存活的天數(shù)或預(yù)測機(jī)械系統(tǒng)將失敗的時間。
命名為survival的R語言包用于進(jìn)行生存分析。 此包包含函數(shù)Surv（），它將輸入數(shù)據(jù)作為R語言公式，并在選擇的變量中創(chuàng)建一個生存對象用于分析。 然后我們使用函數(shù)survfit()創(chuàng)建一個分析圖。
安裝軟件包
install.packages("survival")
語法
在R語言中創(chuàng)建生存分析的基本語法是
Surv(time,event) survfit(formula)
以下是所使用的參數(shù)的描述
- time是直到事件發(fā)生的跟蹤時間。
- event指示預(yù)期事件的發(fā)生的狀態(tài)。
- formula是預(yù)測變量之間的關(guān)系。
例
我們將考慮在上面安裝的生存包中存在的名為“pbc”的數(shù)據(jù)集。 它描述了關(guān)于受肝原發(fā)性膽汁性肝硬化（PBC）影響的人的生存數(shù)據(jù)點(diǎn)。 在數(shù)據(jù)集中存在的許多列中，我們主要關(guān)注字段“time”和“status”。 時間表示在接受肝移植或患者死亡的患者的登記和事件的較早之間的天數(shù)。
# Load the library.
library("survival")
# Print first few rows.
print(head(pbc))
當(dāng)我們執(zhí)行上面的代碼，它產(chǎn)生以下結(jié)果及圖表
id time status trt age sex ascites hepato spiders edema bili chol 1 1 400 2 1 58.76523 f 1 1 1 1.0 14.5 261 2 2 4500 0 1 56.44627 f 0 1 1 0.0 1.1 302 3 3 1012 2 1 70.07255 m 0 0 0 0.5 1.4 176 4 4 1925 2 1 54.74059 f 0 1 1 0.5 1.8 244 5 5 1504 1 2 38.10541 f 0 1 1 0.0 3.4 279 6 6 2503 2 2 66.25873 f 0 1 0 0.0 0.8 248 albumin copper alk.phos ast trig platelet protime stage 1 2.60 156 1718.0 137.95 172 190 12.2 4 2 4.14 54 7394.8 113.52 88 221 10.6 3 3 3.48 210 516.0 96.10 55 151 12.0 4 4 2.54 64 6121.8 60.63 92 183 10.3 4 5 3.53 143 671.0 113.15 72 136 10.9 3 6 3.98 50 944.0 93.00 63 NA 11.0 3
從上述數(shù)據(jù)，我們正在考慮分析的時間和狀態(tài)。
應(yīng)用Surv（）和survfit（）函數(shù)
現(xiàn)在我們繼續(xù)應(yīng)用Surv()函數(shù)到上面的數(shù)據(jù)集，并創(chuàng)建一個將顯示趨勢圖。
# Load the library.
library("survival")
# Create the survival object.
survfit(Surv(pbc$time,pbc$status == 2)~1)
# Give the chart file a name.
png(file = "survival.png")
# Plot the graph.
plot(survfit(Surv(pbc$time,pbc$status == 2)~1))
# Save the file.
dev.off()
當(dāng)我們執(zhí)行上面的代碼，它產(chǎn)生以下結(jié)果及圖表
Call: survfit(formula = Surv(pbc$time, pbc$status == 2) ~ 1)
n events median 0.95LCL 0.95UCL
418 161 3395 3090 3853
上圖中的趨勢有助于我們預(yù)測在特定天數(shù)結(jié)束時的生存概率。
到此這篇關(guān)于R語言關(guān)于生存分析知識點(diǎn)總結(jié)的文章就介紹到這了,更多相關(guān)R語言生存分析內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
本文主要介紹了R語言基礎(chǔ)數(shù)據(jù)類型的具體使用，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-01-01
R語言數(shù)組實(shí)例用法及知識點(diǎn)總結(jié)
在本文里，我們給大家整理的是關(guān)于R語言數(shù)組的相關(guān)知識點(diǎn)，有興趣的朋友們可以跟著學(xué)習(xí)參考下。2021-04-04
在本篇文章里小編給大家整理了一篇關(guān)于R語言繪制折線圖實(shí)例分析內(nèi)容，有需要的小伙伴可以跟著學(xué)習(xí)參考下。2021-04-04
這篇文章主要介紹了R語言數(shù)據(jù)類型深入詳解,內(nèi)容對于數(shù)據(jù)類型講解的很詳細(xì),有對于R語言感興趣的同學(xué)可以多學(xué)習(xí)下2021-03-03
R語言中矩陣matrix和數(shù)據(jù)框data.frame的使用詳解
這篇文章主要介紹了R語言中矩陣matrix和數(shù)據(jù)框data.frame的使用詳解，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2021-03-03
本文介紹了如何用R畫出帶error?bar的分組條形圖，文中為大家分享了繪制的代碼，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步2022-02-02
這篇文章主要介紹了R語言繪制坐標(biāo) 保存圖片的操作，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
在本篇文章里小編給大家整理的是一篇關(guān)于R語言繪制餅狀圖代碼實(shí)例內(nèi)容，有興趣的朋友們可以學(xué)習(xí)下。2021-04-04
R包制作后出現(xiàn)not available for錯誤問題解決解決
這篇文章主要為大家介紹了R包制作后出現(xiàn)not available for...錯誤的問題解決方式，有需要的朋友，可以借鑒參考下，希望能夠有所幫助2021-11-11
這篇文章主要介紹了R語言:排序的應(yīng)用操作，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04