Bayesian Censored Model


/*--------------------------------------------------------------

                    SAS Sample Library

        Name: qliex08.sas
 Description: Example program from SAS/ETS User's Guide,
              The QLIM Procedure
       Title: Bayesian Censored Model
     Product: SAS/ETS Software
        Keys: Bayesian limited dependent variables
        PROC: QLIM
       Notes:

--------------------------------------------------------------*/


title1 'Bayesian Analysis';

ods graphics on;

data test;
   do i=1 to 200;
      e1 = rannor(8726)*2000;
      WinChance = ranuni(8772);
      Price = 10+ranexp(8773)*4;
      y = 48000 + 5000*WinChance - 100 * price +  e1;
      if y>50000 then TicketSales = 50000;
      if y<=50000 then TicketSales = y;
      output;
   end;
   keep WinChance price  y TicketSales;
run;


proc qlim data=test plots(prior)=all;
   model TicketSales = WinChance price;
   endogenous TicketSales ~ censored(lb=0 ub= 50000);
   prior intercept~normal(mean=48000);
   prior WinChance~normal(mean=5000);
   prior Price~normal(mean=-100);
   bayes NBI=10000 NMC=30000 THIN=1 ntrds=1 DIAG=ALL STATS=ALL seed=2;
run;