Statistical analyses of finite sample distributions usually assume that fluctuations are self-averaging, i.e. that they are statistically similar in different regions of the given sample volume. By using the scale-length method, we test whether this assumption is satisfied in several samples of the Sloan Digital Sky Survey Data Release Six. We find that the probability density function (PDF) of conditional fluctuations, filtered on large enough spatial scales (i.e., r>30 Mpc/h), shows relevant systematic variations in different sub-volumes of the survey. Instead for scales r<30 Mpc/h the PDF is statistically stable, and its first moment presents scaling behavior with a negative exponent around one. Thus while up to 30 Mpc/h galaxy structures have well-defined power-law correlations, on larger scales it is not possible to consider whole sample average quantities as meaningful and useful statistical descriptors. This situation is due to the fact that galaxy structures correspond to density fluctuations which are too large in amplitude and too extended in space to be self-averaging on such large scales inside the sample volumes: galaxy distribution is inhomogeneous up to the largest scales, i.e. r ~ 100 Mpc/h, probed by the SDSS samples. We show that cosmological corrections, as K-corrections and standard evolutionary corrections, do not qualitatively change the relevant behaviors. Finally we show that the large amplitude galaxy fluctuations observed in the SDSS samples are at odds with the predictions of the standard LCDM model of structure formation.(Abridged version).

We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the system size. Correspondingly, we find that the density fluctuations follow the Gumbel distribution of extreme value statistics. This distribution is clearly distinguishable from a Gaussian distribution, which would arise for a homogeneous spatial galaxy configuration. We also point out similarities between the galaxy distribution and critical systems of statistical physics.

One of the most striking features predicted by standard models of galaxy formation is the presence of anti-correlations in the matter distribution at large enough scales (r>r_c). Simple arguments show that the location of the length-scale r_c, marking the transition from positive to negative correlations, is the same for any class of objects as for the full matter distribution, i.e. it is invariant under biasing. This scale is predicted by models to be at about the same distance of the scale signaling the baryonic acoustic oscillation scale r_{bao}. We test these predictions in the newest SDSS galaxy samples.We find that, in several MG samples, the correlation function remains positive at scales >250 Mpc/h, while in the concordance LCDM it should be negative beyond r_c\approx 120 Mpc/h. In other samples the correlation function becomes negative at scales <50 Mpc/h. To investigate the origin of these differences we consider in detail the propagation of errors on the sample density into the estimation of the correlation function. We conclude that these are important at large enough separations, and that they are responsible for the observed differences between different estimators and for the measured sample to sample variations of the correlation function. We conclude that, in the newest SDSS samples, the large scale behavior of the galaxy correlation function is affected by intrinsic errors andv olume-dependent systematic effects which make the detection of correlations to be only an estimate of a lower limit of their amplitude, spatial extension and statistical errors. We point out that these results represent an important challenge to LCDM models as they largely differ from its predictions.(Abridged version).

The properties of the galaxy distribution at large scales are usually studied using statistics which are assumed to be self-averaging inside a given sample. We present a new analysis able to quantitatively map galaxy large scale structures while testing for the stability of average statistical quantities in different sample regions. We find that the newest samples of the Sloan Digital Sky Survey provide unambiguous evidence that galaxy structures correspond to large amplitude density fluctuations at all scales limited only by sample sizes. The two-point correlations properties are self-averaging up to approximately 30 Mpc/h and are characterized by a fractal dimension D=2.1 +- 0.1. Then at all larger scales probed density fluctuations are too large in amplitude and too extended in space to be self-averaging inside the considered volumes. These inhomogeneities are compatible with a continuation of fractal correlations but incompatible with: (i) a homogeneity scale smaller than 100 Mpc/h, (ii) predictions of standard theoretical models, (iii) mock galaxy catalogs generated from cosmological Nbody simuations.

We study statistical properties of galaxy structures in several samples extracted from the 2dF Galaxy Redshift Survey. In particular, we measured conditional fluctuations by means of the scale-length method and determined their probability distribution. In this way we find that galaxy distribution in these samples is characterized by large amplitude fluctuations with a large spatial extension, whose size is only limited by the sample's boundaries. These fluctuations are quite typical and persistent in the sample's volumes, and they are detected in two independent regions in the northern and southern galactic caps. We discuss the relation of the scale-length method to several statistical quantities, such as counts of galaxies as a function of redshift and apparent magnitude. We confirm previous results, which have determined by magnitude and redshift counts that there are fluctuations of about 30% between the southern and the northern galactic caps and we relate explicitly these counts to structures in redshift space. We show that the estimation of fluctuation amplitude normalized to the sample density is biased by systematic effects, which we discuss in detail. We consider the type of fluctuations predicted by standard cosmological models of structure formation in the linear regime and, to study nonlinear clustering, we analyze several samples of mock-galaxy catalogs generated from the distribution of dark matter in cosmological N-body simulations. In this way we conclude that the galaxy fluctuations present in these samples are too large in amplitude and too extended in space to be compatible with the predictions of the standard models of structure formation.

