library(freqparcoord)
data(mlb)
m70 <- mlb[mlb$Height >= 70 & mlb$Height <= 77,]
sds <- tapply(m70$Weight,m70$Height,sd)
plot(70:77,sds,xlab='ht')