Distribution of galaxies may be a biased tracer of the dark matter distribu
tion and the relation between the galaxies and the total mass may be stocha
stic, nonlinear and time dependent. Since many observations of galaxy clust
ering will be done at high redshift, the time evolution of nonlinear stocha
stic biasing would play a crucial role for the data analysis of the future
sky surveys. Recently, analytic study of the time evolution induced by grav
ity has been reported in the mildly nonlinear regime. Here, we further deve
lop the nonlinear analysis including the next-to-leading order and attempt
to clarify the nonlinear feature of the stochastic biasing. Employing the p
erturbative approach, we compute the one-loop correction of the power spect
rum for the total mass, galaxies, and their cross correlation. Assuming tha
t the initial distribution of galaxies is given by the local function, we s
pecifically investigate the time evolution of the biasing parameter and the
correlation coefficient deduced from the power spectra. On large scales, w
e find that the time evolution of the biasing parameter could deviate from
the linear theory prediction in presence of the initial skewness, even thou
gh the scale dependence of the biasing is very weak. On the other hand, the
deviation can be reduced if the stochasticity between the galaxies and the
total mass exists. To explore the influence of nonlinear gravity, we focus
on the quasi-linear scales, where the nonlinear growth of the total mass b
ecomes important. It is recognized that the scale dependence of the biasing
dynamically appears and the initial stochasticity could affect the time ev
olution of the scale dependence. The result is compared with the recent N-b
ody simulation that the scale dependence of the halo biasing can appear on
relatively large scales and the biasing parameter takes the lower value on
smaller scales. Qualitatively, our weakly nonlinear results can explain thi
s trend if the halo-mass biasing relation has the large scatter at high red
shift.