We apply the scale-length method to several three dimensional samples of the Two degree Field Galaxy Redshift Survey. This method allows us to map in a quantitative and powerful way large scale structures in the distribution of galaxies controlling systematic effects. By determining the probability density function of conditional fluctuations we show that large scale structures are quite typical and correspond to large fluctuations in the galaxy density field. We do not find a convergence to homogeneity up to the samples sizes, i.e. ~ 75 Mpc/h. We then measure, at scales r <~ 40 Mpc/h, a well defined and statistically stable power-law behavior of the average number of galaxies in spheres, with fractal dimension D=2.2 +- 0.2. We point out that standard models of structure formation are unable to explain the existence of the large fluctuations in the galaxy density field detected in these samples. This conclusion is reached in two ways: by considering the scale, determined by the linear perturbation analysis of a self-gravitating fluid, below which large fluctuations are expected in standard models and through the determination of statistical properties of mock galaxy catalogs generated from cosmological N-body simulations of the Millenium consortitum.

We discuss the large scale properties of standard cold dark matter cosmological models characterizing the main features of the power-spectrum, of the two-point correlation function and of the mass variance. Both the real-space statistics have a very well defined behavior on large enough scales, where their amplitudes become smaller than unity. The correlation function, in the range 0<\xi(r)<1, is characterized by a typical length-scale r_c, at which \xi(r_c)=0, which is fixed by the physics of the early universe: beyond this scale it becomes negative, going to zero with a tail proportional to -(r^{-4}). These anti-correlations represent thus an important observational challenge to verify models in real space. The same length scale r_c characterizes the behavior of the mass variance which decays, for r>r_c, as r^{-4}, the fastest decay for any mass distribution. The length-scale r_c defines the maximum extension of (positively correlated) structures in these models. These are the features expected for the dark matter field: galaxies, which represent a biased field, however may have differences with respect to these behaviors, which we analyze. We then discuss the detectability of these real space features by considering several estimators of the two-point correlation function. By making tests on numerical simulations we emphasize the important role of finite size effects which should always be controlled for careful measurements.

We discuss the estimation of galaxy correlation properties in several volume limited samples, in different sky regions, obtained from the Fourth Data Release of the Sloan Digital Sky Survey. The small scale properties are characterized through the determination of the nearest neighbor probability distribution. By using a very conservative statistical analysis, in the range of scales [0.5,~30] Mpc/h we detect power-law correlations in the conditional density in redshift space, with an exponent \gamma=1.0 \pm 0.1. This behavior is stable in all different samples we considered thus it does not depend on galaxy luminosity. In the range of scales [~30,~100] Mpc/h we find evidences for systematic unaveraged fluctuations and we discuss in detail the problems induced by finite volume effects on the determination of the conditional density. We conclude that in such range of scales there is an evidence for a smaller power-law index of the conditional density. However we cannot distinguish between two possibilities: (i) that a crossover to homogeneity (corresponding to \gamma=0 in the conditional density) occurs before 100 Mpc/h, (ii) that correlations extend to scales of order 100 Mpc/h (with a smaller exponent 0 < \gamma <1). We emphasize that galaxy distributions in these samples present large fluctuations at the largest scales probed, corresponding to the presence of large scale structures extending up to the boundaries of the present survey. Finally we discuss several differences between the behavior of the conditional density in mock galaxy catalogs built from cosmological N-body simulations and real data. We discuss some theoretical implications of such a fact considering also the super-homogeneous features of primordial density fields.

We study galaxy correlations from samples extracted from the 2dFGRS final release. Statistical properties are characterized by studying the nearest neighbor probability density, the conditional density and the reduced two-point correlation function. The result is that the conditional density has a power-law behavior in redshift space described by an exponent \gamma=0.8 \pm 0.2 in the interval from about 1 Mpc/h, the average distance between nearest galaxies, up to about 40 Mpc/h, corresponding to radius of the largest sphere contained in the samples. These results are consistent with other studies of the conditional density and are useful to clarify the subtle role of finite-size effects on the determination of the two-point correlation function in redshift and real space