set.seed(2005)
x <- round(rnorm(1000, 100, 10))
xscale <- 70:130
xtab <- freqtab(xscale, x)
xavg <- freqavg(xtab)
xbump <- freqbump(xtab)
cbind(xtab, avg = xavg,
bump = round(xbump * sum(xtab[, 2]), 6))
Run the code above in your browser using DataLab