Counts-in-cells are measured in the \tau CDM Virgo Hubble Volume simulation . This large N -body experiment has 10 ^ { 9 } particles in a cubic box of size 2000 h ^ { -1 } Mpc . The unprecedented combination of size and resolution allows for the first time a realistic numerical analysis of the cosmic errors and cosmic correlations of statistics related to counts-in-cells measurements , such as the probability distribution function P _ { N } itself , its factorial moments F _ { k } and the related cumulants \overline { \xi } and S _ { N } ’ s . These statistics are extracted from the whole simulation cube , as well as from 4096 sub-cubes of size 125 h ^ { -1 } Mpc , each representing a virtual random realization of the local universe . The measurements and their scatter over the sub-volumes are compared to the theoretical predictions of Colombi , Bouchet & Schaeffer ( 1995 ) for P _ { 0 } , and of Szapudi & Colombi ( 1996 , SC ) and Szapudi , Colombi & Bernardeau ( 1999a , SCB ) for the factorial moments and the cumulants . The general behavior of experimental variance and cross-correlations as functions of scale and order is well described by theoretical predictions , with a few percent accuracy in the weakly non-linear regime for the cosmic error on factorial moments . On highly non-linear scales , however , all variants of the hierarchical model used by SC and SCB to describe clustering appear to become increasingly approximate , which leads to a slight overestimation of the error , by about a factor of two in the worst case . Because of the needed supplementary perturbative approach , the theory is less accurate for non-linear estimators , such as cumulants , than for factorial moments . The cosmic bias is evaluated as well , and , in agreement with SCB , is found to be insignificant compared to the cosmic variance in all regimes investigated . While higher order statistics were previously evaluated in several simulations , this work presents text book quality measurements of S _ { N } ’ s , 3 \leq N \leq 10 , in an unprecedented dynamic range of 0.05 \la \overline { \xi } \la 50 . In the weakly nonlinear regime the results confirm previous findings and agree remarkably well with perturbation theory predictions including the one loop corrections based on spherical collapse by Fosalba & Gaztañaga 1998 . Extended perturbation theory is confirmed on all scales .