rcollyer pointed out in a comment that the the new GroupBy may be substituted for GatherBy in Szabolcs's original to produce the desired function:
cleanPosIdx[x_] := GroupBy[Range @ Length @ x, x[[#]] &]
I shall be using this code until PositionIndex receives and enhancement.