Gravar-mail: Using second-order generalized estimating equations to model heterogeneous intraclass correlation in cluster randomized trials