目的：

使用R語言結合某些包實現nomogram（列線圖）繪制。

緣起：

由于偶然的機會，學到了這種圖的畫法，及其解讀，在此作為筆記記錄。感謝前輩的無所分享，本文學習自：
https://www.cnblogs.com/biostatisc/p/7903160.html
nomogram的原理可在上述鏈接中查看。（如侵犯版權，請留言，鄙人會刪除）

應用范圍：

logistic回歸，cox回歸的結果展示。

自己的理解：

這種圖特別適使用于logistic回歸的結果展示，是“評分系統”的一個替代，可簡便應用于臨床等實踐。

一個可重復的例子：

（引用于https://www.cnblogs.com/biostatisc/p/7903160.html，如侵犯版權，請留言，鄙人會刪除）

    require(rms)    ##調用rms包

    # 建立數據集
    y = c(0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1,
    1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1,
    1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0,
    0, 0, 1, 0, 1, 0, 1, 0, 1)
    age = c(28, 42, 46, 45, 34, 44, 48, 45, 38, 45, 49, 45, 41, 46, 49, 46, 44, 48,
    52, 48, 45, 50, 53, 57, 46, 52, 54, 57, 47, 52, 55, 59, 50, 54, 57, 60,
    51, 55, 46, 63, 51, 59, 48, 35, 53, 59, 57, 37, 55, 32, 60, 43, 59, 37,
    30, 47, 60, 38, 34, 48, 32, 38, 36, 49, 33, 42, 38, 58, 35, 43, 39, 59,
    39, 43, 42, 60, 40, 44)
    sex = c(0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1,
    0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1,
    0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1,
    0, 1, 1, 1, 0, 1)
    ECG = c(0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1,
    0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 2, 1, 0, 0, 2, 2, 0, 0, 2, 2,
    0, 1, 2, 2, 0, 1, 0, 2, 0, 1, 0, 2, 1, 1, 0, 2, 1, 1, 0, 2, 1, 1, 0, 2,
    1, 1, 0, 2, 1, 1)
    # 設定nomogram的參數
    ddist <- datadist(age, sex, ECG)
    options(datadist='ddist')
    # logistic回歸
    f <- lrm(y ~ age + sex + ECG)
    # nomogram
    nom <- nomogram(f, fun=plogis,
    fun.at=c(.001, .01, .05, seq(.1,.9, by=.1), .95, .99, .999),
    lp=F, funlabel="Risk")
    plot(nom)

nomogram for jianshu.png

注意

本文已重復上面的例子，并使用自己的數據做出了很好的結果
還有一個cox回歸的例子本人并未重復出來，也摘錄如下，依然來自：https://www.cnblogs.com/biostatisc/p/7903160.html，（如侵犯版權，請留言，鄙人會刪除）：

require(rms)
require(Hmisc)     ##需要下載安裝

require(survival)   ##R默認自帶的用于做生存分析的包
# 建立數據集（使用rms包example的代碼，未改動）
n <- 1000
set.seed(731)
age <- 50 + 12*rnorm(n)
label(age) <- "Age"
sex <- factor(sample(c('Male','Female'), n,
rep=TRUE, prob=c(.6, .4)))
cens <- 15*runif(n)
h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
dt <- -log(runif(n))/h
label(dt) <- 'Follow-up Time'
e <- ifelse(dt <= cens,1,0)
dt <- pmin(dt, cens)
units(dt) <- "Year"
# 設定nomogram的參數
ddist <- datadist(age, sex)
options(datadist='ddist')
# Cox回歸
S <- Surv(dt,e)
f <- cph(S ~ rcs(age,4) + sex, x=T, y=T)
med <- Quantile(f)
# nomogram
nom <- nomogram(f, fun=function(x) med(x),
fun.at=c(13,12,11,9,8,7,6,5),lp=F, funlabel="Median Survival Time")
plot(nom) ##繪制Nomgram圖