similar to: RDA and trend surface regression

I'm trying to perform a RDA analysis in the vegan package using the "varpart" function. I have a matrix of community structure data in different sites (rows), which I want to explain using 3 matrices of environmental data. The 3 matrices are: water quality parameters; local land use variables; total land use variables. In each matrix, the number of the rows is the same (the

I am not familiar with the vegan package, so I am just making a guess here. If 'na.action=na.omit' is part of the call to varpart, try removing it from the function call and moving it above as follows: options(na.action="na.omit") RDA_Ger <- varpart(comm, x1, x2, x3, transfo="hellinger", scale = FALSE) Maybe that will help. Regards, Eric On Sat, Sep 16, 2017 at

Hello, I apologize for eventual double-posting. I am trying to integrate a 2-dimensional function that already calls the function adapt. More precisely, I am calling adapt(2,lower=c(-100,-100),upper=c(100,100),functn=function(s){1-exp(-50*Unc enteredGauss(c(-10,10,-10,10),60,s)}) where UncenteredGauss is given by the following code in R:

Hello I am using Jari Oksanen's CCA routine from the Vegan package on some estuary data, following a technique applied in (Anderson, M.J. & Gribble, N.A., 1998, Partitioning the variation among spatial, temporal and environmental components in a multivariate data set, Australian Journal of Ecology 23, 158-167). Some steps in the process require that the dependent matrix be constrained by

Hello I try to fit a LDA and RDA model to the same data, which has two classes. The problem now is that the training errors of the LDA model and the training error of the RDA model with alpha=0 are not the same. In my understanding this should be the case. Am I wrong? Can someone explain what the reason for this difference could be? Here my code: LDA model: =========== % x is a dataframe tmp =

Hi, Does anyone know how to export .RDA file (in R program) to a normal text file (readable by any text editor)? Also, how to export an object in R program into a text file (not .RDA file)? Thank you. Jittima

Hi Im using the VEGAN package to do an RDA ordination. In my plot I get my environmental scores as arrows/vectors, but my species scores as points. I would like to get the species scores as arrows as well. I tried the envfit function, and this creates arrows, but the lenght is not proportionate to the strengh of the correlation and the point is still there (so I both have point and arrow). fit

Hi, I have a set of gene expression data in .RDA file. I have downloaded Bioconductor and R program for analyzing these data. Anyway, I am not sure how to open this RDA file in R program (what is the command?) in order to look at these data. And which package should I use for analyzing it e.g. plot the expression image? Thank you. Jittima

Dear developeRs, my package FrF2.catlg128 holds large catalogues and is supposed to gain additional ones. All the catalogues are intended for the user. So far, the catalogues were stored in the data directory, and LazyData was "no". I understand that this is not considered wise any more (if it ever was), so that I want to change to LazyData "yes" with the next release

Hi everyone, sorry if this was discussed before (and in this situation, could you please point me to the discussion in the archive? My search didn't seem to be effective). Is there a way of getting the names of objects in a .rda file without having to load it? Thank you very much, benilton -- PhD Candidate Department of Biostatistics Bloomberg School of Public Health Johns Hopkins

I've tried a few things both with prcomp(), and rda() and its friends in vegan (including biplot.rda and ordiplot), but can't find a solution. I'd like to associate subsets of the points in a resulting biplot ("sites" in the rda object) with different plotting colors/text styles to emphasize certain sets of points. I can't figure out how to keep the arrows (for

Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are respected. What it means? Is it an error in the script? Commands and output: > permutest(rda.ind, perm=999, first=TRUE) Permutation test for rda  Call:

Dear List, assume I want to make a minor local change in a package that is supplied as .rda. For example, I want to get rid of the non-verbose-protected "Iteration" message in GLMM/lme4. Probably I have to load / change / save the package, but could someone help me to get the syntax right? Dieter Menne

Full_Name: Byron Ellis Version: 2.1.0 OS: OS X Submission from: (NULL) (67.124.246.46) > x = charToRaw("12345") > x [1] 31 32 33 34 35 > save(x,file="x.Rda") > rm(x) > load("x.Rda") > x [1] 00 00 00 00 00 Thats no fun.

Hello, I used the lda function from the MASS (VR) package and the rda function from the klaR package. I wanted to compare the result of this two functions by using the same training set. Thus, I used the rda function with lambda=1 an gamma=0, I should emulate the lda function and I should obtain the same result. But this it not the case, the two result are very different. My training set is

Hi, I have seen the answer to this sometime before but I just can't find it again - pointers appreciated. I have a package that contains some data.frames saved as .Rda files in the data/ directory. When the package is loaded I would like to have them be available in the workspace (without the user having to explicitly load them using data(...)). If my package does not use a NAMESPACE, I can

Dear users, I would like to process all the lists from all *.rda files that I have in one folder. Up to now, I can load all the *.rda files without any problem. The problem is when I want to access the list saved within each *.rda file (only one list per rda file). Here is my code: fpath <- "D:/R" listnames <- list.files(path=fpath, pattern=glob2rx("*.rda"),

Hello I'm triying to import data to R What is the problem? > es_rda <- read.table("D:\\fdrtrees\\data\\es.rda",sep=",",header=T) Warning message: incomplete final line found by readTableHeader on 'D:\fdrtrees\data\es.rda' Thank you a lot Rita [[alternative HTML version deleted]]

Hi, Can someone tell me how to create .rda data files in R so that they can be used contributed package? Thanks Deli Wang [[alternative HTML version deleted]]

Hi, to speed up computations in our RobASt family of packages, we use interpolation on a grid of precomputed values which we save together with the interpolating functions (results of splinefun essentially) in sysdata.rda in the R folder of our pkg. After adding grids for some more models, this file has grown considerably, even after application of tools::resaveRdaFiles. At the moment we are at