Ariel THANK YOU for posting this. worked great for me just pasting this part into R Ariel wrote
model.matrix.gls <- function(object, ...) model.matrix(terms(object), data = getData(object), ...) model.frame.gls <- function(object, ...) model.frame(formula(object), data = getData(object), ...) terms.gls <- function(object, ...) terms(model.frame(object),...)
and running ghlt again. it seemed to me like the output was correct. -- View this message in context: http://r.789695.n4.nabble.com/multiple-comparisons-with-generalised-least-squares-tp3441513p4648972.html Sent from the R help mailing list archive at Nabble.com.