Multivariate analysis of globular cluster horizontal branch morphology: searching for the second parameter