Displaying 1 result from an estimated 1 matches for "261ff".
Did you mean:
21ff
2004 Dec 16
0
fitting problems in coxph.fit
...uot;survival")
looks much better
R> coxfit$coef
[1] 0.25123203 0.00000000 -0.42595541 -0.04488913 -0.26061995
0.44426458
[7] -0.38954286 -0.43081374 0.79573107 0.48234405 0.94636357
-0.25193465
[13] 1.15619712 0.32651765 -1.06731019 -0.24249939
I nailed the problem down to lines 261ff of `coxfit2.c' where
zbeta = offset[person];
for (i=0; i<nvar; i++)
zbeta += newbeta[i]*covar[i][person];
risk = exp(zbeta ) * weights[person];
and `zbeta' may become very large and thus `exp' returns `nan'. Of course
one coul...