We review some of the common methods for model selection: the goodness of fit, the likelihood ratio test, Bayesian model selection using Bayes factors, and the classical as well as the Bayesian information theoretic approaches. We illustrate these different approaches by comparing models for the expansion history of the Universe. In the discussion we highlight the premises and objectives entering these different approaches to model selection and finally recommend the information theoretic approach.

Properties of galaxies like their absolute magnitude and their stellar mass content are correlated. These correlations are tighter for close pairs of galaxies, which is called galactic conformity. In hierarchical structure formation scenarios, galaxies form within dark matter halos. To explain the amplitude and the spatial range of galactic conformity two--halo terms or assembly bias become important. With the scale dependent correlation coefficients the amplitude and the spatial range of conformity are determined from galaxy and halo samples. The scale dependent correlation coefficients are introduced as a new descriptive statistic to quantify the correlations between properties of galaxies or halos, depending on the distances to other galaxies or halos. These scale dependent correlation coefficients can be applied to the galaxy distribution directly. Neither a splitting of the sample into subsamples, nor an a priori clustering is needed. This new descriptive statistic is applied to galaxy catalogues derived from the Sloan Digital Sky Survey III and to halo catalogues from the MultiDark simulations. In the galaxy sample the correlations between absolute Magnitude, velocity dispersion, ellipticity, and stellar mass content are investigated. The correlations of mass, spin, and ellipticity are explored in the halo samples. Both for galaxies and halos a scale dependent conformity is confirmed. Moreover the scale dependent correlation coefficients reveal a signal of conformity out to 40Mpc and beyond. The halo and galaxy samples show a differing amplitude and range of conformity.

To study the global morphology of the galaxy distribution in our local neighbourhood we calculate the Minkowski functionals for a sequence of volume limited samples within a sphere of 8 Mpc centred on our galaxy. The well known strong clustering of the galaxies and the dominance of voids and coherent structures on larger scales is clearly visible in the Minkowski functionals. The morphology of the galaxy distribution changes with the limiting absolute magnitude. The samples, encompassing the more luminous galaxies, show emptier voids and more pronounced coherent structures. Indeed there is a prominent peak in the luminosity function of isolated galaxies for MB approx. -14, which at least partly explains these morphological changes. We compare with halo samples from a LambdaCDM simulation. Special care was taken to reproduce the observed local neighbourhood as well as the observed luminosity function in these mock samples. All in all the mock samples render the global morphology of the galaxy distribution quite well. However the detailed morphological analysis reveals that real galaxies cluster stronger, the observed voids are emptier and the structures are more pronounced compared to the mock samples from the LambdaCDM simulation.

In this letter I present results from a correlation analysis of three galaxy redshift catalogs: the SSRS2, the CfA2 and the PSCz. I will focus on the observation that the amplitude of the two--point correlation function rises if the depth of the sample is increased. There are two competing explanations for this observation, one in terms of a fractal scaling, the other based on luminosity segregation. I will show that there is strong evidence that the observed growth is due to a luminosity dependent clustering of the galaxies.

The relative orientation of clusters' major elongation axes and clusters' angular momentum is studied using a large N-body simulation in a box of 500 Mpc/h base length for a standard LambdaCDM model. Employing the technique of mark correlation functions, we successfully separated the correlations in the orientation from the well known clustering signal traced by the two-point correlation function. The correlations in the orientation are highly significant for our sample of 3000 clusters. We found an alignment of neighboring clusters, i.e. an enhanced probability of the major elongation axes of neighboring cluster pairs to be in parallel with each other. At 10 Mpc/h separation the amplitude of this signal is ~10% above the value expected from random orientations, and it vanishes on scales larger than 15 Mpc/h. The ``filamentary'' alignment between clusters' major elongation axes and the lines pointing towards neighboring clusters shows even stronger deviations from random orientation, which can be detected out to scales of 100 Mpc/h, both in 2D and 3D analyses. Similarly, strong correlations of the angular momentum were seen. Also a clear signal in the scalar correlation of the absolute value of the angular momentum, the spin parameter and the mass was found. They extend up to 50 Mpc/h and have an amplitude of 40%, 15%, and 10% above a random distribution at 10 Mpc/h separation, respectively.

