It is well-established that the current distibution along the axis of a self-resonant short vertical helical monopole is cosinusoidal, and not triangular, as in a vertical rod. It is the difference in the current distributions of these two vertical antennas of equal height that cause the differnce in radiation resiatances. Consequently, the gain of a helical monopole compared to a rod the same height is (4/pi)^2 = 1.621, which is about 2dB.
NEC-2: When automatically creating a long, narrow, helix for NEC analysis, care must be taken to avoid overly short segments for modeling the turns, because this can cause numerical instability in the simulation. This is just a suggestion. I don't actually know why the wrong results were obtained. Repeating the analysis with a fatter coil with fewer turns might give more realistic results.
John Kraus: From what I have seen, Kraus was mostly interested in analyzing the helical antenna operating in the axial mode. For helical antennas operating in the normal mode, he mostly deferred to the previous work by Wheeler. To simplify analysis, Kraus assumed additional end loading, or top loading, external to the helix. The additional end (or top) loading changes the current distribution along the axis of the helix, and therefore the calculated radiation resistance. A top-loaded helix has about the same radiation resistance as a top-loaded rod.