Hi everybody,
Im using correlog in package ncf to evaluate spatial autocorrelation of
my data, following these steps:
setwd ("C:/Documents and Settings/
bio <- read.table("bio3.csv", header=TRUE, sep=",")
str(bio)
library(ncf)
wow<-correlog(t1$x, t1$y, t1$z, w = NULL, increment=50, resamp =0, latlon =
TRUE, na.rm = FALSE)
str(wow)
plot(wow)
By the way, bio3.csv contains:
y , x , z
10.0833333 , -72.7166667 , 81.9975
7.5 , -74.8666667 , 87.971
5.6166667 , -67.15 , 75.4197
5.983333 , -61.433333 , 82.4859
3.85 , -65.666667 , 86.6279
1.4333333 , -65.7833333 , 91.7385
3.65 , -65.7666667 , 87.4672
3.9 , -56.1666667 , 75.5698
8.833333 , -70.533333 , 85.1626
3.5 , -65.7166667 , 88
3.6166667 , -65.6833333 , 87.7646
6 , -62.8833333 , 86.6228
5.7166667 , -65.9 , 77.4138
8.1286111 , -71.4586111 , 81.8428
8.2555556 , -71.7955556 , 79.7814
8.0205556 , -71.7419444 , 82.6087
8.7105556 , -71.4291667 , 81.7585
2.8333333 , -59.5166667 , 81.6372
7.87 , -72.4377778 , 82.7694
4.5833333 , -61.45 , 87.2137
1.5833333 , -58.6333333 , 80.7423
9.2422222 , -70.2230556 , 80.2662
4.416667 , -61.583333 , 86.7048
5.1666667 , -59.466667 , 80.6327
1.5333333 , -58.8333333 , 81.2853
9.9833333 , -72.7166667 , 82.0076
4.6666667 , -59.8166667 , 83.1121
8.3630556 , -71.7722222 , 81.8428
2.8333333 , -58.95 , 80.3623
8.5047222 , -71.395 , 82.6823
6.8663889 , -63.3113889 , 82.8629
3.0402778 , -74.2213889 , 84.0164
5.5833333 , -63.5333333 , 85.9463
3.3833333 , -59.4666667 , 81.25
7.3166667 , -71.95 , 85.3675
9.7977778 , -70.0666667 , 82.8776
6.05 , -67.4166667 , 73.642
0.2 , -74.7666667 , 80.4054
9.099951 , -70.674869 , 82.1181
3.95 , -65.85 , 85.4167
6.5933333 , -62.65 , 83.6207
1.7166667 , -64.5833333 , 91.3717
4.5333333 , -67.5833333 , 78.0534
3.35 , -73.9333333 , 81.8871
6.1 , -74.8 , 92.1254
3 , -59.5166667 , 81.6964
4.85 , -73.2666667 , 80.2596
11.3 , -73.9525 , 82.1721
11.333611 , -74.033611 , 80.9722
2.6427778 , -75.9147222 , 90.9357
1.0091111 , -76.3602722 , 87.1142
6.5 , -66.6666667 , 71.9697
2.8166667 , -59.15 , 80.9942
1.4833333 , -65.7333333 , 92.029
3.1666667 , -65.8166667 , 87.6984
2.9333333 , -59.4833333 , 80.826
where y is latitude, x is longitude, and z is a bioclimatic variable.
I can get the correlograms, however their ranges sometimes (like bio3)
exceed the supposed limits of -1 and 1 of the Morans Index.
Can anybody tell me what am I doing wrong?
Thanks a lot!
Paolo
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
