Flux-limited and volume-limited galaxy samples are constructed from SDSS data releases DR4, DR6 and DR7 for statistical analysis. The two-point correlation functions $\xi(s)$, monopole of three-point correlation functions $\zeta_0$, projected two-point correlation function $w_p$ and pairwise velocity dispersion $\sigma_{12}$ are measured to test if galaxy samples are fair for these statistics. We find that with increment of sky coverage of SDSS, $\xi(s)$ of flux-limited sample is extremely robust and insensitive to local structures at low redshift. But for volume-limited samples fainter than $L^*$ at large scales $s>\sim 10\hmpc$, deviation of $\xi(s)$ and $\zeta_0$ of DR7 to those of DR4 and DR6 increases with larger absolute magnitude. In the weakly nonlinear regime, there is no agreement between $\zeta_0$ of different data releases in all luminosity bins. Furthermore, $w_p$ of volume-limited samples of DR7 in luminosity bins fainter than $-M_{r,0.1}=[18.5,19.5]$ are significantly larger, and $\sigma_{12}$ of the two faintest volume-limited samples of DR7 display very different scale dependence than results of DR4 and DR6. Our findings call for cautions in understanding clustering analysis results of SDSS faint galaxy samples, and higher order statistics of SDSS volume-limited samples in the weakly nonlinear regime. The first zero-crossing points of $\xi(s)$ of volume-limited samples are also investigated and discussed.

The zeroth-order component of the cosine expansion of the projected three-point correlation function is proposed for clustering analysis of cosmic large scale structure. These functions are third order statistics but can be measured similarly to the projected two-point correlations. Numerical experiments with N-body simulations indicate that the advocated statistics are redshift distortion free within 10% in the non-linear regime on scales ~0.2-10Mpc/h. Halo model prediction of the zeroth-order component of the projected three-point correlation function agrees with simulations within ~10%. This lays the ground work for using these functions to perform joint analyses with the projected two-point correlation functions, exploring galaxy clustering properties in the framework of the halo model and relevant extensions.

We present detailed clustering analysis of a large K-band selected local galaxy sample, which is constructed from the 2MASS and the SDSS and consists of $82,486$ galaxies with $10 < K < 13.5$ and $0.01 < z < 0.1$. The two-point correlation function of the magnitude-limited sample in real space at small scales is well described by a power law $\xi(r)=(r/6.44\pm0.23)^{-1.81\pm0.02}$. The pairwise velocity dispersion is derived from the anisotropic two-point correlation function and find the dispersion $\sigma_{12}=685\pm 17\kms$ if its scale invariance is assumed, which is larger than values measured in optical bands selected galaxy samples. We further investigate the dependence of the two-point correlation function and the $\sigma_{12}$ on the $g-r$ color and the $K$-band luminosity, obtain similar results to previous works in optical bands. Comparing a mock galaxy sample with our real data indicates that the semi-analytical model can not mimic the $\sigma_{12}$ in observation albeit it can approximate the two-point correlation function within measurement uncertainties.