Displaying 1 result from an estimated 1 matches for "cca_envi_20100815".
2010 Aug 14
1
cca biplot (vegan) failed in matplot
...e'
code
rm(list=ls())
library(vegan)
library(MASS)
# input richness
birdrich <-read.csv("c:/migration/CCA_richness_20100814.csv",header=T,
row.names=1)
# input environment
# GID 1 row sum = 0 therefore deleted GID1 in birdenvi
birdenvi <-read.csv("c:/migration/CCA_envi_20100815.csv",header=T,
row.names=1)
birdrich.cca <- cca(birdrich ~ ., birdenvi)
birdrich.cca
matplot(birdrich.cca)
[[alternative HTML version deleted]]