  We calculate the angular correlation function of galaxies in the Two Micron
All Sky Survey. We minimize the possible contamination by stars, dust, seeing
and sky brightness by studying their cross correlation with galaxy density, and
limiting the galaxy sample accordingly. We measure the correlation function at
scales between 1-18 arcdegs using a half million galaxies. We find a best fit
power law to the correlation function has a slope of 0.76 and an amplitude of
0.11. However, there are statistically significant oscillations around this
power law. The largest oscillation occurs at about 0.8 degrees, corresponding
to 2.8 h^{-1} Mpc at the median redshift of our survey, as expected in halo
occupation distribution descriptions of galaxy clustering.
  We invert the angular correlation function using Singular Value Decomposition
to measure the three-dimensional power spectrum and find that it too is in good
agreement with previous measurements. A dip seen in the power spectrum at small
wavenumber k is statistically consistent with CDM-type power spectra. A fit of
CDM-type power spectra to k < 0.2 h Mpc^{-1} give constraints of
\Gamma_{eff}=0.116 and \sigma_8=0.96. This suggest a K_s-band linear bias of
1.1+/-0.2. This \Gamma_{eff} is different from the WMAP CMB derived value. On
small scales the power-law shape of our power spectrum is shallower than that
derived for the SDSS. These facts together imply a biasing of these different
galaxies that might be nonlinear, that might be either waveband or luminosity
dependent, and that might have a nonlocal origin.
