Mod9P6

library(vegan)

## Warning: package 'vegan' was built under R version 4.3.2

## Loading required package: permute

## Warning: package 'permute' was built under R version 4.3.2

## Loading required package: lattice

## This is vegan 2.6-4

dune.bio<-read.table("dune_bio.txt", sep="\t", header=T, row.names=1)
dune.env<-read.table("dune_env.txt", sep="\t", header=T, row.names=1)
dune.env.euc<-vegdist(decostand(dune.env[,c(1,2,5)], method="standardize"), method="euclidean")
dune.bio.bray<-vegdist(dune.bio, method="bray")
mantel(dune.env.euc, dune.bio.bray)

## 
## Mantel statistic based on Pearson's product-moment correlation 
## 
## Call:
## mantel(xdis = dune.env.euc, ydis = dune.bio.bray) 
## 
## Mantel statistic r: 0.5496 
##       Significance: 0.001 
## 
## Upper quantiles of permutations (null model):
##   90%   95% 97.5%   99% 
## 0.157 0.199 0.241 0.285 
## Permutation: free
## Number of permutations: 999

plot(dune.env.euc, dune.bio.bray)