femalemort_AorB=[.6 .2]'; malemort_AorB=[.2 .2]'; femalemort_CorD=[.2 .2]'; malemort_CorD=[.2 .2]'; % C and D have no mortality differences alpha=1.2; dataA=NaN*ones([2000 6 20]); dataB=dataA; dataC=dataA; dataD=dataA; dataA=NaN*ones([2000 6 20]); dataB=dataA; dataC=dataA; dataD=dataA; done=zeros([4 20]); %save bigsimu