Merge pull request #57 from r-causal/code-examples

Author: malcolmbarrett

slides/pdf/07-using-pscores.pdf

@@ -705,7 +705,7 @@ <h2>Matching in R (ATM)</h2>
 <h2><em>Your Turn 1</em></h2>
 <div class="cell" data-layout-align="center">
 <div class="cell-output-display">
-<div class="countdown" id="timer_f26e5c3e" data-update-every="1" tabindex="0" style="right:0;bottom:0;">
+<div class="countdown" id="timer_3820d331" data-update-every="1" tabindex="0" style="right:0;bottom:0;">
 <div class="countdown-controls"><button class="countdown-bump-down">−</button><button class="countdown-bump-up">+</button></div>
 <code class="countdown-time"><span class="countdown-digits minutes">10</span><span class="countdown-digits colon">:</span><span class="countdown-digits seconds">00</span></code>
 </div>
@@ -725,65 +725,82 @@ <h2>Propensity scores</h2>
 <h2>Target estimands: ATE</h2>
 <p>Average Treatment Effect (ATE)</p>
 <p><span class="math display">\[\Large w_{ATE} = \frac{Z_i}{p_i} + \frac{1-Z_i}{1 - p_i}\]</span></p>
+<div class="cell" data-layout-align="center">
+<div class="sourceCode cell-code" id="cb11"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb11-1"><a href="#cb11-1"></a>(Z <span class="sc">/</span> p) <span class="sc">+</span> ((<span class="dv">1</span> <span class="sc">-</span> Z) <span class="sc">/</span> (<span class="dv">1</span> <span class="sc">-</span> p))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
+</div>
 </section>
 <section id="target-estimands-att-atc" class="slide level2" data-background-color="#23373B">
 <h2>Target estimands: ATT &amp; ATC</h2>
 <div class="fragment">
 <p>Average Treatment Effect Among the Treated (ATT) <span class="math display">\[\Large w_{ATT} = \frac{p_i Z_i}{p_i} + \frac{p_i (1-Z_i)}{1-p_i}\]</span></p>
+<div class="cell" data-layout-align="center">
+<div class="sourceCode cell-code" id="cb12"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb12-1"><a href="#cb12-1"></a>((p <span class="sc">*</span> Z) <span class="sc">/</span> p) <span class="sc">+</span> ((p <span class="sc">*</span> (<span class="dv">1</span> <span class="sc">-</span> Z)) <span class="sc">/</span> (<span class="dv">1</span> <span class="sc">-</span> p))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
 </div>
-<div class="fragment">
+</div>
+</section>
+<section id="target-estimands-att-atc-1" class="slide level2" data-background-color="#23373B">
+<h2>Target estimands: ATT &amp; ATC</h2>
 <p>Average Treatment Effect Among the Controls (ATC) <span class="math display">\[\Large w_{ATC} = \frac{(1-p_i)Z_i}{p_i} + \frac{(1-p_i)(1-Z_i)}{(1-p_i)}\]</span></p>
+<div class="cell" data-layout-align="center">
+<div class="sourceCode cell-code" id="cb13"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb13-1"><a href="#cb13-1"></a>(((<span class="dv">1</span> <span class="sc">-</span> p) <span class="sc">*</span> Z) <span class="sc">/</span> p) <span class="sc">+</span> (((<span class="dv">1</span> <span class="sc">-</span> p) <span class="sc">*</span> (<span class="dv">1</span> <span class="sc">-</span> Z)) <span class="sc">/</span> (<span class="dv">1</span> <span class="sc">-</span> p))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
 </div>
 </section>
 <section id="target-estimands-atm-ato" class="slide level2" data-background-color="#23373B">
 <h2>Target estimands: ATM &amp; ATO</h2>
 <div class="fragment">
 <p>Average Treatment Effect Among the Evenly Matchable (ATM) <span class="math display">\[\Large w_{ATM} = \frac{\min \{p_i, 1-p_i\}}{Z_ip_i + (1-Z_i)(1-p_i)}\]</span></p>
+<div class="cell" data-layout-align="center">
+<div class="sourceCode cell-code" id="cb14"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb14-1"><a href="#cb14-1"></a><span class="fu">pmin</span>(p, <span class="dv">1</span> <span class="sc">-</span> p) <span class="sc">/</span> (Z <span class="sc">*</span> p <span class="sc">+</span> (<span class="dv">1</span> <span class="sc">-</span> Z) <span class="sc">*</span> (<span class="dv">1</span> <span class="sc">-</span> p))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
 </div>
-<div class="fragment">
+</div>
+</section>
+<section id="target-estimands-atm-ato-1" class="slide level2" data-background-color="#23373B">
+<h2>Target estimands: ATM &amp; ATO</h2>
 <p>Average Treatment Effect Among the Overlap Population <span class="math display">\[\Large w_{ATO} = (1-p_i)Z_i + p_i(1-Z_i)\]</span></p>
+<div class="cell" data-layout-align="center">
+<div class="sourceCode cell-code" id="cb15"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb15-1"><a href="#cb15-1"></a>(<span class="dv">1</span> <span class="sc">-</span> p) <span class="sc">*</span> Z <span class="sc">+</span> p <span class="sc">*</span> (<span class="dv">1</span> <span class="sc">-</span> Z)</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
 </div>
 </section>
 <section id="histogram-of-propensity-scores" class="slide level2">
 <h2>Histogram of propensity scores</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-14-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-19-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="ate" class="slide level2">
 <h2>ATE</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-15-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-20-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="att" class="slide level2">
 <h2>ATT</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-16-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-21-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="atc" class="slide level2">
 <h2>ATC</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-17-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-22-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="atm" class="slide level2">
 <h2>ATM</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-18-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-23-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="ato" class="slide level2">
 <h2>ATO</h2>
 
