data(bryceveg) # returns a vegetation data.frame called bryceveg
data(brycesite) # returns an environmental data.frame called brycesite
dis.bc <- dsvdis(bryceveg,'bray/curtis') # returns a Bray/Curtis dissimilarity matrix
dis.elev <- vardist(brycesite$elev) # matrix of pair-wide differences in elevation
cor(as.vector(dis.bc),as.vector(dis.elev))
Run the code above in your browser using DataLab