We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of deterministic linear biasing $\delta_{\rm obj}\simeq b_{\rm obj} \delta_{\rm mass}$ works reasonably well even in the quasi-linear regime $r > 10 h^{-1}$ Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime.

We present a theoretical formalism to predict the two-point clustering statistics (the power spectrum and the two-point correlation function), simultaneously taking account of the linear velocity distortion, the nonlinear velocity distortion (finger-of-god), the cosmological redshift-space distortion and the light-cone effect. To demonstrate the importance of these effects in exploring the clustering of objects at high redshifts, we show several model predictions for magnitude-limited surveys of galaxies and quasars. This methodology provides a quantitative tool to confront theoretical models against the upcoming precision data on clustering in the universe.

We examine the cosmological redshift-space distortion effect on the power spectrum of the objects at high-redshifts, which is an unavoidable observational contamination in general relativistic cosmology. In particular, we consider the nonlinear effects of density and velocity evolution using high-resolution N-body simulations in cold dark matter models. We find that the theoretical modeling on the basis of the fitting formulae of nonlinear density and velocity fields accurately describes the numerical results, especially in quasi-nonlinear regimes. These corrections for nonlinear effects are essential in order to use the the cosmological redshift-space distortion as a cosmological test. We perform a feasibility test to derive constraints from the future catalogues of high-redshift quasars using our theoretical modeling and results of N-body simulations. Applying the present methodology to the future data from on-going surveys of high-redshift galaxies and quasars will provide a useful tool to constrain a geometry of the universe.

The three-dimensional distribution of astronomical objects observed in redshift space significantly differs from the true distribution since the distance to each object cannot be determined by its redshift $z$ only; for $z \ll 1$ the peculiar velocity field contaminates the true recession velocity of the Hubble flow, while the true distance for objects at $z < 1$ sensitively depends on the (unknown and thus assumed) cosmological parameters. This hampers the effort to understand the true distribution of large-scale structure of the universe. In addition, all cosmological observations are carried out on a light-cone, the null hypersurface of an observer at $z=0$. This implies that their intrinsic properties and clustering statistics should change even within the survey volume. Therefore a proper comparison taking account of the light-cone effect is important to extract any cosmological information from redshift catalogues, especially for $z < 1$. We present recent theoretical development on the two effects -- the cosmological light-cone effect and the cosmological redshift-space distortion -- which should play key roles in observational cosmology in the 21st century.