data(columbus)
ind <- rownames(columbus)
GALfile <- tempfile("GAL")
write.nb.gal(col.gal.nb, GALfile, oldstyle=FALSE, ind="ind")
col.queen <- read.gal(GALfile, region.id=ind)
unlink(GALfile)
summary(diffnb(col.queen, col.gal.nb))
Run the code above in your browser using DataLab