Code : Tout sélectionner
Ci <- cuminc(ftime=OS,fstatus=cause.death,group=arm,strata=categ,cencode=0)
En utilisant le package cmprsk, il me présente le message d'erreur suivant:
Error in qr.default(z2[[10]]) :
NA/NaN/Inf in foreign function call (arg 1)
Je me suis rendu compte que ça fonctionne quand j'enlève le "strata".
Donc du coup, j'essaie de relancer l'estimation pour chaque valeur de categ, ça marche pas toujours, principalement pour categ==1.
J'essaie de voir ce que c'est que le z2[[10]], je me rend compte que c'est une fonction étrangère programmée en Fortran:
Code : Tout sélectionner
z2 <- .Fortran("crstm", as.double(d$time), as.integer(causeind),
as.integer(d$group), as.integer(d$strata), as.integer(no),
as.double(rho), as.integer(nst), as.integer(ng),
s, v, as.double(d$time), as.integer(causeind),
as.integer(d$group), vt, s, vt, double((4 + 3 *
ng) * ng), integer(4 * ng), PACKAGE = "cmprsk")
stat[ii] <- -1
a <- qr(z2[[10]])
if (a$rank == ncol(a$qr)) {
b <- diag(dim(a$qr)[1])
stat[ii] <- z2[[9]] %*% qr.coef(a, b) %*% z2[[9]]
}
Du coup, je ne sais plus comment débugger le package pour résoudre mon problème.