libname BIOS6102 "/home/u54386392/BIOS6102"; proc import datafile="/home/u54386392/BIOS6102/cvd_match1.csv" DBMS=CSV OUT=bios6102.cvd; getnames=yes; run; data cvd; set bios6102.cvd; label set='pair_ID' cvd='cardiovascular disease, 1:yes, 0:no' hyperten='hypertension'; run; /*transform data for McNemar's test*/ data cvd1 (keep=set cvd_hyperten); set cvd; if cvd=1; rename hyperten=cvd_hyperten; run; data cvd0 (keep=set nocvd_hyperten); set cvd; if cvd=0; rename hyperten=nocvd_hyperten; run;proc sort data=cvd1; by set; run; proc sort data=cvd0; by set; run; data new; merge cvd1 cvd0; by set; run; proc freq data=new;title "McNemar's test for Paired Samples";tables cvd_hyperten*nocvd_hyperten /agree expected norow nocol nopercent;run; proc logistic data=cvd;strata set;model cvd(event='1')=hyperten