My data is taken from an area where I do expect a certain degree of spatial correlation for each variable. My question is how to account for spatial correlation (defined with X and Y) in the GLM syntax?
Model1 <- glm(INC ~ Var1 + Var2 +Var3 + X:Y ,
data = MyData,
family = binomial(logit))