Criei o seguinte código para calcular a variância através do método jackknife. Porém, ao invés de retirar um termo por vez, agora preciso retirar de 2 em 2. sei que é algo básico, mas não consegui pensar em como fazer isso, alguém poderia me dar alguma ideia?
A<- numeric(1000)
var1<- numeric(1000)
x<- data.frame(numeric(9))
amostra<- rnorm(10,100,10) #gerando amostra
for(i in 1:10){
x[i]<-c(amostra[-i]) #tirando um termo da amostra
}
var2<- numeric(10) #calculando as variancias dos vetores n-1
for(i in 1:10){
var2[i]<- var(x[,i])
}
Tn<- sum(var2) #soma das variancias de 3 termos
var1[j] <- var(amostra)# variancia da amstra de tamanho n
A[j]<-10*var1[j]-((9/10)*Tn) # variancia Jackknife