I'm trying to re-run a script I wrote earlier this year using the function
rda() in the vegan package. The script run fine back then, and I did not
change the dataset, so I was wandering whether there's some problem in a
updated version of the package (I re-installed R + all packages since then).
Thanks for any advice,
So here's the output I get:
rda(X = rel)
Prior probabilities of groups:
In is.na(x$error.rate) :
is.na() applied to non-(list or vector) of type 'NULL'
system x86_64, linux-gnu
svn rev 63987
version.string R version 3.0.2 (2013-09-25)
nickname Frisbee Sailing
To cite package ‘vegan’ in publications use:
Jari Oksanen, F. Guillaume Blanchet, Roeland Kindt, Pierre Legendre, Peter
R. Minchin, R. B. O'Hara, Gavin L. Simpson, Peter Solymos, M. Henry H.
Stevens and Helene Wagner (2013). vegan: Community Ecology Package. R
package version 2.0-10. http://CRAN.R-project.org/package=vegan
> Hello R-list
> Maybe someone knows what's going on here.
> I'm trying to re-run a script I wrote earlier this year using the function
> rda() in the vegan package. The script run fine back then, and I did not
> change the dataset, so I was wandering whether there's some problem in a
> updated version of the package (I re-installed R + all packages since then).
> Thanks for any advice,
> Lukas Kohl
Sorry for the late reply: I don't follow this list regularly.
A couple of points about your question:
(1) vegan indeed has function rda(), but its output has no resemblance
to your example (except for the line repeating Call:). Either you are using
some other package or you have made up your own version of rda() or you
do not show its output.
(2) The rda() function in vegan 2.2-0. This is documented in NEWS. You can
see this by issuing vegandocs("NEWS") command after loading library(vegan).
(3) I have no idea what your script does (id does something different than
vegan::rda()) and I cannot reproduce the problem without that knowledge.
Kind regards, Jari Oksanen
PS. Sorry for not top-posting: Gmane does not allow it.
PS2. Sorry for removing some of your message: Gmane requires this.
> So here's the output I get:
> > rda(rel)
> rda(X = rel)
> Regularization parameters:
> Prior probabilities of groups:
> Misclassification rate:
> apparent: %
> Warning message:
> In is.na(x$error.rate) :
> is.na() applied to non-(list or vector) of type 'NULL'
> There's no NA's in my dataset, which seems ok..
> > sum(is.na(rel))
>  0
> > nrow(rel)
>  59
> > ncol(rel)
>  49
> > head(rel)
> X14.0 i.15.0 ai.15.0 br.16.0 X15.1 n.15.0 i.16.0