*r
foreach type in u adj {
local aa = (corr_`type'_l1_r[1,1]+corr_`type'_l1_r[2,1]+corr_`type'_l1_r[3,1]+corr_`type'_l1_r[4,1])/4
local tr1 = (corr_`type'_l1_r[5,1]+corr_`type'_l1_r[6,1])/2
local bb = (corr_`type'_l1_r[7,1]+corr_`type'_l1_r[8,1]+corr_`type'_l1_r[9,1])/3
local tr2 = (corr_`type'_l1_r[10,1]+corr_`type'_l1_r[11,1])/2
local cc = (corr_`type'_l1_r[12,1]+corr_`type'_l1_r[13,1]+corr_`type'_l1_r[14,1])/3

matrix avg_r_`type' = `aa' \ `tr1' \ `bb' \ `tr2' \ `cc'

//2 lag
local aa = (corr_`type'_l2_r[2,1]+corr_`type'_l2_r[3,1]+corr_`type'_l2_r[4,1])/3
local tr1 = corr_`type'_l2_r[6,1]
local bb = (corr_`type'_l2_r[8,1]+corr_`type'_l2_r[9,1])/2
local tr2 = corr_`type'_l2_r[11,1]
local cc = (corr_`type'_l2_r[13,1]+corr_`type'_l2_r[14,1])/2

matrix avg_r_`type'_l2 = `aa' \ `tr1' \ `bb' \ `tr2' \ `cc'
}

matrix avg_m = avg_m_u , avg_m_adj
matrix avg_r =  avg_r_u, avg_r_adj

matrix avg_m_l2 = avg_m_u_l2 , avg_m_adj_l2
matrix avg_r_l2 =  avg_r_u_l2, avg_r_adj_l2

//matrix list avg_r
//matrix list avg_m

//matrix list avg_r_l2
//matrix list avg_m_l2

matrix rownames avg_r = "Stable 1" "Transition 1" "Stable 2" "Transition 2" "Stable 3"
matrix rownames avg_m = "Stable 1" "Transition 1" "Stable 2" "Transition 2" "Stable 3"
matrix colnames avg_r = Adjusted Unadjusted
matrix colnames avg_m = Adjusted Unadjusted

matrix rownames avg_r_l2 = "Stable 1" "Transition 1" "Stable 2" "Transition 2" "Stable 3"
matrix rownames avg_m_l2 = "Stable 1" "Transition 1" "Stable 2" "Transition 2" "Stable 3"
matrix colnames avg_r_l2 = Adjusted Unadjusted
matrix colnames avg_m_l2 = Adjusted Unadjusted


//PAPER TABLE 3
estout matrix(avg_m, fmt(2)), delimit(",") mlabels(none)
estout matrix(avg_r, fmt(2)), delimit(",") mlabels(none)

//PAPER TABLE 4
estout matrix(avg_m_l2, fmt(2)), delimit(",") mlabels(none)
estout matrix(avg_r_l2, fmt(2)), delimit(",") mlabels(none)