I have several lme objects like the ones shown below and I wish to combine the coefficients and confidence intervals of fixed effects of several models. Is there a function that could do that job? m1 <- lme(mark1 ~ pm10 + temp + + age + gender + bmi + statin + smoke + dow + season , data = df , random = ~ 1 | id,na.action=na.exclude, method="ML") m2 <- lme(mark2 ~ pm10 + temp + + age + gender + bmi + statin + smoke + dow + season , data = df , random = ~ 1 | id,na.action=na.exclude, method="ML")