S1i <- sapply(listw$weights, function(x) sum(x^2))

lx <- lag.listw(listw, x, zero.policy = zero.policy)

res <- (lx - Wi * xibar)

res <- res/sqrt(si2 * (((n - 1) * S1i - Wi^2)/(n - 2)))

程序2-matlab程序

x=input('x=')
weights=input('weights=')
n=length(x)
xibar=(sum(x)*ones(n,1)-x)/(n-1)
si2=((sum(x.^2).*ones(n,1)-x.^2)/(n-1))-xibar.^2
Wi=weights(:)'%这句对应Wi <- sapply(listw$weights, sum)是否正确？
Wi2=Wi.^2
S1i=sum(weights.^2)
lx=Wi*x %这句如何从lx <- lag.listw(listw, x, zero.policy = zero.policy)改写

res=lx-Wi*xibar

res =res./sqrt(si2 .* (((n - 1) * S1i - Wi.^2)/(n - 2)))