Demo example for [ritest] in Stata

ritest-demo.do

// [ritest] example

// Create randomization program ------------------------------------------------
cap prog drop my_randomization
    prog def my_randomization
	
    syntax , [*]
	
    // Must drop the variable for [ritest]. Make sure to keep it saved somewhere
    cap drop treatment
	
    // For demo purpose we do 40%-40-20% randomization
    xtile temp = runiform() , n(5)
    qui recode temp ///
        (1/2=1 "Control")(3/4=2 "Treatment 1")(5=3 "Treatment 2") ///
        , gen(treatment)
		
        lab var treatment "Treatment Assignment"
        drop temp
	
end // -------------------------------------------------------------------------


// Load data and randomize -----------------------------------------------------
sysuse auto, clear

    my_randomization

    ta treatment
	
    ritest treatment _b[treatment] ///
        , samplingprogram(my_randomization) ///
        kdensityplot : reg price treatment

// End of dofile --------------------------------------------------------------