Displaying 1 result from an estimated 1 matches for "stratasizemon".
Did you mean:
stratasizemone
2008 Dec 03
1
help on tapply using sample with differing sample-sizes
...ciqntu)
{
i <- 1
sdatavar <- sort(datavar)
sstratavar <- sort(statavar)
sweight <- 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 <- re...