Displaying 1 result from an estimated 1 matches for "ent2".
Did you mean:
ent
2008 Jul 15
0
implementation of Prentice method in cch()
...h(), the prentice
method is implemented like this:
Prentice <- function(tenter, texit, cc, id, X, ntot,robust){
eps <- 0.00000001
cens <- as.numeric(cc>0) # Censorship indicators
subcoh <- as.numeric(cc<2) # Subcohort indicators
## Calculate Prentice estimate
ent2 <- tenter
ent2[cc==2] <- texit[cc==2]-eps
fit1 <- coxph(Surv(ent2,texit,cens)~X,eps=eps,x=TRUE)
## Calculate Prentice estimate and variance
nd <- sum(cens) # Number of failures
nc <- sum(subcoh) # Number in subcohort
ncd <- sum(cc==1) #Number of failures in...