Block cluster designs involve units that are grouped into clusters which are in turn grouped into blocks.

The clusters are assigned as a group to either treatment or control. When the clusters are of equal size the difference-in-means estimator yields unbiased but possibly imprecise estimates, especially when individuals are very similar within clusters and share characteristics that can potentially influence outcomes.

In addition (optionally) clusters are assigned within blocks. Ideally clusters within a block are similar. For instance blocks might be formed on the basis of pre-treatment covariates. Treatment is then assigned *within* each block. Two main objectives can be achieved with block randomization: first, we can reduce sampling variability when individuals that belong to a specific block have similar potential outcomes. Secondly, we can ensure that certain characteristics that we care about are present at a particular rate in our treatment and control samples. For instance, policymakers might be interested in having an equal number of women and men in treatment and control. By defining one block as men and another as women and conducting the randomization separately for each gender, we can ensure that exactly half of the men and half of the women are treated.

As a rule of thumb having blocks improves precision but having clusters weakens precision. Despite their drawbacks, cluster designs are sometimes unavoidable. For instance, imagine that we want to evaluate the effects of radio advertisements on donations to aid relief campaigns. It is, of course, impossible to target such advertising to a single individual, and the treatment must instead be assigned at a higher level; for instance at the city level.

## Design Declaration

**M**odel: We specify a population of`N`

units partitioned into`N_blocks`

(regions, say). Each block has`N_clusters_in_block`

clusters (cities, say) and within each cluster there are`N_in_cluster`

individuals. The variance of an individual’s outcomes has three components: the block’s variance, the cluster’s variance and the subject’s variance.**I**nquiry: We are interested in the average difference in individuals’ treatment and control outcomes. Formally, \(E[Y_i(Z = 1) - Y_i(Z = 0)]\).**D**ata strategy: The units of randomization are cities (clusters). Cities are blocked by region so that we ensure that equal numbers of cities from each region are included in control and in treatment. All individuals belonging to the same city are either treated or assigned to control.**A**nswer strategy: We estimate the average treatment effect using a linear model that accounts for block structure and clusters.

In code:

```
N_blocks <- 4
N_clusters_in_block <- 4
N_i_in_cluster <- 5
sd_block <- 0.5773
sd_cluster <- 2
sd_i_0 <- 0
sd_i_1 <- 0
rho <- 1
assignment_probs <- 0.5
control_mean <- 0
treatment_mean <- 0
verbose <- TRUE
population <- declare_population(blocks = add_level(N = N_blocks,
u_b = rnorm(N) * sd_block), clusters = add_level(N = N_clusters_in_block,
u_c = rnorm(N) * sd_cluster, cluster_size = N_i_in_cluster),
i = add_level(N = N_i_in_cluster, u_0 = rnorm(N) * sd_i_0,
u_1 = rnorm(n = N, mean = rho * u_0, sd = sqrt(1 -
rho^2)) * sd_i_1))
potential_outcomes <- declare_potential_outcomes(Y ~ (1 -
Z) * (control_mean + u_0 + u_b + u_c) + Z * (treatment_mean +
u_1 + u_b + u_c))
estimand <- declare_estimand(ATE = mean(Y_Z_1 - Y_Z_0))
assignment <- declare_assignment(block_prob = assignment_probs,
blocks = blocks, clusters = clusters)
reveal <- declare_reveal(Y, Z)
estimator <- declare_estimator(Y ~ Z, estimand = estimand,
model = lm_robust, fixed_effects = ~blocks, clusters = clusters)
block_cluster_two_arm_design <- population + potential_outcomes +
estimand + assignment + reveal + estimator
```

## Diagnosis

To better illustrate how different sources of variance differentially affect our power, we compare three designs. Each of them features a design with 50 blocks featuring two clusters each, with 10 individuals per cluster. In each design the effect size is .25 and the standard deviation in the outcome is 1 (when possible the designer defaults to a variance for the individual shock such that total variance =1). In the first, all of the variance derives from differences between individuals, in the second it derives from differences between clusters, and in the third it derives from variance between blocks.

```
individual_var <- block_cluster_two_arm_designer(N_blocks = 50,
N_clusters_in_block = 2,
N_i_in_cluster = 10,
sd_block = 0,
sd_cluster = 0,
ate = .25)
```

```
## [1] "The implied ICC in (control) is 0"
## [1] "The implied ICC in (control) conditional on block is 0"
```

```
cluster_var <- block_cluster_two_arm_designer(N_blocks = 50,
N_clusters_in_block = 2,
N_i_in_cluster = 10,
sd_block = 0,
sd_cluster = .99,
ate = .25)
```

```
## [1] "The implied ICC in (control) is 0.98"
## [1] "The implied ICC in (control) conditional on block is 0.98"
```

```
block_var <- block_cluster_two_arm_designer(N_blocks = 50,
N_clusters_in_block = 2,
N_i_in_cluster = 10,
sd_block = .99,
sd_cluster = 0,
ate = .25)
```

```
## [1] "The implied ICC in (control) is 0.98"
## [1] "The implied ICC in (control) conditional on block is 0"
```

`diagnosis <- diagnose_design(individual_var, cluster_var, block_var)`

Estimator Label | RMSE | Power | Coverage | Mean Estimate | SD Estimate | Mean Se |
---|---|---|---|---|---|---|

estimator | 0.06 | 0.96 | 0.94 | 0.25 | 0.06 | 0.06 |

(0.00) | (0.01) | (0.01) | (0.00) | (0.00) | (0.00) | |

estimator | 0.19 | 0.24 | 0.95 | 0.25 | 0.19 | 0.20 |

(0.01) | (0.02) | (0.01) | (0.01) | (0.01) | (0.00) | |

estimator | 0.00 | 1.00 | 0.99 | 0.25 | 0.01 | 0.01 |

(0.00) | (0.00) | (0.00) | (0.00) | (0.00) | (0.00) |

Diagnosis reveals that the efficiency is greatest when variance comes from blocks, and lowest when it comes from clusters. The reason for this is that our estimator models block-level variance, and can estimate the treatment effect independently from block-level effects because treatment assignment is orthogonal to blocks. However, any variance at the cluster level covaries with our treatment, and thus makes it harder to estimate the treatment effect.

## Using the Block Cluster Designer

In R, you can generate a block_cluster_two_arm design using the template function `block_cluster_two_arm_designer()`

in the `DesignLibrary`

package by running the following lines, which load the package:

`library(DesignLibrary)`

We can then create specific designs by defining values for each argument. For example, we create a design called `my_block_cluster_two_arm_design`

with our own parameter specifications by running the lines below.

```
block_cluster_two_arm_design <- block_cluster_two_arm_designer(N_blocks = 50,
N_clusters_in_block = 30,
N_i_in_cluster = 200,
sd_block = .5,
sd_cluster = .5,
sd_i_0 = .3,
sd_i_1 = .6,
prob = .5,
control_mean = 0,
treatment_mean = .2)
```

You can see more details on the `block_cluster_two_arm_designer()`

function and its arguments by running the following line of code:

`??block_cluster_two_arm_designer`

## Further reading

Gerber, Alan S., and Donald P. Green. 2012. *Field Experiments: Design, Analysis, and Interpretation*. New York: W.W. Norton.