Weighted least squares response surface fitting in factorial paired comparisons