This study integrates longitudinal systolic and diastolic cardiac function measures using Bayesian modeling to identify trajectories associated with distinct heart failure subtypes and proteomic signatures, enhancing risk prediction and therapeutic target discovery.