This package implements the Maximum Likelihood estimator for baseline, placebo, and treatment groups (threegroup) experiments with noncompliance 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 noncompliers (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 threegroup experiments.” Political Analysis 18 (3): 297315.
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.
All documentation is copyright its authors; we didn't write any of that.