Let i = 1, . . . , N index a simple random sample of units drawn from some large population. For each unit we observe the vector of regressors Xi and, for each of the N (N − 1) ordered pairs of units, an outcome Yij . The outcomes Yij and Ykl are independent if their indices are disjoint, but dependent otherwise (i.e., “dyadically dependent”). Let Wij = (X’i, X’j )’ ; using the sampled data we seek to construct a nonparametric estimate of the mean regression function g(Wij)=E[Yij | Xi, Xj].
We present two sets of results. First, we calculate lower bounds on the minimax risk for estimating the regression function at (i) a point and (ii) under the inﬁnity norm. Second, we calculate (i) pointwise and (ii) uniform convergence rates for the dyadic analog of the familiar Nadaraya-Watson (NW) kernel regression estimator. We show that the NW kernel regression estimator achieves the optimal rates suggested by our risk bounds when an appropriate bandwidth sequence is chosen. This optimal rate differs from the one available under iid data: the effective sample size is smaller and dw = dim(Wij ) inﬂuences the rate differently.