/* Monte Carlo study of two independent samples from two Normal populations */ 
/* Macro Variables used to initialize the simulation */

libname midterm 'c:\Stat6250mt';

* Group 1;
%let n1 = 78;
%let mean1 = 56;
%let stdev1 = 3;
%let dist1 ='NORMAL';

* Group 2;
%let n2 = 147;
%let mean2 = 72;
%let stdev2 = 2.5;
%let dist2 = 'NORMAL';

data midterm.data01;
	do i=1 to &n1; drop i;
		y = rand(&dist1, &mean1, &stdev1);
		samp = 'A';
		output midterm.data01;
	end;
run;

data midterm.data02;
	do i=1 to &n2; drop i;
		y = rand(&dist2, &mean2, &stdev2);
		samp = 'B';
		output midterm.data02;
	end;
run;

* merge the data into one dataset with the same variables;
data midterm.alldata02;
	set midterm.data01 midterm.data02;
run;

proc print data=midterm.alldata02;
run;

proc univariate data=midterm.alldata02;
	by samp;
run;

proc boxplot data=midterm.alldata02;
	plot y*samp / boxstyle=schematic boxwidth=15;
run;

proc ttest data=midterm.alldata02;
	class samp;
	var y;
run;