%include "C:\Documents and Settings\pkraft\My Documents\webpage\softetc\wandel.sas";

/* Input data set contains variables zsnp1-zsnp10 that are counts of SNP minor alleles */
proc freq data=UseMe;
   tables d*(zsnp1 zsnp2); /* d is case-control indicator */
run;

%wandel(in=UseMe,vars=zsnp1 zsnp2 zsnp3 zsnp4 zsnp5 zsnp6 zsnp7 zsnp8 zsnp9 zsnp10,nrep=1000,k=3,caco=d);