Simulation: Big and small effects from RCTs

Simulation: Big effects from small RCTs

// Simulation 1: Big effects from small RCTs

cap mat drop results
qui forv iter = 1/1000 {

  clear
    set obs 10000

    gen x = rnormal() > 0
    gen e = rnormal() * 4
    gen y = 5 + 0*x + e

    reg y x
      mat temp = r(table)
      mat result = temp[....,1]
      mat result = result'

      mat results = nullmat(results) ///
        \ result

}

clear
svmat results , n(col)

// End of dofile

Simulation: Small effects from measurement error

// Simulation 2: Small effects from measurement errors

cap mat drop results
qui forv iter = 1/1000 {

  clear
  set obs 10000

  gen x = runiform() * 100
    gen me = rnormal() * 50
    gen temp = x + me
    gen x_obs = max(temp,0)
  gen e = rnormal() * 4
  gen y = 5 + 20*x + e

  reg y x
    mat temp = r(table)
    mat result = temp[....,1]
    mat result = result'

    mat results = nullmat(results) ///
      \ result

}

clear
svmat results , n(col)

// End of dofile