hohenTp<- read.table("C:/CLIM/HohenT.txt",header=T)         ###
attach(hohenTp) 

n<- length(Year); options(digits=3)
Wi<- dcly+jan+feb; Sp<- mar+apr+may  
Su<- jun+jul+aug;  Au<- sep+oct+nov    #no averaging necessary
Quar<- cbind(Wi,Sp,Su,Au)        #binding Wi,Sp,Su,Au together
                                       #Quar is a n x 4 matrix
"----Residuals from polynomials(4)-trend----"
Ja<- Year-1800;  Ja2<- Ja*Ja; Ja3<- Ja2*Ja; Ja4<- Ja2*Ja2
Quares<- Quar-predict(lm(Quar~Ja+Ja2+Ja3+Ja4))      #residuals
 
"1 | 2  refers to preceding | succeeding year"
Quares1<- cbind(Quares[(1:(n-1)),(1:4)])        
Quares2<- cbind(Quares[(2:n),(1:4)])            
 
Quares<- cbind(Quares1,Quares2)  #Quares is a (n-1) x 8 matrix
colnames(Quares)<- c("Wires1","Spres1","Sures1","Aures1",
                     "Wires2","Spres2","Sures2","Aures2")
cor(Quares)         #cross tabulation of pairwise correlations  


detach(hohenTp)                                             ###
rm(list=objects())                                          ###