construct_variance function

a function to create the estimated variance matrix from a PLMM fit