In reply to Wang Yajun.

Thanks for your comment. In the second paragraph there is a link to download to SAS code.