I am performing pairwise spatial correlation of 6 rasters using the raster.modified.ttest()
function in R, as per this very useful post.
I have consistently had the following error with my data;
Error in SpatialPack::modified.ttest(x.var, y.var, sp::coordinates(x[nb[[i]], :
NA/NaN/Inf in foreign function call (arg 6)
and yet no errors with any dummy data. The following will yield a result;
s.pix1 <- as(raster(xmn=0,xmx=10,ymn=0,ymx=10,res=1,vals=sample(1:5,100,replace=T)), "SpatialGridDataFrame")
s.pix1[sample(1:100,25)]<-NA
s.pix2 <- as(raster(xmn=0,xmx=10,ymn=0,ymx=10,res=1,vals=sample(1:5,100,replace=T)), "SpatialGridDataFrame")
s.pix2[sample(1:100,25)]<-NA
corr <- raster.modifed.ttest(x=s.pix1, y=s.pix2)
plot(raster(corr), main="spatially adjusted raster correlation")
Very small sample areas of my data would also yield results. However forcing larger contiguous areas of NA data into the dummy data recreated the error;
s.pix1 <- as(raster(xmn=0,xmx=10,ymn=0,ymx=10,res=1,vals=sample(1:5,100,replace=T)), "SpatialGridDataFrame")
s.pix1[c(1:5,11:15,21:25,31:35,41:45)]<-NA
s.pix2 <- as(raster(xmn=0,xmx=10,ymn=0,ymx=10,res=1,vals=sample(1:5,100,replace=T)), "SpatialGridDataFrame")
s.pix2[c(1:5,11:15,21:25,31:35,41:45)]<-NA
corr <- raster.modifed.ttest(x=s.pix1, y=s.pix2)
Error in SpatialPack::modified.ttest(x.var, y.var, sp::coordinates(x[nb[[i]], :
NA/NaN/Inf in foreign function call (arg 6)
Could be a bug? @JeffreyEvans provided the excellent answer in the linked post
No comments:
Post a Comment