Movatterモバイル変換

Replication data for: Iron Deficiency and Schooling Attainment in Peru (Chong et al, 2016)

Description

The data is taken from Chong et al. (2016), who study the effect of iron deficiency anemia(i.e., anemia caused by a lack of iron) on school-age children’s educational attainment and cognitive ability in Peru.

Usage

data("AEJapp")

Format

A data frame with 215 observations on the 62 variables.

Source

Chong, A., Cohen, I., Field, E., Nakasone, E., and Torero, M. (2016). Replication data for: Iron Deficiency and Schooling Attainment in Peru. Nashville, TN: American Economic Association [publisher], 2016. Ann Arbor, MI: Inter-university Consortium for Political and Social Research [distributor], 2019-10-12.doi:10.3886/E113624V1.

References

Chong, A., Cohen, I., Field, E., Nakasone, E., and Torero, M. (2016). Iron Deficiency and Schooling Attainment in Peru.American Economic Journal: Applied Economics, 8(4), 222–255.doi:10.1257/app.20140494.

Examples

data(AEJapp)

Plot Method for 'sreg' Objects

Description

Visualize estimated ATEs and confidence intervals for objects of classsreg.

Usage

## S3 method for class 'sreg'plot(  x,  treatment_labels = NULL,  title = "Estimated ATEs with Confidence Intervals",  bar_fill = NULL,  point_shape = 23,  point_size = 3,  point_fill = "white",  point_stroke = 1.2,  point_color = "black",  label_color = "black",  label_size = 4,  bg_color = NULL,  grid = TRUE,  zero_line = TRUE,  y_axis_title = NULL,  x_axis_title = NULL,  ...)

Arguments

Value

Invisibly returns the ggplot object. Called for its side effects (i.e., generating a plot).

Printsreg Objects

Description

Print the summary table of estimation results forsreg objects.

Usage

## S3 method for class 'sreg'print(x, ...)

Arguments

Value

No return value, called for side effects.

Examples

data <- sreg.rgen(n = 200, tau.vec = c(0.1), n.strata = 4, cluster = TRUE)Y <- data$YS <- data$SD <- data$DX <- data.frame("x_1" = data$x_1, "x_2" = data$x_2)result <- sreg(Y, S, D, G.id = NULL, Ng = NULL, X)print(result)

Estimate Average Treatment Effects (ATEs) and Corresponding Standard Errors

Description

Estimate the ATE(s) and the corresponding standard error(s) for a (collection of) treatment(s) relative to a control.

Usage

sreg(  Y,  S = NULL,  D,  G.id = NULL,  Ng = NULL,  X = NULL,  HC1 = TRUE,  small.strata = FALSE)

Arguments

Value

An object of classsreg that is a list containing the following elements:

• tau.hat: a1 \times |\mathcal A|vector of ATE estimates, where|\mathcal A| represents the number of treatments

• se.rob: a1 \times |\mathcal A|vector of standard errors estimates, where|\mathcal A| represents the number of treatments

• t.stat: a1 \times |\mathcal A|vector oft-statistics, where|\mathcal A| represents the number of treatments

• p.value: a1 \times |\mathcal A|vector of correspondingp-values, where|\mathcal A| represents the number of treatments

• CI.left: a1 \times |\mathcal A|vector of the left bounds of the 95% as. confidence interval

• CI.right: a1 \times |\mathcal A|vector of the right bounds of the 95% as. confidence interval

• data: an original data of the formdata.frame(Y, S, D, G.id, Ng, X)

• lin.adj: adata.frame representing the covariates that were used in implementing linear adjustments

• small.strata: aTRUE/FALSE logical argument indicating whether the estimators for small strata (e.g., matched pairs or n-tuples) were used

• HC1: aTRUE/FALSE logical argument indicating whether the small sample correction (HC1) was applied to the variance estimator

Author(s)

Authors:

Juri Trifonovjutrifonov@u.northwestern.edu

Yuehao Baiyuehao.bai@usc.edu

Azeem Shaikhamshaikh@uchicago.edu

Max Tabord-Meehanm.tabordmeehan@utoronto.ca

Maintainer:

Juri Trifonovjutrifonov@u.northwestern.edu

References

Bugni, F. A., Canay, I. A., and Shaikh, A. M. (2018). Inference Under Covariate-Adaptive Randomization.Journal of the American Statistical Association, 113(524), 1784–1796,doi:10.1080/01621459.2017.1375934.

