This function will produce the estiamtes comparing to treatment 0 for the LME models that are fit. The Z values will match to SAS, which is what is used in decision making. This version requires glsFit requires glsFit to be a fit for gls
GetLSDiffGLS(glsFit, nTrt, nTime, dMAV, dTV, dLowerCI, dUpperCI, bPlacMinusTrt)