Observational Constraints on Higher Order Clustering up to $z\simeq 1
  Constraints on the validity of the hierarchical gravitational instability theory and the evolution of biasing are presented based upon measurements of higher order clustering statistics in the Deeprange Survey, a catalog of $\sim710,000$ galaxies with $I_{AB} \le 24$ derived from a KPNO 4m CCD imaging survey of a contiguous $4^{\circ} \times 4^{\circ}$ region. We compute the 3-point and 4-point angular correlation functions using a direct estimation for the former and the counts-in-cells technique for both. The skewness $s_3$ decreases by a factor of $\simeq 3-4$ as galaxy magnitude increases over the range $17 \le I \le 22.5$ ($0.1 \lesssim z \lesssim 0.8$). This decrease is consistent with a small {\it increase} of the bias with increasing redshift, but not by more than a factor of 2 for the highest redshifts probed. Our results are strongly inconsistent, at about the $3.5-4 \sigma$ level, with typical cosmic string models in which the initial perturbations follow a non-Gaussian distribution - such models generally predict an opposite trend in the degree of bias as a function of redshift. We also find that the scaling relation between the 3-point and 4-point correlation functions remains approximately invariant over the above magnitude range. The simplest model that is consistent with these constraints is a universe in which an initially Gaussian perturbation spectrum evolves under the influence of gravity combined with a low level of bias between the matter and the galaxies that decreases slightly from $z \sim 0.8$ to the current epoch. 
Istv\'an Szapudi (CITA), Marc Postman (STScI), Tod R. Lauer (NOAO),   and William Oegerle(Johns Hopkins)