bbdaniels

// Part 1: De-biasing a parameter estimate using controls

  // Develop some data generating process for data X’s and for outcome Y, with some (potentially multi-armed) treatment variable and treatment effect. Like last week, you should strongly consider "simulating" data along the lines of your group project.

  // This DGP should include strata groups and continuous covariates, as well as random noise. Make sure that the strata groups affect the outcome Y and are of different sizes, and make the probability that an individual unit receives treatment vary across strata groups. You will want to create the strata groups first, then use a command like expand or merge to add them to an individual-level data set.

cap prog drop runregs
prog def runregs, rclass

syntax anything // sample size goes here (divided by ten)

bbdaniels

bbdaniels

global data "${wd}/ppol_528_GuzmanLidia/data/Data_exp"
local data : dir "${data}/xlsx/" files *
clear
  tempfile allData
  save `allData' , emptyok
cap prog drop infill
prog def infill
  syntax anything
 
  qui foreach var of varlist `anything' {

bbdaniels

  local graphs ""
  local x = 1
  foreach var of varlist ///
    diarrhea_history_duration tb_history_sputum ///
    malaria_history_fevertype pph_history_pph ///
    diabetes_history_numblimb {

      local title :  var lab `var'
      local graphs `"`graphs' "\``var''"  "'
      tempfile `var'

bbdaniels

local nb=`.Graph.plotregion1.barlabels.arrnels'
forval i=1/`nb' {
  di "`.Graph.plotregion1.barlabels[`i'].text[1]'"
  .Graph.plotregion1.barlabels[`i'].text[1]="`.Graph.plotregion1.barlabels[`i'].text[1]'%"
}
.Graph.drawgraph

bbdaniels

\usepackage{pdfpages} % Allow PDF insert for non-TeX elements
\includepdf[noautoscale=true, pages=1]{pdf/15633-Analytics_coverFINAL_print.pdf}
\includepdf[noautoscale=true, pages=2]{pdf/15633-Analytics_coverFINAL_print.pdf}

bbdaniels

sysuse auto.dta, clear

local style "s(0) w(1) lc(none) discrete frac"

tw ///
  (histogram rep78 if foreign == 1 ///
    , fc(black%50) barwidth(0.9) `style') ///
  (histogram rep78 if foreign == 0 ///
    , fc(red%50) barwidth(0.8) `style') 

bbdaniels

// Regression margins visualizer

cap prog drop regmargins
prog def regmargins

  // Syntax command
  syntax anything
    // Expected syntax is like regress
    // Second entry should be factor variable (i.) and we will visualize levels

bbdaniels

sysuse auto.dta , clear

levelsof foreign, local(levels)

tempfile posty
postfile posty str10(a) b c d using `posty' , replace

foreach level in `levels' {
  local l : label (foreign) `level'
  post posty ("`l'") (2) (3) (4)

bbdaniels

Asking for help with errors in statistical programming

No matter how experienced you are in statistical programming, you are certain to encounter errors in your code from time to time. You will recognize some of them and be able to solve them quickly, but you will never stop encountering unexpected new errors that you need help to resolve. Asking for help with statistical programming is hard. You need to give enough information to the person helping you so that they can produce the exact same problem you have