I think you may be looking for something like this:
if male=1 then death=Rand('bernoulli',0.6);
else if male=0 then death=Rand('bernoulli',0.7);
The bernoulli is a single observation, valued 1/0 where the probability of getting a 1 is specified as the parameter P of the function.
You can test this with this code. The Mean result of a 1/0 coded variable is the percent of 1's.
data example;
do male=1,0;
do trial=1 to 10000;
if male=1 then death=Rand('bernoulli',0.6);
else if male=0 then death=Rand('bernoulli',0.7);
output;
end;
end;
run;
proc means data=example;
class male;
var death;
run;
Note that if your data doesn't do enough simulations the random nature means the result may not be really close to .7 or .6. Change the Trial=1 to 10; above and you can get result.
... View more
