150 likes | 402 Views
Seasonal Analysis. UNR * STAT 758 * Fall 2006. What is the period?. R-code. par(bg='yellow') plot(co2,lwd=2,col=4,xlab='Time, years',ylab='CO2 in parts per million (ppm)', main='Mauna Loa Atmospheric CO2 Concentration') abline(v=seq(1950,2006,by=10),lty=2)
E N D
Seasonal Analysis UNR * STAT 758 * Fall 2006
What is the period? R-code par(bg='yellow') plot(co2,lwd=2,col=4,xlab='Time, years',ylab='CO2 in parts per million (ppm)', main='Mauna Loa Atmospheric CO2 Concentration') abline(v=seq(1950,2006,by=10),lty=2) abline(h=seq(320,360,by=10),lty=2) UNR * STAT 758 * Fall 2006
R-code # This line adds yearly grid abline(v=seq(1950,2006,by=1),lty=2) UNR * STAT 758 * Fall 2006
R-code # Filter with 1-year window and plot a trend f<-filter(co2,rep(1,12)/12) points(f,col=2,type='l',lwd=2) UNR * STAT 758 * Fall 2006
R-code # Plot series minus trend plot(co2-f,ldw=2,col=4,lwd=2, xlab='Time, years',ylab='CO2-trend',main='Detrended CO2 series') abline(h=seq(-4,4,by=2),lty=2) abline(v=seq(1950,2006,by=1),lty=2) UNR * STAT 758 * Fall 2006
R-code # Perform seasonal analysis and plot season component m<decompose(co2) plot(m$figure,lwd=2,type='p',pch=19,col=4, xlab='Month',ylab='Seasonal component') axis(1,at=seq(1,12)) grid(col='black',lty=2) UNR * STAT 758 * Fall 2006
R-code m<decompose(co2) plot(m) UNR * STAT 758 * Fall 2006
Moving Average UNR * STAT 758 * Fall 2006
Estimated signal UNR * STAT 758 * Fall 2006
Differencing UNR * STAT 758 * Fall 2006