Dear list,

I have fitted a  gls model with an spatial structure 
in the error. I was trying to model how the number of records is 
determined by factors like roads, rails, presence/absence of airports 
etc, so I have continuous and categorical variables. . I fitted this model using the nlme package. 

I
 would like to do a Partial residual plot. So, plotting the effects of the variables with main effects but 
accounting for the other variables. Is it possible to do it?  
I have 
read the plot.gls in the same package but I'm not really sure how could 
be the proper sintaxis to obtain the plot I want, and if it's possible to use this function to do it.

gls.rec
 <- gls(log(records+1)~airport+rails+roads+pop, data = 
subsam,correlation=corSpher(form=~X_Mol + Y_Mol, nugget=TRUE), 
na.action=na.omit) ## which is my model.

Then I have,
plot 
(gls.rec, log(records+1)~rails |airport+roads+pop) ## which I assume is 
plotting records vs rails accounting for effect of airport,  roads and 
population.

But it doesn't give me the plot I want. Any help or information will be deeply appreciated.

Thanks, cheers

Liliana.


      
	[[alternative HTML version deleted]]

