1whittermore.stat<-function(data, listw, zero.policy=FALSE) 2{ 3 n<-length(data$Observed) 4 r<-scale(data$Observed, center=FALSE, scale=sum(data$Observed) ) 5 6 T<-lag.listw(listw, r, zero.policy = zero.policy) 7 T<-(n/(n-1))*(t(r)%*%T) 8 9 return(T) 10} 11