Difference between revisions of "Randomization in Stata"

Jump to: navigation, search
Line 21: Line 21:


An example of the do-file is as follows:
An example of the do-file is as follows:
<code>
* Set the environment to make randomization replicable</br>
version 12.0  [SETS VERSION] </br>
isid unique_id, sort  [SORTS UNIQUE ID]  </br>
set seed 12345  [SETS THE RANDOM SEED FOR REPLICATION]  </br></br>
* Assign random numbers to the observations and rank them from the smallest to the largest</br>
gen random_number = uniform()  [GENERATES A RANDOM NUMBER BETWEEN 0 AND 1]  </br>
egen ordering = rank(random_number) [ORDERS EACH OBSERVATION FROM SMALLEST TO LARGEST]  </br></br>
* Assign observations to control & treatment group based on their ranks </br>
gen group = .  </br>
replace group = 1 if ordering <= N/2 [ASSIGNS TREATMENT STATUS TO FIRST HALF OF SAMPLE]  </br>
replace group = 0 if ordering > N/2  [ASSIGNS CONTROL STATUS TO SECOND HALF OF SAMPLE]  </br>


 
    * Set the environment to make randomization replicable</br>
</code>
    version 12.0  [SETS VERSION] </br>
    isid unique_id, sort  [SORTS UNIQUE ID]  </br>
    set seed 12345  [SETS THE RANDOM SEED FOR REPLICATION]  </br></br>
    * Assign random numbers to the observations and rank them from the smallest to the largest</br>
    gen random_number = uniform()  [GENERATES A RANDOM NUMBER BETWEEN 0 AND 1]  </br>
    egen ordering = rank(random_number) [ORDERS EACH OBSERVATION FROM SMALLEST TO LARGEST]  </br></br>
    * Assign observations to control & treatment group based on their ranks </br>
    gen group = .  </br>
    replace group = 1 if ordering <= N/2 [ASSIGNS TREATMENT STATUS TO FIRST HALF OF SAMPLE]  </br>
    replace group = 0 if ordering > N/2  [ASSIGNS CONTROL STATUS TO SECOND HALF OF SAMPLE]  </br>


== Randomization with Multiple Treatment Arms ==
== Randomization with Multiple Treatment Arms ==

Revision as of 14:31, 8 April 2019

This page describes how and why to use Stata for randomization of control and treatment assignment in an RCT. Common alternatives to using Stata for randomization include: (i) Using the Excel Rand command; (ii) Randomizing directly within a chosen electronic survey platform such as SurveyCTO; or (iii) randomization through a public lottery.

Why use Stata to randomize

Using Stata to randomize and then preloading the generated data file into the survey software is generally preferred to using Excel or randomizing within the electronic platform. The main advantages of using Stata for randomization are as follows:

  • The process is transparent and reproducible.
  • The researcher has more control of the process, allowing you to check randomization balance and add stratification variables if needed.
  • Since randomization in Stata is done before the survey takes place (as opposed to randomization through the survey platform. This provides an opportunity to double check the result of a randomization and fix bugs before using the software in the field.

Steps needed for replicability when randomizing in Stata

Here are a few steps that should be followed to create a reproducible randomization using Stata:

  • Make sure your dataset includes a unique ID [respondent ID, household number, etc.]. If one doesn't exist yet, you can create one using generate command.
  • While writing a do-file, pay close attention to the following things:
    • Set version. Setting Stata's version in a do file ensures that the randomization algorithm is the same, since the algorithm sometimes changes between Stata versions.
    • Set seed. This makes sure that the same random number is generated for the first observation, for the second observation, and so on, for every time the code is run.
    • Sort the data by the unique ID. The data should be sorted such that observations are in the same order every time the code is run. The most optimal situation is sorting using an ID variable which uniquely and fully identifies each observation.
  • Convert the random numbers into categorical variables for treatment or control status.

The end goal is to have a CSV format file containing the ID variable used for randomization and the categorical variables created from the random numbers generated. This dataset will be preloaded into SurveyCTO, so that after an enumerator enters the respondent ID at the start of a survey, the result of the randomization will be loaded for the form and can be used for various sections of the survey.

An example of the do-file is as follows:

   * Set the environment to make randomization replicable
version 12.0 [SETS VERSION]
isid unique_id, sort [SORTS UNIQUE ID]
set seed 12345 [SETS THE RANDOM SEED FOR REPLICATION]

* Assign random numbers to the observations and rank them from the smallest to the largest
gen random_number = uniform() [GENERATES A RANDOM NUMBER BETWEEN 0 AND 1]
egen ordering = rank(random_number) [ORDERS EACH OBSERVATION FROM SMALLEST TO LARGEST]

* Assign observations to control & treatment group based on their ranks
gen group = .
replace group = 1 if ordering <= N/2 [ASSIGNS TREATMENT STATUS TO FIRST HALF OF SAMPLE]
replace group = 0 if ordering > N/2 [ASSIGNS CONTROL STATUS TO SECOND HALF OF SAMPLE]

Randomization with Multiple Treatment Arms

Randomization with Stratification in Stata

The steps to randomize in Stata with stratification is as follows:



Back to Parent

This article is part of the topic Randomized Control Trials