前言:
建立nomogram時(shí)報(bào)錯(cuò),提示
Error in approx(fu[s], xseq[s], fat, ties = mean) :?
need at least two non-NA values to interpolate
什么是nomogram为鳄,R代碼就不寫了~
直接來到這個(gè)報(bào)錯(cuò)點(diǎn):
在刪除數(shù)據(jù)框中所有含缺失值的情況下芜飘,主要問題在于:“ function(x)surv(1*12,lp=x)”這段代碼中的1*12出現(xiàn)了問題涛贯。假如你的數(shù)據(jù)都沒有達(dá)到1*12或者都超過了1*12,那就會(huì)出錯(cuò),也就沒有預(yù)測(cè)價(jià)值矫钓。毡琉。铁瞒。
啥玩意?人話桅滋?
問題在于這里:surv1<- function(x)surv(3*365,lp=x) # 定義time.inc,n年OS慧耍,
原來的代碼是“*365”的,生存日期算的是day丐谋,但是如果你的生存數(shù)據(jù)是按月算的話芍碧,就會(huì)報(bào)錯(cuò)~
解決辦法:
一年:365(天)變成12(12個(gè)月)
surv1<- function(x)surv(3*12,lp=x) # 定義time.inc,n年OS。
后記:
能幫到你嗎号俐?幫到的話能不能給個(gè)贊嘛泌豆?
1百個(gè)贊,再詳細(xì)寫什么是nomogram吏饿,R代碼了踪危。反正網(wǎng)上都能找到蔬浙,你們看完也不點(diǎn)贊......
哭~~~~