# ===========================================================
# Função: variância do tempo médio com correção do vício       
# Livro: Análise de Sobrevivência Aplicada - 2a edição   	   
# Autora: Vanessa Ferreira Sehaber					   
# ===========================================================

 var_tm <- function(ekm){					
   su <- summary(ekm)
   temp <- su$time
   nj <- su$n.risk      
   dj <- su$n.event
   st <- su$surv         
   k <- length(temp)-1                  
   A <- sapply(1:k, function(i){
        st[i]*(temp[i+1]-temp[i])
        })
   Aj <- sapply(1:k, function(i){
        sum(A[i:k])
        })
   den <- nj*(nj-dj); den <- den[1:k]
   num <- dj[1:k]*(Aj^2)
   m<-sum(dj)
   var <- num/den; var <- sum(var)*(m)/(m-1)
   return(var)
 }

# ========================= END =============================