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)

> library(utils)

>

> #library utils is needed for str function.

>

> #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)

>

> #load and attach data:

> ba <- read.csv(file = "d:/BrasilDep2006/BA BASE PRONTA/ba_base_toda.csv",
header=T)

>

> #Create my .nb weight:

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

>

> #Read the .gwt file into a neighbors object:

> 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]]


