bbdaniels

power oneproportion 0.62 0.5 , n(50(1)150) graph(recast(line) plotopts(lw(thick) lc(red)) xline(100 131,lc(black)) ylab(0 "0%" 0.5 "50%" 0.8 "80%" 0.69 "69%" 1 "100%") yline(.69 .8,lc(black)) note("") title("") subtitle("") ) 

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)