# NOT RUN {
# Load in sample data frame
data(unidata)
# Get data from ID number 21005, Saturday only
counts.saturday <-
unidata[unidata[, "seqn"] == 21005 & unidata[, "paxday"] == 7, "paxinten"]
# Calculate and plot hourly count averages
hourly.averages <- blockaves(x = counts.saturday, window = 60)
plot(hourly.averages)
# }
Run the code above in your browser using DataLab