[Openspace] Re: problems with read.gwt2nb function

Hi everybody,



I'm using Geoda and R in order to calculate Moran I test for my data. I use
GeoDa to create my spatial weight and then R to run a code that allows me
calculate several Moran I test for different attributes of my data at once,
avoiding  tedious repetition. However it seems that the very same code that
worked pretty good for .gal files is not working for .gwt files. In fact,
there is a warning message after the read.gwt2nb command, but I cannot
figure it out. Does anybody can shed any light?

library(spdep)
#write de zero's matrix where the final results will be written.
ba.gwt.results <- matrix(0,nrow=39,ncol=5)
n <- nrow(ba.gwt.results)
ba <- read.csv(file = "d:/BrasilDep2006/BA BASE PRONTA/ba_base_toda.csv",

header=T)

ba.gwt.nb <- (file = "d:/BrasilDep2006/BA BASE PRONTA/ba.gwt")
ba.gwt.nb <- read.gwt2nb(ba.gwt.nb)

Warning messages:
1: region.id not named ID_ in: read.gwt2nb(ba.gwt.nb)

2: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39,
40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58,
59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77,
78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96,
97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112,
113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127,
128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142,
143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172,
173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187,
188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202,
203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217,
218, 219, 220, 221, 22 in: read.gwt2nb(ba.gwt.nb)

3: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39,
40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58,
59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77,
78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96,
97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112,
113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127,
128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142,
143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172,
173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187,
188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202,
203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217,
218, 219, 220, 221, 22 in: read.gwt2nb(ba.gwt.nb)

#Alternatively, one may use the funcion "read.gwt2nb",

# from the same library (spdep) to read .gwt objects.

#Create my weight?s matrix:

ba.gwt.listw <- nb2listw(ba.gwt.nb, style="W")

Error in nb2listw(ba.gwt.nb, style = "W") :

        Empty neighbour sets found

for(i in 1:n)

+ {

+ moran.test(ba[,i+2],ba.gwt.listw)$statistic->ba.gwt.results[i,1]

+ moran.test(ba[,i+2],ba.gwt.listw)$estimate[1]->ba.gwt.results[i,2]

+ moran.test(ba[,i+2],ba.gwt.listw)$estimate[2]->ba.gwt.results[i,3]

+ moran.test(ba[,i+2],ba.gwt.listw)$estimate[3]->ba.gwt.results[i,4]

+ moran.test(ba[,i+2],ba.gwt.listw)$p.value->ba.gwt.results[i,5]

+ }

Error in inherits(listw, "listw") : object "ba.gwt.listw" not found

Thank you very much,



Antonio Pedro Ramos.

	[[alternative HTML version deleted]]