/* Routine to draw Venn diagram for combinations of outcome variables    */
/* November 18 1998  J.M.Lauritsen                          STB-47 gr34            
see help venndiag
 globals returned in _results:
   S_1    venndiagram fail/ok
   S_2     N records in presented datafile
   S_3     N records after exluding [in/if] clause
   S_4     Sum of number of variables with missing
   S_5     N records after excluding missing
   S_6     N records in diagram
   S_7     condition included in Label()
   S_8 .. S_11 Value to count as yes in v1..v4
   S_12 .. S_15    variable v1...v4 (name)
   S_20 .. S_35    counts from areas
   S_36..S_39      counts of variables v1..v4  
   all reset from start
*/
*! Venn Diagram graphs Version 1.00, Jens M.Lauritsen 18.nov 1998

program define venndiag
  version 5.0
    #delimit ;
    local varlist "required existing min(2) max(4)";
    local options "C1(integer 1) C2(integer 1) C3(integer 1) C4(integer 1)
    LIst(string) T1title(string) T2title(string) T3title(string) LAbel(string) MISsing 
    NOFrame NOLabel SAving(string) PRint SHow(string) Gen(string) NOGraph
    R1title(string) R2title(string) R3title(string) 
    R4title(string) R5title(string) R6title(string)*" ;
  #delimit cr 

  local if "opt"
  local in "opt"
  parse "`*'"
  preserve
  if "`label'" == "all" {local label = "acimdft" }
  if "`label'" == "" {local label = "adft" }

  global S_7 = "`label'"
  if "`show'" == "" {local show = "cptf" }
  if "`show'" == "all" {local show = "cvtpfn" }

  if "`t1title'" == "" & index("`label'","x") == 0 {local t1title = "Venn Diagram"}
  if "`gen'" == "" {local gen = "vdx1129"}

confirm new var `gen'
confirm new var vd_id299
confirm new var _merge

local file = "File: " + "$S_FN" + " (Cr:" + substr("$S_FNDATE",1,12) + ")"

local i = 1
while `i' <= 35 {
   global S_`j'
   local i = `i' + 1
   }

* Count number of variables
  parse "`varlist'", parse(" ")
  local varno =0
  while "`1'"~=""  {
     local varno=`varno'+1
     local j = `varno'+11
     global S_`j' = "`1'"
    macro shift
   }

*error if inapropriate arguments:
 vd_exit `varno' "`varlist'"
 if "$S_1" == "fail" {   exit}

*define count value for S_ globals

local i = 1
while `i' <= `varno' {
   local j = `i'+7
   global  S_`j' = `c`i''
   local i = `i' + 1
   }
  
  qui count
  global S_2 = _result(1)
  
  * apply if and in
  capture keep `if'
  capture keep `in'
  qui count
  local before = _result(1)
  global S_3 = _result(1)
  
* save memory and gain speed
  
  keep `varlist' `list'
  

tempfile all miss
qui gen vd_id299 = _n
sort vd_id299
qui gen str6 `gen' = ""
* now we have a pseudo id variable vd_id299 in the datasets

capture save `all'

