VCBART_cs function

Fit a VCBART model with compound symmetry error structure