Structure formation from non-Gaussian initial conditions: Multivariate biasing, statistics, and comparison with N -body simulations