RDA, Error in colMeans(x, na.rm = TRUE) : 'x' must be numeric, when data is numeric?

r vegan rda
user3420443 picture user3420443 · Mar 14, 2014 · Viewed 32.7k times · Source

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?

Answer

IRTFM picture IRTFM · Mar 14, 2014

Instead of using 'mode', you should be testing with 'class'. You probably have a factor column. They are of mode numeric but test FALSE with 'is.numeric'.

Related questions

Plotting envfit vectors (vegan package) in ggplot2

I am working on finalizing a NMDS plot that I created in vegan and ggplot2 but cannot figure out how to add envfit species-loading vectors to the plot. When I try to it says "invalid graphics state". The example below …

Read More
Plotting ordiellipse function from vegan package onto NMDS plot created in ggplot2

Instead of the normal plot function I am using ggplot2 to create NMDS plots. I would like to display groups in the NMDS plot using the function ordiellipse() from the vegan package. Example data: library(vegan) library(ggplot2) data(dune) # …

Read More
Plotting RDA (vegan) in ggplot

I'm still new to R, trying to learn how to use the library vegan, which I can easily plot in R with the normal plot function. The problem arises when I want to plot the data in ggplot. I know …

Read More