plot(visitors, col="lightgray", ylim=c(0,700000), ylab="Visitors")
par(new=TRUE);
plot(SMA(visitors, 30), ylim=c(0,700000), type='l', col="black", ylab="Visitors")

#set up timeseries, ignore first year when visitor counts weren't reliable at all sites
visitors = ts(all$visitor[365:length(all$visitor)], start=1, frequency=365.25);
#add a fictional visitor to each zero day so we can take logarithms
visitors[which(visitors==0)]=1;
#plot decomposition
attribs = decompose(visitors);
plot(attribs);

layout(c(2,1))
acf(attribs$random, na.action=na.pass);
    pacf(attribs$random, na.action=na.pass);

plot(stl(visitors, s.window="periodic"))