Displaying 1 result from an estimated 1 matches for "braycurtis_nohead".
Did you mean:
braycurtis_noheader
2010 Jun 09
1
ordisurf (pkg vegan) gives implausible result
...68 51 31 54 32 58 50 37 35
35 34 52 29 53 24 37 50 62 51
[35] 28 39 47 41 49 83 33 50 55 71 59 50 57 47 46 49 43 30 56 23 49 35 31 27
38 52
The contour lines suggest that there are only two values below 32, when in
fact there are 11.
Any ideas? Thanks!
>>>
> bc <- read.table("braycurtis_noheader.dist",header=F)
> sample.mds <- cmdscale(bc)
> edaphic <- read.table("edaphic.txt", header=T)
> library(vegan)
This is vegan 1.17-2
> ordisurf(sample.mds, edaphic$K, xlab="PC 1", ylab="PC 2", main="K")
Loading required package: mgcv
T...