Displaying 1 result from an estimated 1 matches for "vectobsres".
2008 Dec 03
1
help on tapply using sample with differing sample-sizes
...sort(weight)
sdatavarwght <- sdatavar*sweight
# stramn <- tapply(sdatavar, sstratavar, mean)
meanvect <- rep(0, times = nboot)
stratasize <- tapply(sdatavarwght, sstratavar, length)
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 <-...