calculatePhat function

calculatePhat