[R] Obtain coefficients of several nlme objects

[Kuma Raj]

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")