Script for PSMC-like estimation