Estimation of a two-dimensional distribution function under association