Statistics for Quijote-png: Quasi-maximum likelihood estimation of primordial non-gaussianity in the nonlinear dark matter density field