I assume you know how to model the data to obtain the probabilities, and that you are asking about how to create the box plot of the results. You don't mention what version of SAS you are running, but here's how to create the plot in the current version: data Sim;
label p = "Survival Probability (%)";
call streaminit(1);
do i = 1 to 20;
Genotype = 1; Status = "Perishing";
p = rand("normal", 32, 10) / 100;
output;
Genotype = 1; Status = "Surviving";
p = rand("normal", 72, 5) / 100;
output;
Genotype = 2; Status = "Perishing";
p = rand("normal", 32, 10) / 100;
output;
Genotype = 2; Status = "Surviving";
p = rand("normal", 74, 8) / 100;
output;
end;
run;
proc sgplot data=Sim; format p PERCENT8.;
vbox p / category=Genotype group=Status nofill nooutliers groupdispay=overlay;
scatter y=p x=Genotype / group=Status;
run; One question you might have is when to use CATEGORY= and when to use GROUP=. For an explanation, see the article "What is the difference between categories and groups in PROC SGPLOT?"
... View more
