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