基于Nomogram的公式提取

# 從Nomogram提取個體的得分公式匣屡，開發(fā)一個打分系統(tǒng)
# 參考資料1：https://mp.weixin.qq.com/s?src=11&timestamp=1602050809&ver=2629&signature=EAoh1zmeYnPQMazJcBumM7695mseQ3L45d44zN8P1*9bfPGpZFBX1SNarWynYwc8WlHdiaTiHKfTmQhoqgzcyr0wB1DIQ3vLb-8COSbhSO*BAIktlZP1GHnsOLt3AzZv&new=1
# 參考資料2：https://mp.weixin.qq.com/s?src=11&timestamp=1602050809&ver=2629&signature=2pteahb*mqUZo4jeAJhaLi1l98hNj*fv-OrO7eYPudwq7eelf0C-XBi0aWzCbBHw2vQBSyK5wxwHDJic5mjxWuoyHgHJuOYeSgllLMbk--9nXQkWuuzkGYg0DRbepi08&new=1
## 
library(rms)
n <-1000
age <- rnorm(n,50,10)
sex <- factor(sample(c('female','male'),n,TRUE))
sex <- as.numeric(sex)
ddist <- datadist(age,sex)
options(datadist='ddist')
cens <- 15*runif(n)
time <- -log(runif(n))/0.02*exp(.04*(age-50)+.8*(sex=='Female'))
death <- ifelse(time <= cens,1,0)
time <- pmin(time,cens)
units(time)="month"
f <- cph(formula(Surv(time,death)~sex+age),x=TRUE,y=TRUE,surv=TRUE,time.inc=3)
surv <- Survival(f)
nomo <- nomogram(f, fun=list(function(x) surv(3,x),function(x) surv(6,x)),
                 lp=TRUE,funlabel=c("3-Month Survival Prob","6-Month Survival Prob"))

plot(nomo)

## NomogramEx
#安裝nomogramEx包
#install.packages("nomogramEx")
library(nomogramEx)
nomogramEx(nomo=nomo)

## nomogramFormula
#install.packages("nomogramFormula")
library(nomogramFormula)

## 單個患者列線圖得分
#options(option)
results <- formula_lp(nomogram = nomo)
## 公式提取
results$formula
## 計(jì)算得分
points <- points_cal(formula = results$formula, lp = f$linear.predictors)
head(points)

## 計(jì)算得分
#options(option)
results <- formula_rd(nomogram = nomo)
mydata$points <- points_cal(formula = results$formula,rd=mydata)
head(mydata$points)
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 人面猴
    序言：七十年代末，一起剝皮案震驚了整個濱河市鹅士，隨后出現(xiàn)的幾起案子，更是在濱河造成了極大的恐慌也拜，老刑警劉巖趾痘，帶你破解...
    沈念sama閱讀 216,402評論 6贊 499
  • 死咒
    序言：濱河連續(xù)發(fā)生了三起死亡事件，死亡現(xiàn)場離奇詭異卵贱，居然都是意外死亡侣集，警方通過查閱死者的電腦和手機(jī)，發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,377評論 3贊 392
  • 救了他兩次的神仙讓他今天三更去死
    文/潘曉璐 我一進(jìn)店門编振，熙熙樓的掌柜王于貴愁眉苦臉地迎上來臭埋，“玉大人，你說我怎么就攤上這事斋泄。” “怎么了魁莉？”我有些...
    開封第一講書人閱讀 162,483評論 0贊 353
  • 道士緝兇錄：失蹤的賣姜人
    文/不壞的土叔 我叫張陵募胃，是天一觀的道長。 經(jīng)常有香客問我痹束，道長，這世上最難降的妖魔是什么祷嘶？ 我笑而不...
    開封第一講書人閱讀 58,165評論 1贊 292
  • ?港島之戀（遺憾婚禮）
    正文 為了忘掉前任夺溢，我火速辦了婚禮，結(jié)果婚禮上烛谊，老公的妹妹穿的比我還像新娘。我一直安慰自己丹禀，他們只是感情好，可當(dāng)我...
    茶點(diǎn)故事閱讀 67,176評論 6贊 388
  • 惡毒庶女頂嫁案：這布局不是一般人想出來的
    文/花漫 我一把揭開白布双泪。 她就那樣靜靜地躺著，像睡著了一般焙矛。 火紅的嫁衣襯著肌膚如雪葫盼。 梳的紋絲不亂的頭發(fā)上薄扁，一...
    開封第一講書人閱讀 51,146評論 1贊 297
  • 城市分裂傳說
    那天，我揣著相機(jī)與錄音邓梅，去河邊找鬼。 笑死日缨，一個胖子當(dāng)著我的面吹牛钱反，可吹牛的內(nèi)容都是我干的匣距。 我是一名探鬼主播面哥，決...
    沈念sama閱讀 40,032評論 3贊 417
  • 雙鴛鴦連環(huán)套：你想象不到人心有多黑
    文/蒼蘭香墨 我猛地睜開眼毅待，長吁一口氣：“原來是場噩夢啊……” “哼！你這毒婦竟也來了尸红？” 一聲冷哼從身側(cè)響起吱涉，我...
    開封第一講書人閱讀 38,896評論 0贊 274
  • 萬榮殺人案實(shí)錄
    序言：老撾萬榮一對情侶失蹤外里，失蹤者是張志新（化名）和其女友劉穎，沒想到半個月后盅蝗，有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體鳖链，經(jīng)...
    沈念sama閱讀 45,311評論 1贊 310
  • ?護(hù)林員之死
    正文 獨(dú)居荒郊野嶺守林人離奇死亡墩莫，尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,536評論 2贊 332
  • ?白月光啟示錄
    正文 我和宋清朗相戀三年逞敷，在試婚紗的時候發(fā)現(xiàn)自己被綠了题山。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片兰粉。...
    茶點(diǎn)故事閱讀 39,696評論 1贊 348
  • 活死人
    序言：一個原本活蹦亂跳的男人離奇死亡顶瞳，死狀恐怖，靈堂內(nèi)的尸體忽然破棺而出愕秫，到底是詐尸還是另有隱情，我是刑警寧澤戴甩，帶...
    沈念sama閱讀 35,413評論 5贊 343
  • ?日本核電站爆炸內(nèi)幕
    正文 年R本政府宣布，位于F島的核電站甜孤，受9級特大地震影響，放射性物質(zhì)發(fā)生泄漏缴川。R本人自食惡果不足惜茉稠，卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,008評論 3贊 325
  • 男人毒藥：我在死后第九天來索命
    文/蒙蒙 一把夸、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧恋日，春花似錦膀篮、人聲如沸岂膳。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,659評論 0贊 22
  • 一樁弒父案谈截，背后竟有這般陰謀
    文/蒼蘭香墨 我抬頭看了看天上的太陽筷屡。三九已至，卻和暖如春速蕊，著一層夾襖步出監(jiān)牢的瞬間，已是汗流浹背娘赴。 一陣腳步聲響...
    開封第一講書人閱讀 32,815評論 1贊 269
  • 情欲美人皮
    我被黑心中介騙來泰國打工， 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留诽表，地道東北人隅肥。 一個月前我還...
    沈念sama閱讀 47,698評論 2贊 368
  • 代替公主和親
    正文 我出身青樓袄简，卻偏偏與公主長得像腥放，于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子秃症，可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,592評論 2贊 353