Statistics of h- and p-adaptivity driven by recovery and residual-based error estimators for PHT-splines applied to time-harmonic acoustics