/* Lab 2 */

/* Getting Started */
filename kobe url 'http://www.openintro.org/stat/data/kobe.csv';

proc import
   datafile=kobe
   out=kobe(drop=time)
   dbms=csv
   replace; 
   getnames=yes;
   guessingrows=max;
run;

proc print data=kobe (obs=10);
run;

proc print data=kobe;
   where game=1 and quarter='1';
   var basket;
run; 

filename cstrk url "http://www.openintro.org/stat/data/calc_streak.sas";
%include cstrk;

%calc_streak(dset=kobe, streakvar=basket, outset=kobe_streak);

proc sgplot data=kobe_streak;
   title "Basket Streaks by Kobe Bryant in 2009 (133 attempts)";
   vbar streak;
run;


/* Simulations in SAS */
data coin (keep = outcome);
   flip = rand("Bernoulli",0.5);
   if flip = 1 then outcome="heads"; else outcome="tails";
run;

proc print data=coin;
run;

data sim_fair_coin(keep=outcome);
   do i = 1 to 100;
      flip = rand("Bernoulli",0.5);
	  if flip = 1 then outcome="heads"; else outcome="tails";
	  output;
   end;
run;

proc freq data=sim_fair_coin;
   tables outcome;
run;

/* Simulating the Independent Shooter */
data sim_basket(keep=basket);
   do i = 1 to 100;
      flip = rand("Bernoulli",0.5);
      if flip = 1 then basket="H"; else basket="M";
      output;
   end;
run;

data sim_basket(keep=basket);
   do i = 1 to 133;
      flip = rand("Bernoulli",0.45);
      if flip = 1 then basket="H"; else basket="M";
      output;
   end;
run;

data kobe_sim (keep=kobe sim);
	merge kobe(rename=(basket=kobe)) sim_basket(rename=(basket=sim));
run;

proc print data=kobe_sim;
run;