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

#Number of cases with above-/below-average values
CondFrequ<- function(x,y){     #Conditional frequencies y|x
x0<- x[x > mean(x)]                 
yx0<- y[x > mean(x) & y > mean(y)]
x1<- x[x < mean(x)]
yx1<- y[x < mean(x) & y < mean(y)]
c("x>"=length(x0),"x>&y>"=length(yx0), 
               "y>|x>"= length(yx0)/length(x0),
  "x<"=length(x1),"x<&y<"=length(yx1),
               "y<|x<"=length(yx1)/length(x1))
}
#----------------------------------------------------------
"x=Dec -> y=Jan"
x<- dcly/10; y <- jan/10  
c("mean(x)"=mean(x),"mean(y)"=mean(y),"cor(x,y)"=cor(x,y))
CondFrequ(x,y)

"x=Sep -> y=Oct"
x<- sep/10;  y <- oct/10
c("mean(x)"=mean(x),"mean(y)"=mean(y),"cor(x,y)"=cor(x,y))
CondFrequ(x,y)

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