> # Nonparametric Bootstrap using bootstrap()
> 
> # function to compute the correlation and Fisher's z-transformation
> 
> cor2 <- function(X){
+ 	z1 <- X[,1]
+ 	z2 <- X[,2]
+ 	a <- cor(z1,z2)
+ 	b <- 0.50*log((1+a)/(1-a))
+ 	return(c(a,b))
+ }
> 
> r.boot <- bootstrap(X, cor2)
Forming replications  1  to  100 
Forming replications  101  to  200 
Forming replications  201  to  300 
Forming replications  301  to  400 
Forming replications  401  to  500 
Forming replications  501  to  600 
Forming replications  601  to  700 
Forming replications  701  to  800 
Forming replications  801  to  900 
Forming replications  901  to  1000 

> 
> summary(r.boot)
Call:
bootstrap(data = X, statistic = cor2)

Number of Replications: 1000 

Summary Statistics:
      Observed      Bias   Mean      SE 
cor21   0.8766 -0.003166 0.8735 0.04564
cor22   1.3610  0.018210 1.3792 0.19211

Empirical Percentiles:
        2.5%     5%   95%  97.5% 
cor21 0.7691 0.7931 0.935 0.9402
cor22 1.0180 1.0796 1.697 1.7394

BCa Confidence Limits:
        2.5%     5%   95%  97.5% 
cor21 0.7319 0.7619 0.928 0.9339
cor22 0.9114 0.9950 1.635 1.6779

Correlation of Replicates:
       cor21  cor22 
cor21 1.0000 0.9762
cor22 0.9762 1.0000
> plot(r.boot)
>