Analysis code - Guangfu Hao (2025)