Scripts and code for analyzing Psa screen