[問題] PredectABEL coding

作者: ciuhkimo (羽化新生)   2017-10-08 11:06:27
ctrl + y 可以刪除一整行,請將不需要的內容刪除
使用PredictABEL package時,出現
Error in table(c1, c2, data[, cOutcome]) :
all arguments must have the same length
但是我分別用length(c1), length(c2), length(cOutcome),結果都是相同的
想請問問題可能在哪裡?
[問題類型]:
請把以下不需要的部份刪除
經驗諮詢:IDI_NRI計算作表格
[軟體熟悉度]:
新手入門
[問題敘述]:
使用PredictABEL跑IDI_NRI統計評估model
[程式範例]:
請把以下不需要的部份刪除
coding:
library(PredictABEL)
m1 <- glm(d ~ va1, family=binomial(logit), data=dt)$fitted
m2 <- glm(d ~ va1 + GM, family=binomial(logit), data=dt)$fitted
value<-as.vector(quantile(m1))
cutoff <- c(0,value[2],value[3],value[4],1)
length(m1) # 確定length是否一樣:149
length(m2) # 確定lenght:149
length(dt$va1) #length: 149
reclassification(predrisk1 = m1, predrisk2=m2,cOutcome = 3,
data = dt, cutoff = cutoff) # cOutcome =3, 第3欄為va1
Output:
_________________________________________
Reclassification table
_________________________________________
Error in table(c1, c2, data[, cOutcome]) :
all arguments must have the same length
length不是都一樣?我也有用na.omit()去把可能的NA值去掉了,
但是仍然無法
[關鍵字]:
PredictABEL, IDI, NRI, model
選擇性,也許未來有用
作者: andrew43 (討厭有好心推文後刪文者)   2017-10-11 16:27:00
沒用過,但看help中例子你是不是少了predRisk()的過程?
作者: ciuhkimo (羽化新生)   2017-10-11 21:22:00
它是利用predRisk去建立二組model,但是我用另外的方法建立以及用它的範例建立,會出現如上結果
作者: andrew43 (討厭有好心推文後刪文者)   2017-10-11 22:05:00
其實length(m1)不是你想的那樣。你應該檢查 m1$model 的 nrow,之類的。但你說你看length(m1)是149... 真的嗎?照理說 length(m1) 是 glm/lm 物件的 list 格數

Links booklink

Contact Us: admin [ a t ] ucptt.com