We use a novel statistical tool, the mark correlation functions (MCFs), to study clustering of galaxy-size halos as a function of their properties and environment in a high-resolution numerical simulation of the LambdaCDM cosmology. We applied MCFs using several types of continuous and discrete marks: maximum circular velocity of halos, merger mark indicating whether halos experienced or not a major merger in their evolution history, and a stripping mark indicating whether the halo underwent a tidal stripping. We find that halos which experienced a relatively early (z>1) major merger or mass loss (due to tidal stripping) in their evolution histories are over-abundant in halo pairs with separations < 3 Mpc/h. This result can be interpreted as spatial segregation of halos with different merger histories, qualitatively similar to the morphological segregation in the observed galaxy distribution. The analysis presented in this paper demonstrate that MCFs provide powerful, yet algorithmically simple, quantitative measures of segregation in the spatial distribution of objects with respect to their various properties (marks).

Mark correlations provide a systematic approach to look at objects both distributed in space and bearing intrinsic information, for instance on physical properties. The interplay of the objects' properties (marks) with the spatial clustering is of vivid interest for many applications; are, e.g., galaxies with high luminosities more strongly clustered than dim ones? Do neighbored pores in a sandstone have similar sizes? How does the shape of impact craters on a planet depend on the geological surface properties? In this article, we give an introduction into the appropriate mathematical framework to deal with such questions, i.e. the theory of marked point processes. After having clarified the notion of segregation effects, we define universal test quantities applicable to realizations of a marked point processes. We show their power using concrete data sets in analyzing the luminosity-dependence of the galaxy clustering, the alignment of dark matter halos in gravitational $N$-body simulations, the morphology- and diameter-dependence of the Martian crater distribution and the size correlations of pores in sandstone. In order to understand our data in more detail, we discuss the Boolean depletion model, the random field model and the Cox random field model. The first model describes depletion effects in the distribution of Martian craters and pores in sandstone, whereas the last one accounts at least qualitatively for the observed luminosity-dependence of the galaxy clustering.

The definition and the properties of a Gaussian point distribution, in contrast to the well-known properties of a Gaussian random field are discussed. Constraints for the number density and the two-point correlation function arise. A simple method for the simulation of this so-called Gauss-Poisson point process is given and illustrated with an example. The comparison of the distribution of galaxies in the PSCz catalogue with the Gauss-Poisson process underlines the importance of higher-order correlation functions for the description for the galaxy distribution. The construction of the Gauss-Poisson point process is extended to the n-point Poisson cluster process, now incorporating correlation functions up to the nth-order. The simulation methods and constraints on the correlation functions are discussed for the n-point case and detailed for the three-point case. As another approach, well suited for strongly clustered systems, the generalized halo-model is discussed. The influence of substructure inside the halos on the two- and three-point correlation functions is calculated in this model.

The redshift dependence of the abundance of collapsed objects places strong constraints on cosmological models of structure formation. We apply a recently proposed model describing the anisotropic collapse of inhomogeneous spatial domains. Compared with the spherical top-hat model, this generic model leads to significantly more collapsed objects at high redshifts: at redshift one and on the scale of rich clusters a factor of 65. Furthermore, for a fixed normalization of the initial fluctuation spectrum (\sigma_8=1), we predict four times as much presently collapsed objects on the mass-scale of rich clusters within the standard CDM cosmogony, compared to the spherical collapse.

In order to quantify higher-order correlations of the galaxy cluster distribution we use a complete family of additive measures which give scale-dependent morphological information. Minkowski functionals can be expressed analytically in terms of integrals of n-point correlation functions. They can be compared with measured Minkowski functionals of volume limited samples extracted from the REFLEX survey. We find significant non-Gaussian features in the large-scale spatial distribution of galaxy clusters. A Gauss-Poisson process can be excluded as a viable model for the distribution of galaxy clusters at the significance level of 95%.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

In a follow-up study to a previous analysis of the IRAS 1.2Jy catalogue, we quantify the morphological fluctuations in the PSCz survey. We use a variety of measures, among them the family of scalar Minkowski functionals. We confirm the existence of significant fluctuations that are discernible in volume-limited samples out to 200Mpc/h. In contrast to earlier findings, comparisons with cosmological N-body simulations reveal that the observed fluctuations roughly agree with the cosmic variance found in corresponding mock samples. While two-point measures, e.g. the variance of count-in-cells, fluctuate only mildly, the fluctuations in the morphology on large scales indicate the presence of coherent structures that are at least as large as the sample.

The results are still valid, however this letter is superceeded by a significantly extended version which is available at astro-ph/0003358 and is scheduled for publication in the ApJ Dec. 10, 2000. We have withdrawn this letter to avoid confusion.

