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