This article is dedicated to the design of a linear-circular regression technique and to its application to ground-based GNSS-Reflectometry (GNSS-R) altimetry. The altimetric estimation is based on the observation of the phase delay between a GNSS signal sensed directly and after a reflection off of the Earth's surface. This delay evolves linearly with the sine of the emitting satellite elevation, with a slope proportional to the height between the reflecting surface and the receiving antenna. However, GNSS-R phase delay observations are angular and affected by a noise assumed to follow the von Mises distribution. In order to estimate the phase delay slope, a linear-circular regression estimator is thus defined in the maximum likelihood sense. The proposed estimator is able to fuse phase observations obtained from several satellite signals. Moreover, unlike the usual unwrapping approach, the proposed estimator allows the sea-surface height to be estimated from datasets with large data gaps. The proposed regression technique and altimeter performances are studied theoretically, with further assessment on both synthetic and real data.