###################################################################################
# Date: 2014-04-30
# Time series analysis
# listening to Climate history
###################################################################################
rm(list=ls())
library(audio)

# interpolated 10-year resolution
# CO2 concentrations
buffer<-read.table("data/timeseries_co2_800000.asc")
time<-buffer[,1]
x<-buffer[,2]
n<-length(x)

plot(time,x,xlab="time [years]",
     ylab="CO2 [parts per million]",
     typ='l')

wait<-readline("press enter to listen to the time series")  
rate<-22050


play(x,rate=rate)