[R-br] bootstrap em correlação

[geovane barbosa]

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 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140213/38514a0b/attachment.html>