Extreme precipitation probability estimation is vital for hazard protection design but has high uncertainty. We tested six statistical models using 2000 years of climate data. Our Bayesian hierarchical duration-dependent Generalized Extreme Value model shows the highest accuracy and robustness for sample sizes between 30 and 100 years, making it highly promising for use with limited observational records.