SAS code from stat 135