similar to: Problem installing package

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)),

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 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)

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 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 After trying a simple decorana analysis (from the vegan package) on a simple data frame which contains no NA's the following error was returned: > tt_decorana(covN) Error in decorana(covN) : NA/NaN/Inf in foreign function call (arg 1) Have any vegan users come across this error and know what can be done about it? Cheers, J

Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan documentation p.24: library(vegan) data(varespec) ## Bray-Curtis distances between samples dis <- vegdist(varespec) ## First 16 sites grazed, remaining 8 sites

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

Hello, I am interested in using the capscale function of vegan package of R. I already have a dissimilarity matrix and I am intended to use it as 'distance' argument. But then, I don't know what kind of data must be in 'comm' argument. I don't understand what type of data must be referred as 'species scores' and 'community data frame' since my data refer to

Hi, I have used Vegan to construct an NMDS ordination plot. I plotted sites of three forest types with the site number in it. My reviewer has asked me to use different symbols for each of the forest types. Can anyone send me how I can do this in R in simple steps. I have used the options like ordiplot, sel and pl syntaxes that are not working for the question that I asked for. Best, Sinu --

Hi, I'm using nmds command (library vegan) to analyze some fishing data. I'd like to plot not only points, but also the names of species and stations in a specified position. I used the command text(nmds$points[,1], nmds $points[,2],labels=row.names(nmds $points),pos=3,cex=0.5) But the labels are sometimes overlapped. Is there any way to use identify, or a similar command, to plot 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 been using princomp() recently, its very useful indeed, but I have a question about how to specify the rows of data you want it to choose. I have a set of variables relating to bird characteristics and I have been using princomp to produce PC scores from these. However since I have multiple duplicate entries per individual (each bird had a varying number of chicks), I only want

Dear all, Does anyone knows why the results of a BIOENV (PRIMER v. 6.1.15) are diferent of the bioenv() + mantel() in vegan? Not the spearman correlation, indeed the pseudo-p value. I know that the approach bioenv() + mantel() is biased. So, how the BIOENV (PRIMER) ends with larger p values (permutated). Acctualy how the permutation test in BIOENV (PRIMER) is conducted. The user guide does not

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

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

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:

Hi all, ok, i know i can cut a dendrogram, which i did. all i get is three objects that a dendrograms itself. for example: myd$upper, myd$lower[[1]], myd$lower[[2]] and so on. of course i can plot them seperately now. but the lower parts still have hundreds of branches. i?ll need a 30 " widescreen to watch the whole picture. what i?d like to is group the lower branches , so that i get a

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) :

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