*count number of records with missing included in the diagram (when missing option present):
  local nmiss = 0   
  if "`missing'" != "" {
      parse "`varlist'", parse(" ")
      while "`1'"~=""  {
        quietly count if `1' == .
        local nmiss = `nmiss' + _result(1)
        macro shift
      }
   }
  global S_4 = `nmiss'

* drop all missing
   parse "`varlist'", parse(" ")
      while "`1'"~=""  {
        quietly drop if `1' == .
        macro shift
      }

capture save `miss'

if "`missing'" != "" {capture use `all',clear}
  qui count 
  global S_5 = $S_3 - _result(1)


*present documentation:
   vd_doc "`varlist'" `varno' "`if'"  "`in'"  "`missing'" "`file'" "`gen'" 
tempfile ok
capture save `ok'

* draw graphs 
if "`nograph'" == "" {
  if "`t2title'" == "" & index("`label'","x") == 0 {local t2title = "N = " + "$S_6"}
  if "`saving'" == "" {gph open}
  if "`saving'" != "" {gph open,saving("`saving'",replace)}
*  if "$S_grid" == "y" {vd_grid}   /* for testing and lookup: Global S_grid = "y" */  

  gph pen 1

  * draw boxes
    vd_box `varno' `noframe'

   
  if index("`label'","f") > 0 {
   vd_font 0.48 
   gph text 22000 31500 0 1 `file'
   }

  * add labels:
     vd_lbl `varno' "`show'" "`t1title'" "`t2title'" "`t3title'"  /*
             */ "`nolabel'" "`varlist"' "`all'" "`miss'"


 if index("`label'","x") > 0 {     /*show user specified information */
    local line = 6000
    local left = 18500
    local height = 1000
    vd_font 0.82 
    gph text `line' `left' 0 -1  `r1title' 
    vd_font 0.65
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `r2title' 
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `r3title'
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `r4title' 
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `r5title' 
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `r6title'
  }

  * add variables in boxes:
  vd_count "`show'" `varno'

  vd_font 1.0
  gph close
  * check for Mach and Win95
  if ("$S_MACH" == "PC" & "$S_OSDTL" == "95"){
    if "`saving'" != "" {gphprint ,nologo saving("`saving'.wmf",replace)}
    if "`print'" != "" {gphprint , nologo} 
  }
}

capture use `ok', clear
if "`gen'" != "vdx1129" { /* add variable to original file */
  drop `varlist' 
  qui sort vd_id299
  capture save `all',replace
  restore
  capture drop vd_id299
  capture drop `gen' 
  gen vd_id299 = _n
  qui  sort vd_id299 
  merge vd_id299 using `all'
  qui  replace `gen' = "miss" if _merge == 1
  di _newline
  drop vd_id299 _merge
  local vdlbl = "`varlist'" + "(vd)"
  label var `gen' "`vdlbl'"  
  local t = "$S_12" + ":"+ "`c1'" + " "+ "$S_13:" + "`c2'" + " "
  if `varno' > 2 {local t = "`t'" + " " +"$S_14" + ":"+ "`c3'" + "  " }
  if `varno' > 3 {local t = "`t'"+ " " +"$S_15" + ":"+ "`c4'"}
  di  "New variable created. Name: " in green "`gen'" in yel "  Label: " in green  "`vdlbl'"
  di  "notes added: "
  note `gen': generated by venndiag.ado on TS  . Variables and values were:
  note `gen': `t'
  local t = "`if'"+" "+ "`in'"
  if length("`t'") > 0 {
    local t = "`gen'"+"=miss for "+ string(($S_2-$S_3)) + " records excluded by:"+ "`t'"
    note `gen': `t'
    }

  if $S_5 > 0 {
    local t = "`gen'"+"=miss for " + string($S_5) + " records with missing values (.)"
    note `gen': `t'
    }

  if length("`missing'") > 0 & $S_4 > 0 {
      local t = "note that some records ("+ string($S_4) + " variables) had missing values (.)"
    note `gen': "`t'"
    }

  capture compress `gen'
  note `gen'
  }
  else { 
  drop vdx1129
  local gen = "" 
}

if "`list'" != "" {list `varlist' `list' `gen', nodisplay}
exit

end

* ------------------------------------------------------------------------------
* Some subroutines for various purposes

program define vd_lbl
   local varno = `1'
   local label = "$S_7"
   local show  = "`2'"
   local t1title = "`3'"
   local t2title = "`4'"
   local t3title = "`5'"
   local nolabel = "`6'"
   local varlist = "`7'"
   local all = "`8'"
   local miss = "`9'"

   local top   = 6000
   local left  = 18500
   local height = 750
   local line = `top' - 4000
 
  if "`nolabel'" != "" {exit}

  if index("`label'","a") > 0 {
     local fn = substr("$S_DATE",1,12) 
     vd_font 0.48
     gph text 21000 31500 0 1 `fn'
     }

  if index("`label'","t") > 0 {     /*show information of records */
    local line = `top' - 3000
    vd_font 0.82 
    local l1 = cond(`varno' == 2,3000,1500)
*    if `varno' == 4  {local l1 = 1000}
    gph text 1000 `l1' 0 -1  `t1title'
    vd_font 0.65
    gph text `line' `left' 0 -1  `t2title' 
    local line = `line' + int(1.5*`height') 
    gph text `line' `left' 0 -1  `t3title'

  }


 if `varno' < 4 {local line = `line' + 0.5*`height'}

 if index("`label'","d") > 0 {           /*show description*/
    local line = `line' + int(1.5*`height')
    local i = 1
    parse "`varlist'", parse(" ")
    while "`1'"~=""  {
      local varlbl : variable label `1'   
      if "`varlbl'"=="" {local varlbl="`1'"}  
      local t = substr("ABCD",`i',1) + " `varlbl'"
      gph text `line' `left' 0 -1 `t'
       local i = `i' + 1
       local line = `line' + 1.25*`height'     
       mac shift
     }    
  local line = `line' + 1.25*`height'
  }  

  if `varno' < 4 {local line = `line' + 0.5*`height'}

  if index("`label'","c") > 0 {     /*show counting */
    gph text `line' `left' 0 -1 Value indicators:
    local line = `line' + int(1.25*`height')
    local i = 1
    parse "`varlist'", parse(" ")
    while "`1'"~=""  {
        local t = ""
        if `i' == 1  {local t = "A: (" + "`1'" + "=" + string($S_8)+")" }       
        if `i' == 2  {local t = "B: (" + "`1'" + "=" + string($S_9)+")"}        
        if `i' == 3  {local t = "C: (" + "`1'" + "=" + string($S_10)+")"}        
        if `i' == 4  {local t = "D: (" + "`1'" + "=" + string($S_11)+")"}        
       gph text `line' `left' 0 -1 `t'
       local i = `i' + 1
       local line = `line' + 1.25*`height'     
       mac shift
     }
  local line = `line' + `height'
  }

  if `varno' < 4 {local line = `line' + 0.5*`height'}

  if index("`label'","i") > 0 {     /*show information of records */
     local t = " Records in file: " 
     gph text `line' `left' 0 -1 `t'
     local t = string($S_2)
     local j = `left' + 12000
     gph text `line' `j' 0 1 `t'
     local t = " Excluded: Miss " + string($S_5) + " In/if: " + string($S_2-$S_3)
     vd_font 0.5
     if length("`t'") > 35 {vd_font 0.4}
     local line = `line' + `height'
     gph text `line' `left' 0 -1 `t'
     vd_font 0.65 
     local t = string($S_5+$S_2-$S_3)
     gph text `line' `j' 0 1 `t'
     local line = `line' + 1.1*`height'
     gph text `line' `left' 0 -1  Total Records in graph:
     local t = string($S_6)
     gph text `line' `j' 0 1 `t'
   }

  local line = `line' + 1.4*`height'
  if `varno' < 4 {local line = `line' + 0.5*`height'}

  if index("`label'","m") > 0 {/*show missing in variables */
    capture use `all',clear
    vd_font 0.48
    gph text `line' `left' 0 -1 Missing values (Records)
    vd_font 0.45 
    local l1 =  "(EXcluded from graph)"
    if $S_5 == 0  {local l1 =  "(Included in graph)"}
    gph text `line' 31000 0 1 `l1'
    vd_font 0.48
    local line = `line' + int(1.20*`height')
    local left = `left' + 500
    local i = 1
    parse "`varlist'", parse(" ")
    while "`1'"~=""  {
        local t = ""
        quietly count if `1' == .
        local t = substr("ABCD",`i',1) + ":(" + "`1'" + "=" + string(_result(1))+")"
*        if `i' == 4  {local t = "D: (" + "`1'" + "=" + string(_result(1))+")"}        
       gph text `line' `left' 0 -1 `t'
       local i = `i' + 1
       local line = `line' + `height'    
       mac shift
     }
   capture use `ok', clear
*    if "`missing'" != "" {capture use `miss',clear}     /*use missing in variables */
  }
  
 vd_font 0.65
