output(T) rho.txt
pon
smpl 94.4 99.3
ls d(food) sar(12) sma(12)
!n=@regobs

@ RHO 1 ON D(FOOD)
genr dfood=d(food)
stom(dfood,vdfood)
genr dfood1=d(food(-1))
stom(dfood1,vdfood1)
!x=0
!z=0
for !i=1 to !n
        !x=!x+vdfood(!i)
        !z=!z+vdfood1(!i)
next
!x=!x/!n
!z=!z/!n
!x1=0
!z1=0
!c=0
!c2=0
!c3=0

for !i=1 to !n
        !x1=(vdfood(!i)-!x)^2+!x1
        !z1=(vdfood1(!i)-!z)^2+!z1
        !c=(vdfood(!i)-!x)*(vdfood1(!i)-!z)+!c
        !c2=(vdfood(!i)-!x)*(vdfood1(!i)-!x)+!c2
        !c3=(vdfood(!i)-!z)*(vdfood1(!i)-!z)+!c3
next
vector(2) var
var(1)=!x1/!n
var(2)=!z1/!n
show var
vector(6) rho
rho(1)=(!c/!n)/var(1)
rho(2)=(!c/!n)/var(2)
rho(3)=(!c2/!n)/var(1)
rho(4)=(!c2/!n)/var(2)
rho(5)=(!c3/!n)/var(1)
rho(6)=(!c3/!n)/var(2)
show rho
cor dfood dfood1



@ RHO2 ON D(FOOD)
genr dfood=d(food)
stom(dfood,vdfood)
genr dfood1=d(food(-2))
stom(dfood1,vdfood1)
!x=0
!z=0
for !i=1 to !n
        !x=!x+vdfood(!i)
        !z=!z+vdfood1(!i)
next
!x=!x/!n
!z=!z/!n
!x1=0
!z1=0
!c=0
!c2=0
!c3=0
for !i=1 to !n
        !x1=(vdfood(!i)-!x)^2+!x1
        !z1=(vdfood1(!i)-!z)^2+!z1
        !c=(vdfood(!i)-!x)*(vdfood1(!i)-!z)+!c
        !c2=(vdfood(!i)-!x)*(vdfood1(!i)-!x)+!c2
        !c3=(vdfood(!i)-!z)*(vdfood1(!i)-!z)+!c3
next
vector(2) var
var(1)=!x1/!n
var(2)=!z1/!n
show var
vector(6) rho
rho(1)=(!c/!n)/var(1)
rho(2)=(!c/!n)/var(2)
rho(3)=(!c2/!n)/var(1)
rho(4)=(!c2/!n)/var(2)
rho(5)=(!c3/!n)/var(1)
rho(6)=(!c3/!n)/var(2)
show rho
cor dfood dfood1

 