  In this work, we have performed numerical simulations of primordial black
hole (PBH) formation in the Friedman-Lema\^itre-Robertson-Walker universe
filled by radiation fluid, introducing the local-type non-Gaussianity to the
primordial curvature fluctuation. We have compared the numerical results from
simulations with previous analytical estimations on the threshold value for PBH
formation done in the previous paper arXiv:2109.00791, particularly for
negative values of the non-linearity parameter $f_{\rm NL}$. Our numerical
results show the existence of PBH formation of (the so-called) type I also in
the case $f_{\rm NL} \lesssim -0.336$, which was not found in the previous
analytical expectations using the critical averaged compaction function. In
particular, although the universal value for the averaged critical compaction
function $\bar{\mathcal{C}}_{c}=2/5$ found previously in the literature is not
satisfied for all the profiles considered in this work, an alternative direct
analytical estimate has been found to be roughly accurate to estimate the
thresholds, which gives the value of the critical averaged density with a few
$\%$ deviation from the numerical one for $f_{\rm NL}\gtrsim -1$.
