/* GENDATA.GPR: program to generate data for prisoner's dilemma problem

                John Rust, University of Maryland, November, 2005 */


   #include setup.gpr;  /* setup.gpr generates nobs x_a and x_b observations */

   p_a=zeros(nobs,1);
   p_b=zeros(nobs,1);
   y=zeros(nobs,2);
   seeds=rndu(nobs,2);

   i=1; do until i > nobs;

      {p_a[i],p_b[i]}=equil(x_a[i],x_b[i]);

      y[i,1]=(seeds[i,1] < p_a[i]);
      y[i,2]=(seeds[i,2] < p_b[i]);

   i=i+1; endo;

   "summary statistics";

   " Mean equilibrium probability of confessing:       "; meanc(p_a)~meanc(p_b);
   " Mean of a,b's actions (1=confess, 0 don't confess)"; meanc(y)';
   " Mean of a,b's observed types                      "; meanc(x_a)~meanc(x_b);

   save p_a,p_b,x_a,x_b,y,seeds;