Data transformation and model selection in bivariate allometry

Students of biological allometry have used the logarithmic transformation for over a century to linearize bivariate distributions that are curvilinear on the arithmetic scale. When the distribution is linear, the equation for a straight line fitted to the distribution can be back-transformed to form a two-parameter power function for describing the original observations. However, many of the data in contemporary studies of allometry fail to meet the requirement for log-linearity, thereby precluding the use of the aforementioned protocol. Even when data are linear in logarithmic form, the two-parameter power equation estimated by back-transformation may yield a misleading or erroneous perception of pattern in the original distribution. A better approach to bivariate allometry would be to forego transformation altogether and to fit multiple models to untransformed observations by nonlinear regression, thereby creating a pool of candidate models with different functional form and different assumptions regarding random error. The best model in the pool of candidate models could then be identified by a selection procedure based on maximum likelihood. Two examples are presented to illustrate the power and versatility of newer methods for studying allometric variation. It always is better to examine the original data when it is possible to do so.