-<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-19-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
+<img data-src="07-using-pscores_files/figure-revealjs/unnamed-chunk-24-1.png" style="width:80.0%" class="r-stretch quarto-figure-center"></section>
 <section id="ate-in-r" class="slide level2">
 <h2>ATE in R <img data-src="img/hex/propensity.png" class="absolute" style="top: 0px; right: 0px; width: 140px; "></h2>
 <p><br></p>
 <p>Average Treatment Effect (ATE) <span class="math inline">\(w_{ATE} = \frac{Z_i}{p_i} + \frac{1-Z_i}{1 - p_i}\)</span></p>
 <div class="cell" data-layout-align="center">
-<div class="sourceCode cell-code" id="cb11" data-code-line-numbers="|4"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb11-1"><a href="#cb11-1"></a><span class="fu">library</span>(propensity)</span>
-<span id="cb11-2"><a href="#cb11-2"></a>df <span class="ot">&lt;-</span> propensity_model <span class="sc">|&gt;</span></span>
-<span id="cb11-3"><a href="#cb11-3"></a>  <span class="fu">augment</span>(<span class="at">type.predict =</span> <span class="st">"response"</span>, <span class="at">data =</span> nhefs_complete) <span class="sc">|&gt;</span></span>
-<span id="cb11-4"><a href="#cb11-4"></a>  <span class="fu">mutate</span>(<span class="at">w_ate =</span> <span class="fu">wt_ate</span>(.fitted, qsmk)) </span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
+<div class="sourceCode cell-code" id="cb16" data-code-line-numbers="|4"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb16-1"><a href="#cb16-1"></a><span class="fu">library</span>(propensity)</span>
+<span id="cb16-2"><a href="#cb16-2"></a>df <span class="ot">&lt;-</span> propensity_model <span class="sc">|&gt;</span></span>
+<span id="cb16-3"><a href="#cb16-3"></a>  <span class="fu">augment</span>(<span class="at">type.predict =</span> <span class="st">"response"</span>, <span class="at">data =</span> nhefs_complete) <span class="sc">|&gt;</span></span>
+<span id="cb16-4"><a href="#cb16-4"></a>  <span class="fu">mutate</span>(<span class="at">w_ate =</span> <span class="fu">wt_ate</span>(.fitted, qsmk)) </span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
 </div>
 </section>
 <section id="your-turn" class="slide level2">
 <h2><em>Your Turn</em></h2>
 <div class="cell" data-layout-align="center">
 <div class="cell-output-display">
-<div class="countdown" id="timer_7fd186af" data-update-every="1" tabindex="0" style="right:0;bottom:0;">
+<div class="countdown" id="timer_5ab3ca30" data-update-every="1" tabindex="0" style="right:0;bottom:0;">
 <div class="countdown-controls"><button class="countdown-bump-down">−</button><button class="countdown-bump-up">+</button></div>
 <code class="countdown-time"><span class="countdown-digits minutes">10</span><span class="countdown-digits colon">:</span><span class="countdown-digits seconds">00</span></code>
 </div>

slides/raw/07-using-pscores.html

@@ -242,6 +242,12 @@ Direct Adjustment
 Average Treatment Effect (ATE)
 
 $$\Large w_{ATE} = \frac{Z_i}{p_i} + \frac{1-Z_i}{1 - p_i}$$
+```{r}
+#| eval: false
+(Z / p) + ((1 - Z) / (1 - p))
+```
+
+
 
 ## Target estimands: ATT & ATC {background-color="#23373B"}
 
@@ -250,11 +256,20 @@ $$\Large w_{ATE} = \frac{Z_i}{p_i} + \frac{1-Z_i}{1 - p_i}$$
 Average Treatment Effect Among the Treated (ATT)
 $$\Large w_{ATT} = \frac{p_i Z_i}{p_i} + \frac{p_i (1-Z_i)}{1-p_i}$$
 
-. . .
+```{r}
+#| eval: false
+((p * Z) / p) + ((p * (1 - Z)) / (1 - p))
+```
+
+## Target estimands: ATT & ATC {background-color="#23373B"}
 
 Average Treatment Effect Among the Controls (ATC)
 $$\Large w_{ATC} = \frac{(1-p_i)Z_i}{p_i} + \frac{(1-p_i)(1-Z_i)}{(1-p_i)}$$
 
+```{r}
+#| eval: false
+(((1 - p) * Z) / p) + (((1 - p) * (1 - Z)) / (1 - p))
+```
 
 ## Target estimands: ATM & ATO {background-color="#23373B"}
 
@@ -263,11 +278,23 @@ $$\Large w_{ATC} = \frac{(1-p_i)Z_i}{p_i} + \frac{(1-p_i)(1-Z_i)}{(1-p_i)}$$
 Average Treatment Effect Among the Evenly Matchable (ATM)
 $$\Large w_{ATM} = \frac{\min \{p_i, 1-p_i\}}{Z_ip_i + (1-Z_i)(1-p_i)}$$
 
-. . .
+
+```{r}
+#| eval: false
+pmin(p, 1 - p) / (Z * p + (1 - Z) * (1 - p))
+```
+
+## Target estimands: ATM & ATO {background-color="#23373B"}
 
 Average Treatment Effect Among the Overlap Population
 $$\Large w_{ATO} = (1-p_i)Z_i + p_i(1-Z_i)$$
 
+
+```{r}
+#| eval: false
+(1 - p) * Z + p * (1 - Z)
+```
+
 ## Histogram of propensity scores
 
 ```{r}