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