Towards ‘h-p adaptive’ generalized ANOVA