[R-br] bootstrap em correlação

Olá pessoal boa tarde Tenho a seguinte situação.  Possuo um conjunto de dados caracterizados por x1,x2,x3. e gostaria de realizar o bootstrap para calcular a correlação e efetuar intervalos de confiança bootsrap. A linha de pensamento é essa. Suponhamos que  x1<-rnorm(100) x2<-rnorm(100) x3<-rnorm(100) dados<-data.frame(x1,x2,x3) rx1x2<-cor(x1,x2) rx1x3<-cor(x1,x3) rx2x3<-cor(x2,x3) theta1<- sqrt[(rx1x2+rx1x3)/rx2x3] theta2<- sqrt[(rx1x3+rx2x3)/rx2x1] Gostaria de ter como resultados final os intervalos bootstrap para theta1, theta2.  obrigado