Dear all Sorry to bother you. I have a GLS model: modela<-gls(pop2~pop1-1 + ccs + pop1:ccs,data=data2,corr = corSpher(c(5, 0.5 ),form = ~ latitude + longitude, nugget = TRUE), method = "ML") I was wondering how I can work out how much variance is explained by the model AND by each term. Any help would be much appreciated. Debbie Russell ******************************************************************************* Debbie Russell (PhD student) Ecology and Evolution group 10.14 Miall Building Faculty of Biological Sciences University of Leeds Clarendon Way Leeds, LS2 9JT +44(0)113 3433107 +44(0)7842 702 561 d.j.russell05 at leeds.ac.uk