  For the univariate current status and, more generally, the interval censoring
model, distribution theory has been developed for the maximum likelihood
estimator (MLE) and smoothed maximum likelihood estimator (SMLE) of the unknown
distribution function, see, e.g., [12], [7], [4], [5], [6], [10], [11] and [8].
For the bivariate current status and interval censoring models distribution
theory of this type is still absent and even the rate at which we can expect
reasonable estimators to converge is unknown. We define a purely discrete
plug-in estimator of the distribution function which locally converges at rate
n^{1/3} and derive its (normal) limit distribution. Unlike the MLE or SMLE,
this estimator is not a proper distribution function. Since the estimator is
purely discrete, it demonstrates that the n^{1/3} convergence rate is in
principle possible for the MLE, but whether this actually holds for the MLE is
still an open problem. If the cube root n rate holds for the MLE, this would
mean that the local 1-dimensional rate of the MLE continues to hold in
dimension 2, a (perhaps) somewhat surprising result. The simulation results do
not seem to be in contradiction with this assumption, however. We compare the
behavior of the plug-in estimator with the behavior of the MLE on a sieve and
the SMLE in a simulation study. This indicates that the plug-in estimator and
the SMLE have a smaller variance but a larger bias than the sieved MLE. The
SMLE is conjectured to have a n^{1/3}-rate of convergence if we use bandwidths
of order n^{-1/6}. We derive its (normal) limit distribution, using this
assumption. Finally, we demonstrate the behavior of the MLE and SMLE for the
bivariate interval censored data of [1], which have been discussed by many
authors, see e.g., [18], [3], [2] and [15].
