clear set obs 100 forvalues i = 1/100 { generate y`i' = runiform() generate rng`i' = c(rngstate) } clear set obs 100 forvalues i = 1/100 { set seed `i' generate y`i' = runiform() *generate rng`i' = c(rngstate) } quietly pwcorr y1-y100 matrix A = r(C) svmat A