Maximum likelihood method for fitting the Fundamental Plane of the 6dF Galaxy Survey

We have used over 10,000 early-type galaxies from the 6dF Galaxy Survey (6dFGS) to construct the Fundamental Plane across the optical and near-infrared passbands. We demonstrate that a maximum likelihood fit to a multivariate Gaussian model for the distribution of galaxies in size, surface brightness and velocity dispersion can properly account for selection effects, censoring and observational errors, leading to precise and unbiased parameters for the Fundamental Plane and its intrinsic scatter. This method allows an accurate and robust determination of the dependencies of the Fundamental Plane on variations in the stellar populations and environment of early-type galaxies.