Re: [R-sig-Geo] spplot(), size of plotting symbol in legend

2006-09-19 Thread karl . sommer
Hi Edzer, your suggestion worked well. I saved the plot as an object named z1. In the legend variable I found the settings for symbol size. In my case there were 4 entries which I changed from 0.1 to 1 using the following statement. z1$legend$inside$args$key$points$cex <- c(1,1,1,1) plot(z1) R

Re: [R-sig-Geo] Moran test on residuals of a logistic regression

2006-09-19 Thread Roger Bivand
On Tue, 19 Sep 2006, [EMAIL PROTECTED] wrote: > Hi > > I have a logistic regression model, and I would like to see if the > residuals of this model are spatially autocorrelated. I saw that the > lm.morantest function could be used to test for autocorrelation in > residuals from a linear model (lm

[R-sig-Geo] Moran test on residuals of a logistic regression

2006-09-19 Thread [EMAIL PROTECTED]
Hi I have a logistic regression model, and I would like to see if the residuals of this model are spatially autocorrelated. I saw that the lm.morantest function could be used to test for autocorrelation in residuals from a linear model (lm). Does a similar function exist for residuals of a logisti

Re: [R-sig-Geo] Sampling from a SAR error model

2006-09-19 Thread Roger Bivand
Sam, Indeed perplexing. I've tried to slim things down to: library(spdep) nsamp <- 100 rho <- 0.4 nb7rt <- cell2nb(7, 7, torus=TRUE) set.seed(20060918) e <- matrix(rnorm(nsamp*length(nb7rt)), nrow=length(nb7rt)) listw <- nb2listw(nb7rt) invmat <- invIrW(listw, rho=rho) y <- invmat %*% e lag_res0