#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

# temperature measurements: 
Z = c( 200, 202, 208, 204, 204, 207, 207, 204, 202, 199, 201, 198, 200,
       202, 203, 205, 207, 211, 204, 206, 203, 203, 201, 198, 200, 206,
       207, 206, 200, 203, 203, 200, 200, 195, 202, 204 )
Z = as.double(Z)
N = length(Z)

# Extract the numerical values:
# 
z_covariance  = acf( Z, type="covariance" )
z_correlation = acf( Z, type="correlation" ) 

# Plot the autocorrelation:
# 
postscript("../../WriteUp/Graphics/Chapter2/prob_4_plot_autocorrelation.eps", onefile=FALSE, horizontal=FALSE)
acf( Z, type="correlation" ) 
dev.off()