data(battingLabels)
str(battingLabels)
require(plyr)
# find and plot maximum number of homers per year
batHR <- ddply(subset(Batting, !is.na(HR)), .(yearID),
summarise, max=max(HR))
with(batHR, {
plot(yearID, max,
xlab=Label("yearID"), ylab=paste("Maximum", Label("HR")),
cex=0.8)
lines(lowess(yearID, max), col="blue", lwd=2)
abline(lm(max ~ yearID), col="red", lwd=2)
})
Run the code above in your browser using DataLab