I want to perform a rda in R, using vegan.
My code looks like this:
species<- read.delim("springspecies1.txt", header=T)
envdata<- read.delim("springenv1.txt", header=T)
RDA <- rda(species~Temperature + Salinity + O2 + Phosphate + Nitrate + Silica, envdata, scale=T, na.action=na.omit)
and I get the error message:
Error in colMeans(x, na.rm = TRUE) : 'x' must be numeric
when I check my data I get:
sapply(species, mode)
Station Year Month S.marinoi C.tripos
"numeric" "numeric" "numeric" "numeric" "numeric"
P.alata P.seriata R.setigera C.pelagica D.confervacea
"numeric" "numeric" "numeric" "numeric" "numeric"
C.decipiens P.farcimen C.furca
"numeric" "numeric" "numeric"
There are no NA or blanks in my data set. But it seems that the species data set is the problem. I compiled a new data set with the species, but I get the same problem again. Any ideas?