1. INTRODUCTION

2. DATA

3. METHODOLOGY
3.1. The Models
3.2. Continuous Fits
3.3. Implementation
3.4. Model Fitting

4. THE BIVARIATE HI MASS FUNCTION FOR MASSIVE GALAXIES
4.1. Model Comparison
4.2. Covariance of Parameters
4.3. Model fit at high HI masses and comparison with ALFALFA
4.4. Model Selection and Validation

5. DISCUSSION
5.1. The Shape of the HIMF
5.2. The Dependence of the HI Mass Function on Stellar Mass and Star Formation Rate
5.3. Physical Drivers of the HIMF

6. SUMMARY AND CONCLUSIONS