xv=m[1]+s[1]*2
ycx=y[x>xv-0.01 & x<xv+0.01]
mu_ycx = m[2]+r*s[2]/s[1]*(xv-m[1])
var_ycx = s[2]^2*(1-r^2)
ybarcx = mean(ycx)
ycxs2 = var(ycx)

length(ycx)

mu_ycx
var_ycx
ybarcx
ycxs2