ar function

Set up AR(p) correlation structures