SAS code replicates Dinh2013's logistic regression