Wearable biosensors and smartphone applications can measure physiological variables over multiple days in free-living conditions, revealing circadian rhythms and responses to external stressors such as meals and physical activity. Here we develop a probabilistic Bayesian framework to learn interpretable, personal parameters from wearable time series data. We measure food and drink ingestion, glucose dynamics, physical activity, heart rate (HR) and heart rate variability (HRV) in 25 healthy participants over 14 days. Modelling ingestion events with glucose reveals that slow glucose decreases are associated with large postprandial glucose spikes, and we uncover a circadian baseline rhythm in glucose levels with high amplitudes in some individuals. Physical activity and circadian rhythms explain 40- 65% of HR variance, whereas the variance explained for HRV is more heterogeneous across individuals (20-80%). Finally, incorporating activity, HR and HRV in the modelled glucose explains an additional 10% glucose variability in some individuals, highlighting the relevance of integrating multiple physiological signals for a complete and predictive understanding of glucose dynamics.