How does the clustering of galaxies depend on their inner properties like morphological type and luminosity? We address this question in the mathematical framework of marked point processes and clarify the notion of luminosity and morphological segregation. A number of test quantities such as conditional mark-weighted two-point correlation functions are introduced. These descriptors allow for a scale-dependent analysis of luminosity and morphology segregation. Moreover, they break the degeneracy between an inhomogeneous fractal point set and actual present luminosity segregation. Using the Southern Sky Redshift Survey~2 (da Costa et al. 1998, SSRS2) we find both luminosity and morphological segregation at a high level of significance, confirming claims by previous works using these data (Benoist et al. 1996, Willmer et al. 1998). Specifically, the average luminosity and the fluctuations in the luminosity of pairs of galaxies are enhanced out to separations of 15Mpc/h. On scales smaller than 3Mpc/h the luminosities on galaxy pairs show a tight correlation. A comparison with the random-field model indicates that galaxy luminosities depend on the spatial distribution and galaxy-galaxy interactions. Early-type galaxies are also more strongly correlated, indicating morphological segregation. The galaxies in the PSCz catalog (Saunders et al. 2000) do not show significant luminosity segregation. This again illustrates that mainly early-type galaxies contribute to luminosity segregation. However, based on several independent investigations we show that the observed luminosity segregation can not be explained by the morphology-density relation alone.

Averaging and evolving inhomogeneities are non-commuting operations. This implies the existence of deviations of an averaged model from the standard Friedmann-Lemaitre cosmologies. We quantify these deviations, encoded in a backreaction parameter, in the framework of Newtonian cosmology. We employ the linear theory of gravitational instability in the Eulerian and Lagrangian approaches, as well as the spherically- and plane-symmetric solutions as standards of reference. We propose a model for the evolution of the average characteristics of a spatial domain for generic initial conditions that contains the spherical top-hat model and the planar collapse model as exact sub cases. A central result is that the backreaction term itself, calculated on sufficiently large domains, is small but, still, its presence can drive the cosmological parameters on the averaging domain far away from their global values of the standard model. We quantify the variations of these parameters in terms of the fluctuations in the initial data as derived from the power spectrum of initial cold dark matter density fluctuations. E.g. in a domain with a radius of 100Mpc today and initially one-sigma fluctuations, the density parameters deviate from their homogeneous values by 15%; three-sigma fluctuations lead to deviations larger than 100%.

Nine of the most important estimators known for the two-point correlation function are compared using a predetermined, rigorous criterion. The indicators were extracted from over 500 subsamples of the Virgo Hubble Volume simulation cluster catalog. The ``real'' correlation function was determined from the full survey in a 3000Mpc/h periodic cube. The estimators were ranked by the cumulative probability of returning a value within a certain tolerance of the real correlation function. This criterion takes into account bias and variance, and it is independent of the possibly non-Gaussian nature of the error statistics. As a result for astrophysical applications a clear recommendation has emerged: the Landy & Szalay (1993) estimator, in its original or grid version Szapudi & Szalay (1998), are preferred in comparison to the other indicators examined, with a performance almost indistinguishable from the Hamilton (1993) estimator.

Methods for the statistical characterization of the large-scale structure in the Universe will be the main topic of the present text. The focus is on geometrical methods, mainly Minkowski functionals and the J-function. Their relations to standard methods used in cosmology and spatial statistics and their application to cosmological datasets will be discussed. This work is not only meant as a short review for comologist, but also attempts to illustrate these morphological methods and to make them accessible to scientists from other fields. Consequently, a short introduction to the standard picture of cosmology is given.

