oberon% splus
S-PLUS : Copyright (c) 1988, 1996 MathSoft, Inc.
S : Copyright AT&T.
Version 3.4 Release 1 for HP 9000 Series, HP-UX 9.0.5 : 1996 
Working data will be in /u/disk13/stat/mai/.Data 
> source("BlindSimu.r")
> BlindSimu
function(X1 = rexp(100), C1 = rexp(100), X2 = rexp(100, rate = 1.15), C2 = rexp(
        100, rate = 1.3))
{
        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)
}
> BlindSimu()
$Tvariance:
[1] 0.007033941

$Bvariance:
[1] 0.007076371

$TnelsonVar:
[1] 0.03418437

$BnelsonVar:
[1] 0.03143522

$Best:
[1] 0.4716073

$BGW:
[1] 0.04206058

> set.seed(153)
> set.seed(153)
> result2 <- matrix(NA, nrow=5000, ncol=6)
> for(i in 1:1000) result2[i,] <- unlist(BlindSimu())
> for(i in 1001:2000) result2[i,] <- unlist(BlindSimu())
> for(i in 2001:3000) result2[i,] <- unlist(BlindSimu())
> for(i in 3001:4000) result2[i,] <- unlist(BlindSimu())
> for(i in 4001:5000) result2[i,] <- unlist(BlindSimu())  
Syntax error: unknown expression ("\033") used illegally at this point:
for(i in 4001:5000) result2[i,] <- unlist(BlindSimu())
> 
> 
> for(i in 4001:5000) result2[i,] <- unlist(BlindSimu())
> 
> mean(result2[,1])
[1] 0.006676352
> mean(result2[,2])
[1] 0.006689299
> mean(result2[,3])
[1] 0.02013357
> mean(result2[,4])
[1] 0.01972328
> var( result2[,5])
[1] 0.001620016
> mean( (result2[,6])^2 )
[1] 0.001672325
> 

> X11()
> plot(result2[,2]-result2[,1])
> Starting to make postscript file.
Finished postscript file, executing command "lpr -Pstat ps.out.0001.ps &".

> q()