Hi,

I try to plot the prediction values from two fitted models and make a comparison between it, but sadly, it doesn't give me the expected graph when i try to run the plot command in part 2.  

##################
# Part 1
# This part works fine

> library(nlme)
> library(lattice)
> tmp <- groupedData(Temp ~ t2| Depth, data2, FUN = mean, 
> order.groups = TRUE,
> labels = list(x = "Month", y = "Temperature"),
> units = list(y = "(0C)"))

> lm.list1 <- lmList(tmp)
> lm.lme1 <- lme(lm.lis1)

####################

# Part 2
# The plot doesn't work as it gives the messege below this command.

plot(comparePred(lm.lis1, lm.lme1, length.out = 2), layout = c(4,3))


Error in augPred.lmList(object = list("-51" = list(coefficients = c(6.59358559883991,  : 
  augPred.lmList without "primary" can only be used with fits of groupedData objects

####################

Could someone give some ideas on how to sort out on this problem, please?

Thank you


      
	[[alternative HTML version deleted]]


