Significance values for different distances
(moran.i.sign<-cbind(bin.dist, mean.dist, count, site.dist.i.p, vasc.pa.i.p, vasc.abd.i.p, vasc.pa.phy.i.p,
vasc.abd.phy.i.p,angio.pa.i.p, angio.abd.i.p, angio.pa.phy.i.p,
angio.abd.phy.i.p))
## bin.dist mean.dist count site.dist.i.p vasc.pa.i.p vasc.abd.i.p
## [1,] 96.74171 0.00 1926 0.000 0.000 0.000
## [2,] 198.61975 0.00 1924 0.000 0.000 0.152
## [3,] 272.69169 0.00 1926 0.000 0.000 0.542
## [4,] 336.51455 0.00 1924 0.000 0.000 0.356
## [5,] 395.40082 0.00 1926 0.000 0.006 0.923
## [6,] 456.59676 0.02 1924 0.023 0.072 0.717
## [7,] 517.76554 0.10 1926 0.097 0.065 0.215
## [8,] 583.27411 0.16 1924 0.156 0.399 0.988
## [9,] 655.11692 0.56 1926 0.557 0.995 0.575
## [10,] 733.56936 0.99 1924 0.988 0.551 0.819
## [11,] 821.20522 0.18 1926 0.181 0.413 0.449
## [12,] 921.78696 0.00 1924 0.002 0.015 0.952
## [13,] 1040.00858 0.00 1926 0.000 0.000 0.332
## [14,] 1185.66599 0.00 1924 0.000 0.000 0.335
## [15,] 1388.76471 0.00 1926 0.000 0.000 0.008
## [16,] 2079.99684 0.00 1924 0.000 0.000 0.000
## vasc.pa.phy.i.p vasc.abd.phy.i.p angio.pa.i.p angio.abd.i.p
## [1,] 0.000 0.074 0.000 0.622
## [2,] 0.112 0.013 0.021 0.865
## [3,] 0.112 0.107 0.010 0.959
## [4,] 0.086 0.839 0.443 0.235
## [5,] 0.006 0.058 0.572 0.233
## [6,] 0.151 0.048 0.257 0.799
## [7,] 0.606 0.456 0.575 0.932
## [8,] 0.983 0.810 0.376 0.507
## [9,] 0.404 0.365 0.036 0.192
## [10,] 0.542 0.207 0.409 0.209
## [11,] 0.008 0.200 0.875 0.963
## [12,] 0.796 0.097 0.897 0.664
## [13,] 0.000 0.032 0.984 0.700
## [14,] 0.207 0.004 0.033 0.910
## [15,] 0.000 0.001 0.000 0.897
## [16,] 0.002 0.073 0.000 0.776
## angio.pa.phy.i.p angio.abd.phy.i.p
## [1,] 0.858 0.001
## [2,] 0.808 0.017
## [3,] 0.669 0.528
## [4,] 0.699 0.367
## [5,] 0.299 0.150
## [6,] 0.131 0.070
## [7,] 0.224 0.062
## [8,] 0.515 0.545
## [9,] 0.948 0.791
## [10,] 0.641 0.735
## [11,] 0.045 0.233
## [12,] 0.544 0.201
## [13,] 0.791 0.601
## [14,] 0.591 0.674
## [15,] 0.987 0.150
## [16,] 0.677 0.041
Plots for spatial autocorrelation
Moran’s I for presence-absence data
#plots showing Moran's I for vascular
#dev.new(width=11.8, height=8)
par(mfrow=c(1,2))
plot(bin.dist, site.dist.i, type="l",axes=TRUE, cex=1, pch=16, col="grey15",ylab=expression(paste("Moran's ", italic("I"))), xlab="Distance (m)", las=1, cex.axis=1, cex.lab=1.2, ylim=c(-1,1),
cex.main=0.85, bty="l", lwd=3)
points(bin.dist,vasc.pa.i, type="l",col="gray65", lwd=3)
points(bin.dist,vasc.pa.phy.i, type="l",col="gray30", lwd=3)
abline(h=0, lwd=1, lty=2)
legend("bottomleft", c("SiteEnv", "VascBD", "VascPhBD"), col = c("grey15","grey65", "grey30" ), cex=1.1,
lty = c(1, 1, 1), pch = c(NA, NA, NA), lwd=3 , bg = "white", bty="n")
box(bty="l", lwd=3)
plot(bin.dist, site.dist.i, type="l",axes=TRUE, cex=1, pch=16, col="grey15",ylab=expression(paste("Moran's ", italic("I"))), xlab="Distance (m)", las=1, cex.axis=1, cex.lab=1.2, ylim=c(-1,1),
cex.main=0.85, bty="l", lwd=3)
points(bin.dist,angio.pa.i, type="l",col="gray65", lwd=3)
points(bin.dist,angio.pa.phy.i, type="l",col="gray30", lwd=3)
abline(h=0, lwd=1, lty=2)
legend("bottomleft", c("SiteEnv", "AngioBD", "AngioPhBD"), col = c("grey15","grey65", "grey30" ), cex=1.1,
lty = c(1, 1, 1), pch = c(NA, NA, NA), lwd=3 , bg = "white", bty="n")
box(bty="l", lwd=3)
#textClick("(A)", cex=1.45)
#textClick("(B)", cex=1.45)
Moran’s I for abundance data
#plots showing Moran's I for vascular
#dev.new(width=11.8, height=8)
par(mfrow=c(1,2))
plot(bin.dist, site.dist.i, type="l",axes=TRUE, cex=1, pch=16, col="grey15",ylab=expression(paste("Moran's ", italic("I"))), xlab="Distance (m)", las=1, cex.axis=1, cex.lab=1.2, ylim=c(-1,1),
cex.main=0.85, bty="l", lwd=3)
points(bin.dist,vasc.abd.i, type="l",col="gray80", lwd=3)
points(bin.dist,vasc.abd.phy.i, type="l",col="gray50", lwd=3)
abline(h=0, lwd=1, lty=2)
legend("bottomleft", c("SiteEnv", "VascBD", "VascPhBD"), col = c("grey15","grey80", "grey50" ), cex=1.1,
lty = c(1, 1, 1), pch = c(NA, NA, NA), lwd=3 , bg = "white", bty="n")
box(bty="l", lwd=3)
plot(bin.dist, site.dist.i, type="l",axes=TRUE, cex=1, pch=16, col="grey15",ylab=expression(paste("Moran's ", italic("I"))), xlab="Distance (m)", las=1, cex.axis=1, cex.lab=1.2, ylim=c(-1,1),
cex.main=0.85, bty="l", lwd=3)
points(bin.dist,angio.abd.i, type="l",col="gray80", lwd=3)
points(bin.dist,angio.abd.phy.i, type="l",col="gray50", lwd=3)
abline(h=0, lwd=1, lty=2)
legend("bottomleft", c("SiteEnv", "AngioBD", "AngioPhBD"), col = c("grey15","grey80", "grey50" ), cex=1.1,
lty = c(1, 1, 1), pch = c(NA, NA, NA), lwd=3 , bg = "white", bty="n")
box(bty="l", lwd=3)
#textClick("(A)", cex=1.45)
#textClick("(B)", cex=1.45)