Displaying 1 result from an estimated 1 matches for "repvectboot".
2008 Dec 03
1
help on tapply using sample with differing sample-sizes
...th)
stratasizemone <- stratasize -1
while (i < nboot)
{ #vector of resampled observations
vectobsrestemp <- tapply(sdatavarwght, sstratavar, sample, replace = T, size = stratasizemone)
vectobsres <- unlist(vectobsrestemp)
meanvect[i] <- mean(vectobsres)
i <- i + 1
}
repvectboot <- rep(mean(meanvect), times = i)
vb <- sum((repvectboot - meanvect)^2)/(i -1)
lwrbndmnp <- quantile(meanvect, ciqntl)
uppbndmnp <- quantile(meanvect, ciqntu)
vb
}