Estimated ${\widehat{f}}_{\mathrm{NL}}^{\mathrm{local}}$ (left) and its error bar $\sigma ({\widehat{f}}_{\mathrm{NL}}^{\mathrm{local}})$ (right) using different subsets of the PR4 simulations for linear correction and error bars, as a function of the lensing amplitude in the corresponding subset, determined using the integrated bispectrum estimator. Each subset consists of 100 maps drawn randomly from the 600 SEVEM PR4 simulations. Note that the expected lensing bias of 4.3 on the estimated ${\widehat{f}}_{\mathrm{NL}}^{\mathrm{local}}$ has not been subtracted.