Most cosmological models studied today are based on the assumption of homogeneity and isotropy. Observationally one can find evidence that supports these assumptions on very large scales, the strongest being the almost isotropy of the Cosmic Microwave Background radiation after assigning the whole dipole to our proper motion relative to this background. However, on small and on intermediate scales up to several hundreds of Mpcs, there are strong deviations from homogeneity and isotropy. Here the problem arises how to relate the observations with the homogeneous and isotropic models. The usual proposal for solving this problem is to assume that Friedmann-Lemaitre models describe the mean observables. Such mean values may be identified with spatial averages. For Newtonian fluid dynamics the averaging procedure has been discussed in detail in Buchert and Ehlers (1997), leading to an additional backreaction term in the Friedmann equation. We use the Eulerian linear approximation and the `Zel'dovich approximation' to estimate the effect of the backreaction term on the expansion. Our results indicate that even for domains matching the background density in the mean, the evolution of the scale factor strongly deviates from the Friedmann solution, critically depending on the velocity field inside.

Second-order measures, such as the two-point correlation function, are geometrical quantities describing the clustering properties of a point distribution. In this article well-known estimators for the correlation integral are reviewed and their relation to geometrical estimators for the two-point correlation function is put forward. Simulations illustrate the range of applicability of these estimators. The interpretation of the two-point correlation function as the excess of clustering with respect to Poisson distributed points has led to biases in common estimators. Comparing with the approximately unbiased geometrical estimators, we show how biases enter the estimators introduced by Davis and Peebles, Landy and Szalay, and Hamilton. We give recommendations for the application of the estimators, including details of the numerical implementation. The properties of the estimators of the correlation integral are illustrated in an application to a sample of IRAS galaxies. It is found that, due to the limitations of current galaxy catalogues in number and depth, no reliable determination of the correlation integral on large scales is possible. In the sample of IRAS galaxies considered, several estimators using different finite-size corrections yield different results on scales larger than 20 Mpc/h, while all of them agree on smaller scales.

We present the function J as a morphological descriptor for point patterns formed by the distribution of galaxies in the Universe. This function was recently introduced in the field of spatial statistics, and is based on the nearest neighbor distribution and the void probability function. The J descriptor allows to distinguish clustered (i.e. correlated) from ``regular'' (i.e. anti-correlated) point distributions. We outline the theoretical foundations of the method, perform tests with a Matern cluster process as an idealised model of galaxy clustering, and apply the descriptor to galaxies and loose groups in the Perseus-Pisces Survey. A comparison with mock-samples extracted from a mixed dark matter simulation shows that the J descriptor can be profitably used to constrain (in this case reject) viable models of cosmic structure formation.

We use a measure of clustering derived from the nearest neighbour distribution and the void probability function to distinguish between regular and clustered structures. This measure offers a succinct way to incorporate additional information beyond the two-point correlation function. Application to a supercluster catalogue by Einasto et al. (1997) reveals no clustering in the distribution of superclusters. However, we show that this supercluster catalogue is severely affected by construction effects. Taking these biases into account we still find some indications for regularity on the largest scales, but the significance is only at the one-sigma level.

An analysis of the IRAS 1.2 Jy redshift catalogue with emphasis on the separate examination of northern and southern parts (in galactic coordinates) is performed using a complete set of morphological descriptors (Minkowski functionals), nearest neighbour distributions and the variance of the galaxy counts. We find large fluctuations in the clustering properties as seen in a large difference between the northern and southern parts of the catalogue on scales of 100Mpc/h. These fluctuations remain discernible even on the scale of 200Mpc/h. We also identify sparse sampling as a major source of ``apparent homogenization''. Tests on observational selection effects concerning luminosity, colour and redshift-space distortion support the significance of these large-scale fluctuations.

Using a measure of clustering derived from the nearest neighbour distribution and the void probability function we are able to distinguish between regular and clustered structures. With an example we show that regularity is a property of a point set, which may be invisible in the two point correlation function. Applying this measure to a supercluster catalogue we conclude that there is some evidence for regular structures on large scales.

In an analysis of the IRAS 1.2 Jy redshift catalogue, with emphasis on the separate examination of northern and southern parts (in galactic coordinates), we found that the clustering of galaxies differs significantly between north and south, showing fluctuations in the clustering properties at least on scales of 100 Mpc/h (Kerscher et al. 1997, astro-ph/9704028). We give a brief description of our morphological method which is based on Minkowski functionals and show the results obtained from the IRAS 1.2 Jy galaxy catalogue. We discuss several error estimates and select different subsamples from the 1.2 Jy catalogue according to flux and colour and validate the results for the whole sample. Furthermore we look closer at the spatial origin of the fluctuation and we compare with optically selected galaxies from the CfA1 survey.

We determine the Minkowski functionals for a sample of Abell/ACO clusters, 401 with measured and 16 with estimated redshifts. The four Minkowski functionals (including the void probability function and the mean genus) deliver a global description of the spatial distribution of clusters on scales from $10$ to $60\hMpc$ with a clear geometric interpretation. Comparisons with mock catalogues of N--body simulations using different variants of the CDM model demonstrate the discriminative power of the description. The standard CDM model and the model with tilted perturbation spectrum cannot generate the Minkowski functionals of the cluster data, while a model with a cosmological constant and a model with breaking of the scale invariance of perturbations (BSI) yield compatible results.

We discuss the application of Minkowski functionals to galaxy catalogues with emphasis on the boundary corrections and point out the relevance for other statistical methods. The analysis of a cubic subsample of the CfA--1 data serves as an example.

Minkowski functionals provide a novel tool to characterize the large-scale galaxy distribution in the Universe. Here we give a brief tutorial on the basic features of these morphological measures and indicate their practical application for simulation data and galaxy redshift catalogues as examples.