I'm trying to use Moran.test on a SpatialPolygonDataFrame consisting of 7194 elements in R. I know that there is around 150 polygons with NA values.
First I generate a spatial weights matrix:
W_cont_msoa moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T, na.rm=TRUE)Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T, : unused argument(s) (na.rm = TRUE)I've tried to go back to the original data and remove missing data using complete cases but I received the following:
> all.data[complete.cases(all.data)]Error in all.data[complete.cases(all.data)] : error in evaluating the argument 'i' in selecting a method for function '[': Error in complete.cases(all.data) : invalid 'type' (S4) of argumentThe structure of the data is as follows:
`Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots..@ data :'data.frame': 7194 obs. of 23 variables:.. ..$ MSOA01CD : Factor w/ 7194 levels "E02000001","E02000002",..: 1 2 3 4 5 6 7 8 9 10 ... .. ..$ MSOA01NM : Factor w/ 7194 levels "Adur 001","Adur 002",..: 1468 125 126 127 128 129 130 131 132 133 ... .. ..$ averagesp : num [1:7194] 11.6 20.9 20.3 21.4 19.1 ... .. ..$ nga : int [1:7194] 41 112 177 97 120 NA 142 161 122 113 ... .. ..$ mediansp : num [1:7194] 11.5 13.1 12.6 13.5 10 ... .. ..$ mbps2 : int [1:7194] 125 74 104 48 83 NA 64 67 76 25 ... .. ..$ msoa.total : int [1:7194] 2372 114 191 97 126 NA 157 208 130 119 ... .. ..$ p.nga : num [1:7194] 1.73 98.25 92.67 100 95.24 ... .. ..$ p.mbps2 : num [1:7194] 5.27 64.91 54.45 49.48 65.87 ... .. ..$ residents : Factor w/ 4057 levels "10,003","10,004",..: 2222 1748 27 1261 3216 NA 3380 403 3093 3256 ... .. ..$ area : Factor w/ 6908 levels "1,002.23","1,002.58",..: 3707 2869 2840 3304 839 NA 1899 1327 6811 748 ... .. ..$ density : num [1:7194] 25.5 31.3 46.9 24.8 72.1 NA 50.6 81.5 87.4 76.8 ... .. ..$ micro.bc : num [1:7194] 36 32 83 33 39 NA 74 60 44 32 ... .. ..$ small.bc : num [1:7194] 49 14 64 3 14 NA 50 14 14 73 ... .. ..$ medium.bc : num [1:7194] 850 0 15 0 0 NA 10 0 0 0 ... .. ..$ large.bc : num [1:7194] 245 0 0 0 0 NA 0 0 0 0 ... .. ..$ primary : num [1:7194] 60 0 0 0 0 NA 0 0 0 0 ... .. ..$ secondary : num [1:7194] 135 0 20 0 0 NA 0 0 5 5 ... .. ..$ tertiary : num [1:7194] 16600 95 335 100 145 NA 300 200 160 80 ... .. ..$ kibs : num [1:7194] 11805 20 80 20 30 ... .. ..$ edu : num [1:7194] 2205 2681 811 2658 116 ... .. ..$ p.non.white: num [1:7194] 21.4 35 45.4 19 33.7 NA 32.8 30.4 38.8 30.8 ... .. ..$ claimants : num [1:7194] 117 245 316 254 295 ... ..@ polygons :List of 7194 Can anyone suggest how I overcome this problem please?
Thanks,
Ed
أكثر...
First I generate a spatial weights matrix:
W_cont_msoa moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T, na.rm=TRUE)Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T, : unused argument(s) (na.rm = TRUE)I've tried to go back to the original data and remove missing data using complete cases but I received the following:
> all.data[complete.cases(all.data)]Error in all.data[complete.cases(all.data)] : error in evaluating the argument 'i' in selecting a method for function '[': Error in complete.cases(all.data) : invalid 'type' (S4) of argumentThe structure of the data is as follows:
`Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots..@ data :'data.frame': 7194 obs. of 23 variables:.. ..$ MSOA01CD : Factor w/ 7194 levels "E02000001","E02000002",..: 1 2 3 4 5 6 7 8 9 10 ... .. ..$ MSOA01NM : Factor w/ 7194 levels "Adur 001","Adur 002",..: 1468 125 126 127 128 129 130 131 132 133 ... .. ..$ averagesp : num [1:7194] 11.6 20.9 20.3 21.4 19.1 ... .. ..$ nga : int [1:7194] 41 112 177 97 120 NA 142 161 122 113 ... .. ..$ mediansp : num [1:7194] 11.5 13.1 12.6 13.5 10 ... .. ..$ mbps2 : int [1:7194] 125 74 104 48 83 NA 64 67 76 25 ... .. ..$ msoa.total : int [1:7194] 2372 114 191 97 126 NA 157 208 130 119 ... .. ..$ p.nga : num [1:7194] 1.73 98.25 92.67 100 95.24 ... .. ..$ p.mbps2 : num [1:7194] 5.27 64.91 54.45 49.48 65.87 ... .. ..$ residents : Factor w/ 4057 levels "10,003","10,004",..: 2222 1748 27 1261 3216 NA 3380 403 3093 3256 ... .. ..$ area : Factor w/ 6908 levels "1,002.23","1,002.58",..: 3707 2869 2840 3304 839 NA 1899 1327 6811 748 ... .. ..$ density : num [1:7194] 25.5 31.3 46.9 24.8 72.1 NA 50.6 81.5 87.4 76.8 ... .. ..$ micro.bc : num [1:7194] 36 32 83 33 39 NA 74 60 44 32 ... .. ..$ small.bc : num [1:7194] 49 14 64 3 14 NA 50 14 14 73 ... .. ..$ medium.bc : num [1:7194] 850 0 15 0 0 NA 10 0 0 0 ... .. ..$ large.bc : num [1:7194] 245 0 0 0 0 NA 0 0 0 0 ... .. ..$ primary : num [1:7194] 60 0 0 0 0 NA 0 0 0 0 ... .. ..$ secondary : num [1:7194] 135 0 20 0 0 NA 0 0 5 5 ... .. ..$ tertiary : num [1:7194] 16600 95 335 100 145 NA 300 200 160 80 ... .. ..$ kibs : num [1:7194] 11805 20 80 20 30 ... .. ..$ edu : num [1:7194] 2205 2681 811 2658 116 ... .. ..$ p.non.white: num [1:7194] 21.4 35 45.4 19 33.7 NA 32.8 30.4 38.8 30.8 ... .. ..$ claimants : num [1:7194] 117 245 316 254 295 ... ..@ polygons :List of 7194 Can anyone suggest how I overcome this problem please?
Thanks,
Ed
أكثر...