thselect.samsee function

Threshold selection via SAMSEE