add function for computing probability scale fitted values
improve documentation/examples