*Generate the parametric bootstrap data;
DATA ParBootstraps (KEEP=Replicate waiting);
  ARRAY MLE	{5} p mu sigma nu tau (0.361 54.615 5.871 80.091 5.868);
  DO Replicate=1 TO 10000;
    DO i=1 TO 272;
      b=RANBIN(0,1,p);
      waiting=b*(mu+sigma*RANNOR(0))+(1-b)*(nu+tau*RANNOR(0));
      OUTPUT;
    END;
  END;
RUN;

ODS OUTPUT ParameterEstimates=PBootPars;
PROC NLMIXED DATA=ParBootstraps;
  BY Replicate;
  PARMS p=0.5 mu=55 sigma=5 nu=80 tau=5;
  BOUNDS 0<p<1, 0<sigma, 0<tau;
  ll=log( p*PDF("NORMAL",waiting,mu,sigma)+
          (1-p)*PDF("NORMAL",waiting,nu,tau) );
  MODEL waiting ~ GENERAL(ll);
RUN;

PROC UNIVARIATE DATA=PBootPars NOPRINT;
  CLASS parameter (ORDER=DATA);
  VAR estimate;
  OUTPUT OUT=BootStrapCIs STD=Std_Error PCTLPTS=2.5 97.5 PCTLPRE=Q; 
RUN;

PROC PRINT DATA=BootStrapCIs NOOBS LABEL; 
  LABEL Q2_5="CI lower limit" Q97_5="CI upper limit";
RUN;