This one lost X connection AFTER computation complete. I did set seed to 153 twice > BlindSimu function(X1 = rexp(100, rate = 1.6), C1 = runif(100, min = 2, max = 7), X2 = rexp(100, rate = 1.3), C2 = runif(100, min = 2, max = 9)) { Z1 <- pmin(X1, C1) D1 <- as.numeric(C1 > X1) Z2 <- pmin(X2, C2) D2 <- as.numeric(C2 > X2) temp1 <- KMNAatT(Z1, D1, 0.5) temp2 <- KMNAatT(Z2, D2, 0.5) Tvar <- (temp1$std.err)^2 + (temp2$std.err)^2 temp3 <- KMNAatT(c(Z1, Z2), c(D1, D2), 0.5) Bvar <- 4 * (temp3$std.err)^2 list(Tvariance = Tvar, Bvariance = Bvar, TnelsonVar = (temp1$ NelsonAalenVar + temp2$NelsonAalenVar), BnelsonVar = 4 * (temp3$ NelsonAalenVar), Best = temp3$surv, BGW = temp3$std.err) } > mean(result7[, 1]) [1] 0.004921417 > mean(result7[, 2]) [1] 0.00497129 > mean(result7[, 3]) [1] 0.02148737 > mean(result7[, 4]) [1] 0.02120381 > var(result7[, 5]) [1] 0.001240745 > mean((result7[, 6])^2) [1] 0.001242823 >