Package 'ThreeGroups'

Title: ML Estimator for Baseline-Placebo-Treatment (Three-Group) Experiments
Description: 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).
Authors: Holger L. Kern [aut]
Maintainer: Holger L. Kern <[email protected]>
License: GPL-2
Version: 0.21
Built: 2024-12-19 06:43:34 UTC
Source: CRAN

Help Index


ML Estimator for Baseline-Placebo-Treatment (Three-Group) Experiments

Description

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).

Details

Package: ThreeGroups
Type: Package
Version: 0.21
Date: 2015-09-09
License: GPL-2
LazyLoad: yes

Author(s)

Holger L. Kern, [email protected], http://www.holgerkern.com

References

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.


ML Estimator for Baseline-Placebo-Treatment (Three-Group) Experiments

Description

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).

Usage

ThreeG(data, starting = NULL)

Arguments

data

a numeric vector or list of length 10, containing the data in the form of counts:
the number of successes (observations with Y=1Y=1) in the baseline group (s_b),
the number of failures (observations with Y=0Y=0) in the baseline group (f_b),
the number of successes in the treatment group among compliers (s_t_c),
the number of failures in the treatment group among compliers (f_t_c),
the number of successes in the treatment group among non-compliers (s_t_nc),
the number of failures in the treatment group among non-compliers (f_t_nc),
the number of successes in the placebo group among compliers (s_p_c),
the number of failures in the placebo group among compliers (f_p_c),
the number of successes in the placebo group among non-compliers (s_p_nc),
the number of failures in the placebo group among non-compliers (f_p_nc).

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.

Details

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, [email protected]) with questions related to the article.

Value

ThreeG() returns a list with the following components:

starting

contains the starting values used by optim().

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 optim(), including the Hessian.

Author(s)

Holger L. Kern, [email protected], http://www.holgerkern.com

References

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.

Examples

#
# 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)