Displaying 1 result from an estimated 1 matches for "tpseudo".
Did you mean:
pseudo
2013 Apr 09
2
R crash
...93,9.11,11.8,15.3,20.1,23.7,27.3,30.9,34.3,37.6,41.5,46.5,51.1,57,63,73.1,93.68)
a$icr <- REL
library(pseudo)
pseudo <- pseudoci(time=a$TIME,event=a$icr,tmax=cutoffs)
#rearranging data
b <- NULL
for(it in 1:length(pseudo$time)){
b <- rbind(b,cbind(a,pseudo = pseudo$pseudo[[1]][,it],
tpseudo = pseudo$time[it],id=1:nrow(a)))
}
b <- b[order(b$id),]
library(geepack)
fit3 <- geese(pseudo ~ as.factor(AGE) + as.factor(WEIGHT) +
as.factor(HEIGHT) ,
data =b, id=id, jack = TRUE, scale.fix=TRUE,
mean.link = "logit", corstr="unstructured")
#The results using the...