Tarefa para vc: entender os detalhes do que acontece no snippet abaixo.

f = function(dat, i){
cc = cor(dat[i,])
c(t1=sum(cc[1, -1])/cc[2, 3], t2=sum(cc[3, -3])/cc[2, 1])
}

set.seed(1)

dados = data.frame(x1=rnorm(100), x2=rnorm(100), x3=rnorm(100))

library(boot)

br = boot(dados, f, 10000)

boot.ci(br, index=1)

boot.ci(br, index=2)

Em 13 de fevereiro de 2014 15:19, geovane barbosa <geovanecb@yahoo.com.br> escreveu:

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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.