similar to: discerning species by color in cca biplot

Dear List, I am trying to plot the result of cca using matplot but failed. Pls kindly help and thanks. Elaine The error message was error in xy.coords(x, y, xlabel, ylabel, log = log) : (list) object cannot be coerced to type 'double' code rm(list=ls()) library(vegan) library(MASS) # input richness birdrich

Dear List, I used spec and envi for cca in ade4. (both are data.frame) However, there is a message telling that error in if (nf > rank) nf <- rank R missing value in TRUE/FALSE Please kindly help how to modify the code below. Thank you. Elaine code rm(list=ls()) spec <-read.csv("c:/migration/M_R_20100718_winterM_spec_vegan.csv",header=T, row.names=1) dim(spec) spec[1,]

Dear List, I would like to obtain biplot containing arrows and variable labels. However, after checking the previous message and package CCA and function anacor, my attempt remains an open question. Please kindly advise code or package I might miss to achieve it. Thank you. Elaine [[alternative HTML version deleted]]

Dear List, I tried to do cca based on species data and environmental variables (formula instead of community data). However, there was an error saying row sums must be >0. I searched the previous related messages but found few solutions. Please kindly help and thank you in advance. code This is vegan 1.17-3 Warning message: package 'vegan' was built under R version 2.10.1

Dear List, I tried to conduct cca using csv data but failed. The message said that data.frame is required. Please kindly share how to convert a csv-imported file to a data.frame. Thank you. Elaine code rm(list=ls()) spec <-read.csv("c:/migration/M_R_20100718_winterM_spec_vegan.csv",header=T, row.names=1) dim(spec) spec[1,] envi

Hi, I?m doing cca for a community data set in R and I have made a biplot for my data. Otherwise everything seems to be allright but the biplot is so messy I can?t read it well enough or publish it. I would like to get the row numbers out of the plot: I want the species position and the environmental variables in it as vectors but not the station numbers. How do I get them out? I have a raw

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

Dear List, I am using the command plot to present the relationship between bird richness (Y axis) and elevation (X axis). However, I would like to observe the distributions of bird richness in different administrative areas (A, B, C, …., G) in this plot. For example, the dots in area A might fall in the upper right part of the plot, while those in area B might appear in the middle of

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

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 All I am not clear quite how one could use cca from package vegan and the associated predict.cca to predict species abundance from environmental data (or if this is possible in a generalised way). In other words, can one derive a cca object based on known community data and use that to predict e.g. species abundances in a different number of samples based on environmental data? The help

He everybody, Well I have a biplot CCA-like origined from plot.cca (vegan package). I need to rotate on y axis the lines and symbols of constrained and sites representation. If I do that on an image editor, I rotate everything, including titles, axes labels and positions. I just need to rotate the inner par and keep the variables names (constrained) and symbols in the new positions but with the

Hi everybody, I'm trying, without success, to avoid some variables labels to overlap on a capscale (vegan package) object graph. I was able to move the labels with the coding below, but it moves the entire vector... CAPobj<-CAPpotiFT for(i in c(2:8))CAPobj$CCA[[i]][,1]<-CAPpotiFTI$CCA[[i]][,1] CAPpotiFTI$CCA[[i]][,1][2]<-0.85 And here the graphic script

Hi, I've just started using R and am having some problems with CCA using vegan. I'm looking at abundance p/m2 (hence decimals) vs environmental variables and have been using http://ecology.msu.montana.edu/labdsv/R/labs/lab12/lab12.html to guide me through. My organism data looks like this: Sample "Species_1" "Species_2" "Species_3" etc Sample_1

Dear list, I am trying to predict species volume from bioclimatic data, I have various sites and I have a data frame with species volume and the corresponding bioclimatic data for each site. I read on a discussion forum that you can use ordination to predict species abundance (in my case volume) from 'new' climate data for sites where you do not know the abundance. Unfortunately I

Hi All, I am using canonical correspondence analysis to compare a community composition matrix to a matrix of sample spatial relationships and environmental variables. In order to parse out how much variance is explained purely by space (S/E) or the environment (E/S) I am using a conditional (partial) CCA. I want to test significance via Monte Carlo but I can not find a way to do this with a

Dear all, I would like to calculate a cca (package vegan) with species and environmental data. One of these environmental variables is cos(EXPOSURE). The problem: for flat releves there is no exposure. The value is missing and I can't call it 0 as 0 stands for east and west. The cca does not run with missing values. What can I do to make vegan cca ignoring these missing values? Thanks a lot,

Hello, This is Elaine. I am using R 3.0 to download package vegan but failed. The warning message is package ‘vegan’ successfully unpacked and MD5 sums checked Warning: unable to move temporary installation ‘C:\Users\elaine\Documents\R\win-library\3.0\file16c82da53b1b\vegan’ to ‘C:\Users\elaine\Documents\R\win-library\3.0\vegan’ I cannot find the folder \file16c82da53b1b\ below

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

Dear List, This is Elaine, a postgraduate studying in bird distributions in East Asia. I want to calculate Simpson dissimilarity index, based on a presence/absence matrix of bird species in islands in East Asia. (matrix row: 36 islands/matrix column: species ID) (R package vegan to make NMDS and R package betapart) In most papers using vegan for NMDS and betapart for dissimilarity