data canada; infile 'Macintosh HD:work:626:sex_can.dat'; /* change to appropriate path */ input year fc fnf fpei fns fnb fq fo fm fs fa fbc mc mnf mpei mns mnb mq mo mm ms ma mbc ; mfc = mc+fc; ratio=mc/fc; proc means n min max mean sum ; var mc fc mfc; proc genmod ; model mc/mfc = / link = identity dist = binomial; proc genmod ; model mc/mfc = / link = logit dist = binomial; proc genmod ; model mc/mfc = / link = log dist = binomial; run;