  We compare higher moments of baryon numbers measured at the RHIC heavy ion
collision experiments with those by the lattice QCD calculations. We employ the
canonical approach, in which we can access the real chemical potential regions
avoiding the sign problem. In the lattice QCD simulations, we study several
fits of the number density in the pure imaginary chemical potential, and
analyze how these fits affects behaviors at the real chemical potential. In the
energy regions between $\sqrt{s}_{NN}$=19.6 and 200 GeV, the susceptibility
calculated at $T/T_c=0.93$ is consistent with experimental data at $0 \le
\mu_B/T < 1.5$, while the kurtosis shows similar behavior with that of the
experimental data in the small $\mu_B/T$ regions $0 \le \mu_B/T < 0.3$. The
experimental data at $\sqrt{s}_{NN}=$ 11.5 shows quite different behavior. The
lattice result in the deconfinement region,$T/T_c=1.35$, is far from
experimental data.
