similar to: question about capscale (vegan)

Hi list, I am using the capscale function in vegan_2.0-7 to do a constrained principal coordinates analysis, and I kept getting the following error message: Error in Y.r[, oo, drop = FALSE] : subscript out of bounds I googled but I couldn't find an answer. Could anyone tell me why this error msg and what to do? Here is the command I used:

Dear r-helpers, I just read in an article by Virtanen et al. (2006) where vegetation-environment relationships are studied by fitting smoothed surfaces on an NMDS ordination using GAMs (Wood 2000). The authors describe, that they used R? as goodness-of-fit statistic, which they compare to the R? of fitted vectors. Calculations were carried out using the package vegan (Oksanen). I know that I can

Dear all, I realised a correspondence analysis with function cca() of vegan library. Just like in Okansen (2010) in the example of R help: library(vegan) data(varespec) data(varechem) vare.cca<-cca(varespec~ Al + P + K, varechem) With plot.cca() function I represented the species matrix in the next way: plot(vare.cca,display="species") Being similar to: plot((c(-2,2)),(c(-2,2)),

I performed canonical correspondence analysis (cca) with the example data of vegan, but I'm not able to obtain a table like scores() for the constraining variables. I can see them in the summary() mode, but it would be great to have in a separate table. Any suggestion?, thanx Gianandrea require(vegan) data(varespec) data(varechem) vare.cca<-cca(varespec,varechem) scores(vare.cca)

Hi, I have a couple questions about fitting environmental (land use factors, plant species presence-absence, and soil variables) constraints to my CCA biplot. 1. After successfully plotting species and site scores in my CCA, I have been trying to insert the biplot arrows of the environmental constraints in my data set using the text() function. When I do that, the plot changes completely. Is there

Estimados, Buenas Tardes, Estoy teniendo problema para testar la significancia de los ejes del RDA. NO se cual seria el error. Alguien me podría ayudar? Desde ya muchas gracias. Saludos, Luis # Tests of all canonical axes anova.cca(ssp.rda.hel, by="axis", step=1000)#Para saber la significancia de cada eje Error in anova.cca(ssp.rda.hel, by = "axis", step = 1000) :

Hello This is a simple question but I couldn't google an answer. In the procrustes function of the vegan package, one uses plot(procrustes_object, kind=2) to obtain a plot of the residual differences. For instance: data(varespec) vare.dist <- vegdist(wisconsin(varespec)) library(MASS) mds.null <- isoMDS(vare.dist, tol=1e-7) mds.alt <- isoMDS(vare.dist,

Dear R-help, We are conducting a distance-based redundancy analysis using capscale and then testing for statistical significance for six terms in the model for the constrained ordination using anova.cca in the vegan package. The significance test is sequential, i.e., testing for significance of a term only after accounting for all preceding terms. Could someone please provide us with either the

Hello, I am using vegan to do an NMDS plot and I would like to suppress the labels for the loading vectors. Is this possible? Alternatively, how can I avoid overlap? Many thanks for the help. Example code: #perform NMDS using metaMDS() function spe.nmds<-metaMDS(data, distance='bray',k=2 , engine = "isoMDS", autotransform=F, trymax=1000) #calculate the loading (i.e.,

Hi all, I have site data with plant species cover and am looking for trends. I'm kind of new to this, but have done lots of reading and can't find an answer. I tried decorana (I know it's been replaced by ca.) and see a trend, but I'm not sure what it means. Is there a way to get the loadings/eigenvectors of the axes (like in PCA)? Is there a way to do this with rda() too? How

Dear All I am trying to do a partial canonical analysis of principal coordinates using Bray-Curtis distances. The capscale addin to R appears to be the only way of doing it, however, when I try and calculate a Bray-Curtis distance matrix either using Capscale or Vegedist (capscale I understand uses Vegedist anyway to calculate its distance matrix), R uses up all available memory on the computer,

Hello, I'm new to the installation R on Linux and I've followed the instructions on the R Installation and Administration Manual and on the FAQs. I also have changed the path at the .bash_profile file but when I try to type R at the shell prompt, it says: Fatal error: R home directory is not defined I've tried many things but all the time it says the same. Any help would be

Hello R experts, I've used metaMDS to run NMDS on some fish abundance data, and am also working on correlating environmental data to the NMDS coordinates. I'm fairly new to metaMDS and NMDS in general, so I have what are probably some very basic questions. My fish abundance data consists of 66 sites for which up to 20 species of fish were identified and counted. I ran metaMDS on this data

If you submit these lines, you end up with variable "vare.dis". I want to export vare.dis to csv. Stuck I am. library(vegan,logical.return = TRUE) #return=true verifies package is available library(MASS,logical.return=TRUE) #return=true verifies package is available data(varespec) #varespec is an example data file in the vegan package vare.dis <- vegdist(varespec)

I've been unable to find a R package that provides the means of performing Clarke & Ainsworth's BIO-ENV procedure or something comparable. Briefly, they describe a method for comparing two separate sample ordinations, one from species data and the second from environmental data. The analysis includes selection of the 'best' subset of environmental variables for explaining

Hi, I'm a beginner with R, but I'm getting excellent results with it. Well I've got an capscale object (vegan package) and I want to made a biplot with symbols representing six groups of areas. With the plot.cca function and some par attributes (like 'labels') I was able to substitute the samples names for keyboard symbols ('x','o','#',ect), but it's

Was there ever an answer to this question about how to do amova with factorial data? t hundley [[alternative HTML version deleted]]

Hello! I am eager to learn if I only have a data about the relationship between otu and sample, could I use the function capscale, and final make a point plot that x-axis is CAP1 and y-axis is CAP2? Besides, what function could I use to make the different rarefaction curve with different color in the a plot in R? Appreciate very much.

Dear List, I am writing a custom panel function and xyplot method to plot the results of a procrustes analysis from the vegan package. I am having trouble getting the call to panel.arrows to work as I wish when conditioning. The attached file contains the function definitions for the xyplot method and the custom panel and prepanel functions I am using. This example, using data and functions from

Hi, I'm trying to install de RSPerl module and i have some problems trying to test it. I've tried to contact with the author but the e-mail address seems not to exist. Hope that anyone in this list could explain me a little about this problem. I have a bash shell and what i do is the following: #PERLLIB=/usr/lib/R/library/RSPerl/share/blib/arch:/usr/lib/R/library/RSPerl/share/blib/lib #