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