めも：qplot(1/RESOL**3, CORR, data=m[m$SIGIOBS>0], color=IOBS/SIGIOBS) + geom_density2d()