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