%let DATA_DIR = C:\Users\baileraj\Desktop;
options nodate nocenter nonumber;

filename FROMWEB url "http://www.users.muohio.edu/baileraj/classes/sta402/Data/SMSA-DASL-2space-sep.txt ";

data SMSA_from_txt;
  infile FROMWEB firstobs=2;
* infile "C:\Users\baileraj\Desktop\SMSA-DASL-2space-sep.txt"
firstobs=2;
  length city $ 39;
  input city & JanTemp JulyTemp RelHum Rain Mortality
        Education PopDensity pct_NonWhite pct_WC pop
        pop_per_house income HCPot NOxPot S02Pot NOx;
run;

/* DEBUG option - print out N=60 observations in the SMSA data set */
/*
proc print data=SMSA_from_txt;
var city JanTemp JulyTemp RelHum NOxPot S02Pot NOx;
run;
*/

proc surveyselect data=SMSA_from_txt
                  method=srs n=6 out=sample_SMSA_cities;
run;

options fromdlim="-";
ods rtf bodytitle file="&DATA_DIR\Display-3.32.rtf";
proc print data=sample_SMSA_cities;
   var city JanTemp JulyTemp RelHum NOxPot S02Pot NOx;
run;
ods rtf close;



proc plan seed=776655;
  title "generate randomization/allocation";
  title2 "assign: 4 treatments to 20 EUs";
  factors EU=20;
  output out=design_crd;
run;

proc print data=design_crd;
run;

data design_crd2;
  set design_crd;
  ID = _n_;   * _N_ corresponds to observation # in dataset;
run;
proc print data=design_crd2 noobs;
run;

data design_crd3;
  set design_crd2;
  if ID <= 5 then TRT=1; * assign first 5 observations to trt 1;
  else if ID <= 10 then TRT=2; * assign next 5 obs. to trt 2;
  else if ID <= 15 then TRT=3; * assign next 5 obs. to trt 3;
  else TRT=4; * assign last 5 obs. to trt 4;
run;

proc print data=design_crd3;
run;

proc transpose data=design_crd3 out=t_design_crd prefix=C;
run; * transpose the data for printing;

ods rtf bodytitle file="&DATA_DIR\Display-3.34.rtf";
proc print data=t_design_crd;
run;
ods rtf close;


%let DATA_DIR=C:\Users\baileraj\BAILERAJ\Classes\Web-Classes\STA402\FALL-2007\data;
%let OUT_DIR=C:\Users\baileraj\Desktop;

filename FROMWEB url "http://www.users.muohio.edu/baileraj/classes/sta402/Data/ch2-dat.txt";

data nitrofen;
  infile FROMWEB expandtabs missover pad firstobs=16;
* reading from WEB vs. fixed folder;
* infile "&DATA_DIR\ch2-dat.txt" firstobs=16 expandtabs missover pad ;
input @9 animal 2. @17 conc 3. @25 brood1 2.
      @33 brood2 2. @41 brood3 2. @49 total 2.;
run;

ods rtf BODYTITLE file="&OUT_DIR\Display-4.1.rtf";;
proc tabulate data=nitrofen format=3.; * include default format for;
  class conc;   * classification variable;
  var brood1 brood2 brood3 total;  * numeric summaries;
  table (brood1 brood2 brood3 total)*conc,
                   n min q1 median*f=4.1 q3 max;
run;


proc tabulate data=nitrofen format=3.; * include default format for;
  class conc;   * classification variable;
  var brood1 brood2 brood3 total;  * numeric summaries;
  table conc*(brood1 brood2 brood3 total),
                   n min q1 median*f=4.1 q3 max;
run;


ods rtf close;


proc tabulate data=nitrofen;
  var brood1 brood2 brood3 total;
  table brood1 brood2 brood3 total, n min q1 median q3 max;
run;

options formdlim="-";
title;
proc tabulate data=nitrofen format=3.;
  class conc;
  var brood1 brood2 brood3 total;
  table (brood1 brood2 brood3 total)*conc,
         n min q1 median*f=4.1 q3 max;
run;

proc tabulate data=nitrofen format=3.;
  class conc;
  var brood1 brood2 brood3 total;
  table  brood1 brood2 brood3 total,conc,
         n min q1 median*f=4.1 q3 max;
run;