This package implements the Maximum Likelihood estimator for baseline, placebo, and treatment groups (three-group) experiments with non-compliance proposed by Gerber, Green, Kaplan, and Kern (2010).

1 |

`data` |
a numeric vector or list of length 10, containing the data in the form of counts: |

`starting` |
a numeric vector of length 4, containing starting values for the 4 parameters to be estimated. If missing, suitable starting values are generated internally. |

The data, whether given as a list or a numeric vector, must be sorted in the order given above,
i.e., `(s_b)`, `(f_b)`, ... , `(f_p_nc)`.

In most cases there should be no need to specify starting values.

Note that the notation used here is slightly different from that presented in section 3.3 of Gerber, Green, Kaplan, and Kern (2010).

Please contact the maintainer of this package with questions related to this package and the corresponding author of the article (Alan Gerber, alan.gerber@yale.edu) with questions related to the article.

`ThreeG()`

returns a list with the following components:

`starting` |
contains the starting values used by |

`est` |
a 4x2 matrix containing estimates (first column) and standard errors (second column) of the proportion of compliers (row 1), the probability of success given treatment for compliers (row 2), the probability of success for non-compliers (row 3), and the complier average treatment effect (row 4). |

`optim.out` |
a list containing the output from |

Holger L. Kern, hkern@fsu.edu, http://www.holgerkern.com

Alan S. Gerber, Donald P. Green, Edward H. Kaplan, and Holger L. Kern. 2010. “Baseline, placebo, and treatment: Efficient estimation for three-group experiments.” *Political Analysis* 18 (3): 297-315.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | ```
#
# This example replicates the empirical analysis presented in
# Gerber, Green, Kaplan, and Kern (2010), Table 5, first column.
#
# input data
sstat <- list(
s_b = 51766,
f_b = 253182,
s_t_c = 858,
f_t_c = 3128,
s_t_nc = 650,
f_t_nc = 3812,
s_p_c = 752,
f_p_c = 3154,
s_p_nc = 672,
f_p_nc = 3779)
#
# let's estimate the complier average treatment effect
#
out <- ThreeG(sstat)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.