Bugni, F., Canay, I., Shaikh, A., and Tabord-Meehan, M. (2024+). Inference for Cluster Randomized Experiments with Non-ignorable Cluster Sizes.Forthcoming in the Journal of Political Economy: Microeconomics,doi:10.48550/arXiv.2204.08356.

Jiang, L., Linton, O. B., Tang, H., and Zhang, Y. (2023+). Improving Estimation Efficiency via Regression-Adjustment in Covariate-Adaptive Randomizations with Imperfect Compliance.Forthcoming in Review of Economics and Statistics,doi:10.48550/arXiv.2204.08356.

Bai, Y., Jiang, L., Romano, J. P., Shaikh, A. M., and Zhang, Y. (2024). Covariate adjustment in experiments with matched pairs.Journal of Econometrics, 241(1),doi:10.1016/j.jeconom.2024.105740.

Bai, Y. (2022). Optimality of Matched-Pair Designs in Randomized Controlled Trials.American Economic Review, 112(12),doi:10.1257/aer.20201856.

Bai, Y., Romano, J. P., and Shaikh, A. M. (2022). Inference in Experiments With Matched Pairs.Journal of the American Statistical Association, 117(540),doi:10.1080/01621459.2021.1883437.

Liu, J. (2024). Inference for Two-stage Experiments under Covariate-Adaptive Randomization.doi:10.48550/arXiv.2301.09016.

Cytrynbaum, M. (2024). Covariate Adjustment in Stratified Experiments.Quantitative Economics, 15(4), 971–998,doi:10.3982/QE2475.

Examples

library("sreg")library("dplyr")library("haven")### Example 1. Simulated Data.data <- sreg.rgen(n = 1000, tau.vec = c(0), n.strata = 4, cluster = FALSE)Y <- data$YS <- data$SD <- data$DX <- data.frame("x_1" = data$x_1, "x_2" = data$x_2)result <- sreg(Y, S, D, G.id = NULL, Ng = NULL, X)print(result)### Example 2. Empirical Data.?AEJappdata("AEJapp")data <- AEJapphead(data)Y <- data$gradesq34D <- data$treatmentS <- data$class_leveldata.clean <- data.frame(Y, D, S)data.clean <- data.clean %>%  mutate(D = ifelse(D == 3, 0, D))Y <- data.clean$YD <- data.clean$DS <- data.clean$Stable(D = data.clean$D, S = data.clean$S)result <- sreg(Y, S, D)print(result)pills <- data$pills_takenage <- data$age_monthsdata.clean <- data.frame(Y, D, S, pills, age)data.clean <- data.clean %>%  mutate(D = ifelse(D == 3, 0, D))Y <- data.clean$YD <- data.clean$DS <- data.clean$SX <- data.frame("pills" = data.clean$pills, "age" = data.clean$age)result <- sreg(Y, S, D, G.id = NULL, X = X)print(result)### Example 3. Matched Pairs (small strata).data <- sreg.rgen(  n = 1000, tau.vec = c(1.2), cluster = FALSE,  small.strata = TRUE, k = 2, treat.sizes = c(1, 1))Y <- data$YS <- data$SD <- data$DX <- data.frame("x_1" = data$x_1, "x_2" = data$x_2)result <- sreg(Y = Y, S = S, D = D, X = X, small.strata = TRUE)print(result)

Generate a Pseudo-Random Sample under the Stratified Block Randomization Design

Description

The function generates the observed outcomes, treatment assignments, strata indicators, cluster indicators, cluster sizes, and covariates for estimating the treatment effect within the context of a stratified block randomization design under the covariate-adaptive randomization (CAR).

Usage

sreg.rgen(  n,  Nmax = 50,  n.strata = 10,  tau.vec = c(0),  gamma.vec = c(0.4, 0.2, 1),  cluster = TRUE,  is.cov = TRUE,  small.strata = FALSE,  k = 3,  treat.sizes = c(1, 1, 1))

Arguments

Value

An object that is a 'data.frame' withn observations containing the generated values of the following variables:

• Y: a numericn \times 1vector of observed outcomes

• S: a numericn \times 1vector of strata indicators

• D: a numericn \times 1vector of treatments indexed by\{0, 1, 2, \ldots\}, where\code{D} = 0 denotes the control

• G.id: a numericn \times 1vector of cluster indicators

• X: adata.frame with columns representing the covariate values for every observation

Examples

data <- sreg.rgen(n = 1000, tau.vec = c(0), n.strata = 4, cluster = TRUE)