作者iterate (碎形)
看板Statistics
标题[程式] R coxph model 的 error
时间Tue May 15 14:21:36 2012
各位好
目前以 R 在执行存活分析（coxph）
程式如下：
coxpp <- NULL
for(i in 1:100000){
result <- coxph(Surv(da[,22],da[,20])~da[,14]+t(car[i,])),na.action = na.omit)
coe <- summary(result)
coee <- coe$coefficients
coxpp <- rbind(coxpp, c(i, coee[2,5]))
}
在某些 car 的变项，会出现如下的 error
错误在fitter(X, Y, strats, offset, init, control, weights = weights, :
外部函数呼叫时不能有 NA/NaN/Inf (引数 6)
此外: Warning message:
In fitter(X, Y, strats, offset, init, control, weights = weights, :
Ran out of iterations and did not converge
导致程式无法继续执行，想请教各位问题所在，及可以如何处理这样的问题！
谢谢各位先进。
da[,22]:月数（连续资料）
da[,20]:状态（类别资料 0,1）
da[,14]:类别资料（0,1）
car[i,]:类别资料（0,1,2）
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.129.5
1F：→ moonbat:你的data有NA吧，试试看在回圈前面加上na.omit(car) 05/15 16:18