end

program define vd_pct
      if `3' == 0 & `4' != 1 {exit}
      local l1 = round(100*(`3'/$S_6),1)
      local t = "(" + "`l1'" + " %)"
      vd_font 0.50
      local y = `1' + 850
      gph text `y' `2' 0 0 `t'
      vd_font 0.65
end

program define vd_text
  * place text on graph
  if `4' == 0 & `5' != 1 {exit}
  local t = string(`4')
  gph text `1' `2' 0 `3' `t'
end

program define vd_font
   * multiply size
     local vds: set textsize   
     local vds = `vds'/100
       local x = int(`vds'*923*(`1'))
       local y = int(`vds'*444*(`1'))
     gph font  `x' `y'
end

program define vd_box
* draw boxes on graph
  local varno `1' 
  local frame `2'
  local leftm = 2500

  if `varno' < 4 {
      vd_gphsq 2000 (`leftm'+3500) 12000 6500 
      vd_gphsq 8000 (`leftm'+7000) 12000 7000 
      if `varno' == 3 {
          vd_gphsq 11000 `leftm' 6500 11000}
   }

  if "`frame'" == "" {
   if `varno' == 2 {vd_gphsq 1600 (`leftm'+1000) (`leftm'+17000) 13500}
   if `varno' == 3 {vd_gphsq 1600 (`leftm'-1000) (`leftm'+17000) 15500}
   if `varno' == 4 {vd_gphsq 1600 (`leftm'-1500) (`leftm'+17000) 17250}
   }

 if `varno' == 4 {
  local leftm = `leftm' - 1500
  vd_gphsq 2000 (`leftm'+5900) 17250 6800  /* a*/
  vd_gphsq 6500 (`leftm'+9300) 14250 7300 
  vd_gphsq 9500 (`leftm'+300) 7000 15900   /*c*/
  vd_gphsq 6250 (`leftm'+700) 6500 16250 /* d */
   }
end

*draw grid on screen
program define vd_grid
local i = 5050
gph pen 4
while `i' < 20000 {
*   local txt = string(`i')
   local txt = substr(string(`i'),1,2)
   gph text `i' 500 0 -1 `txt'
   gph line  `i' 500  `i' 31500
   local i = `i' + 2500
   }
local i = 7550
while `i' < 25000 {
   local txt = substr(string(`i'),1,2)
   gph text  500 `i' 0 -1 `txt'
   gph line  500 `i' 22000 `i'
   local i = `i' + 2500
   }
gph pen 1
end


* draw square on graph 
program define vd_gphsq
  local top = `1'      /* top row  */
  local left = `2'     /* top column */
  local bottom = `1' + `3'     /* rt + height   */
  local right  = `2' + `4'     /* ct+ width     */
  gph line `top' `left' `bottom'  `left'
  gph line `top' `left' `top'  `right'
  gph line `top' `right' `bottom'  `right'
  gph line `bottom' `left' `bottom'  `right'
end

program define vd_exit
local varno   = `1'
local varlist = "`2'"

  parse "`varlist'", parse(" ")
  local string = 0
  local vrs = ""
   while "`1'"~=""  {
     local fmt : format `1'
     if index("`fmt'","s") > 0 {local string = `string' + 1}
     if index("`vrs'","`1'") > 0 {local string = `string' + 1}    /*repeating var*/
     local vrs = "`vrs'" + "`1'"
    macro shift
   }

*error if inapropriate arguments:

 if `string' > 0 {
     di in red _dup(70) "_" _newline
     disp in red "Errors in command, therefore quitting" _newline         
     disp _newline in red "Variables must be numeric and different:" _newline
     d 
     di in red _dup(70) "_" _newline _newline _newline
     global S_1 = "fail"
 }
end

program define vd_doc
* display basis for diagram in results window
   local varlist = "`1'"
   local varno   = `2'  
   local if      = "`3'"
   local in      = "`4'"
   local missing = "`5'"
   local file    = "`6'"
   local vdvar   = "`7'"

   qui count
   global S_6 = _result(1)
   local after = _result(1)
   local before = $S_3
   local nmiss =  $S_4
   local total = $S_2
   
   di in green _dup(70) "_" _newline
   
   di in yellow "Venn diagram of variables: " in green "`varlist'"
   di in yellow "`file'" _newline
  * and variables: 
  local i = 1
  di in yellow "     Outcome    Variable and label "
  parse "`varlist'", parse(" ")
  local varno =0
  while "`1'"~=""  {
     local varno=`varno'+1
     local var`varno' = "`1'"
     local varlbl`varno' : variable label `1'
     local j = `i' + 7
     if "`varlbl`varno''"=="" {local varlbl`varno'="`1'"}
     * can't get macro substition to work - saving time
     if `i' == 1  {di in yellow substr("ABCD",`i',1) ": " _col(10) $S_8  _col(22) "`var`i''"  "    "  _column(22) "`varlbl`i''"}
     if `i' == 2  {di in yellow substr("ABCD",`i',1) ": " _col(10) $S_9  _col(22) "`var`i''"  "    "  _column(22) "`varlbl`i''"}
     if `i' == 3  {di in yellow substr("ABCD",`i',1) ": " _col(10) $S_10 _col(22) "`var`i''"  "    "  _column(22) "`varlbl`i''"}
     if `i' == 4  {di in yellow substr("ABCD",`i',1) ": " _col(10) $S_11 _col(22) "`var`i''"  "    "  _column(22) "`varlbl`i''"}
     local i = `i' + 1
     mac shift
     }
 
   local siz = length(string(`total'))+8
   di _newline in green _column(5) `total'  _column(`siz') "Records in file"
   
   if "`if'" != "" | "`in'" != "" {
          local col = length(string(`total'))+5 -length(string(`total'-`before'))
          di in green _column(`col') `total'-`before' in yellow _column(`siz') "Records excluded by [in/if] clauses"
          }

   if "`missing'" == "" {
         local col = length(string(`total'))+5 -length(string(`before'-`after'))
         di  in green _column(`col') `before'-`after' in yellow _column(`siz') "Records excluded by missing values"
         }

   local col = length(string(`total'))
   di in green _column(5) _dup(`col') "_"

   local col = length(string(`total'))+5 -length(string(`after'))
   di in green _column(`col') `after' in yellow _column(`siz') "Records in Diagram: "   

   if "`missing'" != ""  {
        local col = length(string(`total'))+5 -length(string(`nmiss'))
        di _newline in red _column(`col') `nmiss' in yellow _column(`siz') "variables in all records contain missing values"
        parse "`varlist'", parse(" ")
        while "`1'"~=""  {
          quietly count if `1' == .
          local col = `siz' + 13 
          di in yellow _column(`siz') "`1'" _column(`col') ":" _continue
          local col = 8 - length(string(_result(1)))
          di in red _column(`col') _result(1) 
          macro shift
       }
   }  
  vd_tabl "`vdvar'" `varno'

  di in green _dup(70) "_"
end

program define vd_tpct
      local l1 = round(100*(`1'/$S_6),1)
      local t = "`l1'" + " %"
      vd_font 0.60
      gph text `2' `3' 0 0 `t'
      vd_font 0.65
end

program define vd_count
local show = "`1'"
local varno = `2'

local zero = 1

if index("`show'","f") > 0 & (index("`show'","x") == 0 | index("`show'","t") > 0) {vd_font 0.50
         gph text 22000 10000 0 1 % of total
         }

 if index("`show'","f") > 0 & index("`show'","p") > 0 {vd_font 0.50
         gph text 22000 10400 0 -1 (% in area of total)
         }

vd_font 0.65

if index("`show'","n") == 0 {local zero = 0}

if `varno' < 4 {
 
   if index("`show'","v") == 0 {
          vd_font 0.70
          local t = "A"}
       else {local t = "$S_12"}       
       gph text 3300 7000 0 -1 `t'
   if index("`show'","t") > 0 {vd_tpct $S_36 4200 8000}
   if index("`show'","v") == 0 {
          vd_font 0.70
          local t = "B"}
       else {local t = "$S_13"}
       gph text 9350 12700 0 -1 `t'
    if index("`show'","t") > 0 {vd_tpct $S_37 10350 14000}
    }


if `varno' == 3 {                                 /* add C */
    if index("`show'","v") == 0 {
          vd_font 0.70
          local t = "C"}
      else { local t = "$S_14"}
       gph text 12500 2750 0 -1 `t'
    if index("`show'","t") > 0 {vd_tpct $S_38 13500 4250}
  }


 if `varno' == 4 {               
                                  /* add letters */
   if index("`show'","v") == 0 { 
          vd_font 0.70
        gph text 3300 7500 0 -1 A
        gph text 20250 11000 0 -1 B
        gph text 14500 1500 0 -1 C
        gph text 7500 2000 0 -1 D
     }
     else {vd_font 0.65
       gph text 3300 7500 0 -1 $S_12
       gph text 20250 10500 0 -1 $S_13
       gph text 14500 1500 0 -1 $S_14
       gph text 7500 2000 0 -1 $S_15
     }
   vd_font 0.65
   if index("`show'","t") > 0 {
      vd_tpct $S_36 4200 8500
      vd_tpct $S_37 20250 15000
      vd_tpct $S_38 15500 3000
      vd_tpct $S_39 8500 3000
   }
  }

/* add counts and percentages */


if `varno' == 2 {
  if index("`show'","c") {vd_text 10000 11000 0 $S_25  `zero'}     /* ab */
  if index("`show'","p") > 0 {vd_pct 10000 11000 $S_25 `zero'}

  if index("`show'","c") {vd_text 16500 11000 0 $S_22 `zero'}     /*    b */
  if index("`show'","p") > 0 {vd_pct 16500 11000 $S_22 `zero'}

  if index("`show'","c") {vd_text 6000 8500 0 $S_21 `zero'}       /* a */ 
  if index("`show'","p") > 0 {vd_pct 6000 8500 $S_21 `zero'}

  if index("`show'","x") == 0 & index("`show'","c") {vd_text 4000 14500 0 $S_20 `zero'}
  if index("`show'","x") == 0 & (index("`show'","t") > 0 | index("`show'","p")){vd_tpct $S_20 5000 14500 `zero'}
}

if `varno' == 3 {
  if index("`show'","c") {vd_text 12500 11000 0 $S_31  `zero'}   /* a b c */
  if index("`show'","p") > 0 {vd_pct 12500 11000 $S_31 `zero'}

  if index("`show'","c") {vd_text 6000 8000 0 $S_21  `zero'}  /* a      */
  if index("`show'","p") > 0 {vd_pct 6000 8000 $S_21 `zero'}

  if index("`show'","c") {vd_text 18500 13500 0 $S_22  `zero'}  /*    b   */
  if index("`show'","p") > 0 {vd_pct 18500 13500 $S_22 `zero'}

  if index("`show'","c") {vd_text 15500 4500 0 $S_23  `zero'}  /*     c */
  if index("`show'","p") > 0 {vd_pct 15500 4500 $S_23 `zero'}

  if index("`show'","c") {vd_text 9250 11000 0 $S_25  `zero'}  /* a b  */
  if index("`show'","p") > 0 {vd_pct 9250 11000 $S_25 `zero'}

  if index("`show'","c") {vd_text 12500 7500 0 $S_26  `zero'}  /* a   c*/
  if index("`show'","p") > 0 {vd_pct 12500 7500 $S_26 `zero'}
 
  if index("`show'","c") {vd_text 15500 11000 0 $S_28  `zero'} /*   b c*/
  if index("`show'","p") > 0 {vd_pct 15500 11000 $S_28 `zero'}

  if index("`show'","x") == 0 & index("`show'","c") {vd_text 4000 14500 0 $S_20  `zero'}  /*       */ 
  if index("`show'","x") == 0 & (index("`show'","t") > 0 | index("`show'","p")) {vd_tpct $S_20 5000 14500 `zero'}
}

if `varno' == 4 {
  if index("`show'","c") {vd_text 4000 12000 0 $S_21  `zero'}  /* a      */
  if index("`show'","p") > 0 {vd_pct 4000 12000 $S_21 `zero'}

  if index("`show'","c") {vd_text 17500 15500 0 $S_22  `zero'}  /*    b   */
  if index("`show'","p") > 0 {vd_pct 17500 15500 $S_22 `zero'}

  if index("`show'","c") {vd_text 14500 5500 0 $S_23  `zero'}  /*     c */
  if index("`show'","p") > 0 {vd_pct 14500 5500 $S_23 `zero'}

  if index("`show'","c") {vd_text 7500 5500 0 $S_24  `zero'}  /*     d */
  if index("`show'","p") > 0 {vd_pct 7500 5500 $S_24 `zero'}

  if index("`show'","c") {vd_text 17500 11500 0 $S_25  `zero'}  /* ab*/
  if index("`show'","p") > 0 {vd_pct 17500 11500 $S_25 `zero'}
  if "$S_grid" == "y" {vd_grid}   /* for testing and lookup: Global S_grid = "y" */  

  if index("`show'","c") {vd_text 14500 8500 0 $S_26  `zero'}  /* ac*/
  if index("`show'","p") > 0 {vd_pct 14500 8500 $S_26 `zero'}

  if index("`show'","c") {vd_text 7500 8500 0 $S_27  `zero'}  /* ad */
  if index("`show'","p") > 0 {vd_pct 7500 8500 $S_27 `zero'}

  if index("`show'","c") {vd_text 14500 15500 0 $S_28  `zero'} /* bc*/
  if index("`show'","p") > 0 {vd_pct 14550 15500 $S_28 `zero'}

  if index("`show'","c") {vd_text 7500 15500 0 $S_29  `zero'} /* bd*/
  if index("`show'","p") > 0 {vd_pct 7500 15500 $S_29 `zero'}

  if index("`show'","c") {vd_text 11000 5500 0 $S_30 `zero'}  /* cd*/
  if index("`show'","p") > 0 {vd_pct 11000 5500 $S_30 `zero'}

  if index("`show'","p") > 0 {vd_pct 14500 12000 $S_31 `zero'}
  if index("`show'","c") {vd_text 14500 12000 0 $S_31  `zero'} /*abc*/

  if index("`show'","p") > 0 {vd_pct 7500 12000 $S_32 `zero'}
  if index("`show'","c") {vd_text 7500 12000 0 $S_32 `zero'}   /* abd */

  if index("`show'","p") > 0 {vd_pct 11000 8500 $S_33 `zero'}
  if index("`show'","c") {vd_text 11000 8500 0 $S_33 `zero' }   /* acd */

  if index("`show'","p") > 0 {vd_pct 11000 15500 $S_34 `zero'}
  if index("`show'","c") {vd_text 11000 15500 0 $S_34 `zero' }   /* bcd */

  if index("`show'","c") {vd_text 11000 12000 0 $S_35  `zero'}   /* abcd*/
  if index("`show'","p") > 0 {vd_pct 11000 12000 $S_35 `zero'}

  if index("`show'","x") == 0 & index("`show'","c") {vd_text 4000 15500 0 $S_20  `zero'}  /*       */ 
  if index("`show'","x") == 0 & (index("`show'","t") > 0 | index("`show'","p")) {vd_tpct $S_20 5000 15500 `zero'}
}

/* end counts and percentages */

end

program define vd_tabl
 local vdvar = "`1'"
 local varno = `2'


qui count if ($S_12 == $S_8)
  global S_36 = _result(1)

  qui count if ($S_13 == $S_9)
  global S_37 = _result(1)

if `varno' > 2 {
 qui count if ($S_14 == $S_10)
 global S_38 = _result(1)
  if `varno' > 3 {
    qui count if ($S_15 == $S_11)
    global S_39 = _result(1)
  }
 }

qui replace `vdvar' = "error"

set d l 132
di "Counts for combined variables:"
di _dup(70) "-"

if `varno' == 2 {
  vd_val "A"  2 `vdvar' 21
  vd_val "B"  2 `vdvar' 22
  vd_val "AB" 2 `vdvar' 25
  vd_val "--"  2 `vdvar' 20
}

if `varno' == 3 {
  vd_val "A"  3 `vdvar' 21
  vd_val "B"  3 `vdvar' 22
  vd_val "C" 3 `vdvar'  23
  vd_val "AB" 3 `vdvar' 25
  vd_val "AC" 3 `vdvar' 26
  vd_val "BC" 3 `vdvar' 28
  vd_val "ABC" 3 `vdvar' 31 
  vd_val "---"  3 `vdvar' 20
  }

if `varno' == 4 {
  vd_val "A"  4 `vdvar' 21
  vd_val "B"  4 `vdvar' 22
  vd_val "C" 4 `vdvar'  23
  vd_val "D" 4 `vdvar'  24
  vd_val "AB" 4 `vdvar' 25
  vd_val "AC" 4 `vdvar' 26
  vd_val "AD" 4 `vdvar' 27
  vd_val "BC" 4 `vdvar' 28
  vd_val "BD" 4 `vdvar' 29
  vd_val "CD" 4 `vdvar' 30
  vd_val "ABC" 4 `vdvar' 31
  vd_val "ABD" 4 `vdvar' 32
  vd_val "ACD" 4 `vdvar' 33
  vd_val "BCD" 4 `vdvar' 34
  vd_val "ABCD" 4 `vdvar' 35
  vd_val "----"  4 `vdvar' 20
}
di _dup(70) "-"

exit
end

program define vd_val
 local gen = "`3'"
 local a = "($S_12 == $S_8)"
 local b = "&($S_13 == $S_9)"
 local c = "&($S_14 == $S_10)"
 local d = "&($S_15 == $S_11)"

 local na = "($S_12 != $S_8)"
 local nb = "&($S_13 != $S_9)"
 local nc = "&($S_14 != $S_10)"
 local nd = "&($S_15 != $S_11)"

* build logical if:
 local cnt = ""
 if index("`1'","A") > 0 {local cnt = "`a'"}
                    else {local cnt = "`na'"}
 if index("`1'","B") > 0 {local cnt = "`cnt'" + "`b'"}
                    else {local cnt = "`cnt'" + "`nb'"}

  if `2' > 2 {
    if index("`1'","C") > 0 {local cnt = "`cnt'" + "`c'"}
                    else {local cnt = "`cnt'" + "`nc'"}
  }

 if `2' > 3  {
    if index("`1'","D") > 0 {local cnt = "`cnt'" + "`d'"}
                    else {local cnt = "`cnt'" + "`nd'"}
  }

 qui  replace `gen' = "`1'" if `cnt'
 qui  count if `cnt'
 global S_`4' = _result(1) 
 local siz = 16 - length(string(_result(1)))
 local t = ""
 di "`1'" _col(8) "|" _col(`siz') _result(1) _continue
 local p =  round((_result(1)/$S_6)*100,1)
 local siz = 6 - length(string(`p'))
 di _col(`siz') "`p'" " %  " "`cnt'" 
end