We reconstruct the large-scale initial density field from the distribution of galaxies observed by the Sloan Digital Sky Survey (SDSS). After adding the small-scale fluctuations to match the power spectrum to that of the standard LCDM model, we make a cosmological N-body simulation of formation of structures from the initial conditions. Properties of the objects formed in the simulation can be statistically compared with those of the observed SDSS galaxies. The simulation makes it possible to know the past history of evolution of objects located in different environments, and also gives us information on the environmental parameters that cannot be directly obtained observationally. It is hoped that this comparative study leads us to better understanding of formation and evolution of galaxies in conjunction with large-scale structures in the universe.