similar to: Ask help

Hi Folks... Œbeen playing with this for a while, with no luck, so I¹m hoping someone knows it off the top of their head... Difficult to find this nuance in the archives, as so many msgs deal with read.csv! I¹m trying to read a data file with the following structure (a little piece of the actual data, they are actually csv just didn¹t paste with the commas): wavelength SampleA SampleB SampleC

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

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,

I am trying to use the cca/rda/capscale functions in vegan to analyse genetic distance data ( provided as a dist object calculated using dist.genpop in package adegenet) with geographic distance partialled out ( provided as a distance object using dist function in veganthis method is attempting to follow the method used by Geffen et al 2004 as suggested by Legendre and . FORTIN (2010). I

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

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

How should I analysis it in R ???? all the resposes variables are ordinal from 0 to 10. and the explanatory variable is a factor ... -- View this message in context: http://old.nabble.com/Multivariate-problems-.-.-.-with-200-resposes-variables-and-1-explanatory-variable-tp26522912p26522912.html Sent from the R help mailing list archive at Nabble.com.

Hi everybody, i'm trying to draw rarefaction curves to estimate a population size from genotyped faeces. I used the Gimlet software which gave me a script and a "rarefaction.txt" file. I've copied both files in the work directory of R. I changed library(nls) by library(stats) in the script. But now, i'm still unable to run it. If i ask to show error messages, the software

Hi List, I’m using the vegan function specaccum to produce a rarefaction curve. In the function’s help it says: “Function ‘specaccum’ finds species accumulation curves or the number of species for a certain number of sampled sites or individuals”. Well, I would like to finds this curve for individuals, but when I compute it the function (using the ‘rarefaction’ method) gives me Sites, Richness

Hi all, I used R (vegan package) to make rarefaction curves and I calculated the Chao index for each curve. However, the plateau is far from reached. What I want to do now is the following: Based on the Chao index, I want to extrapolate the curve so I get an x-value which gives me an estimation of the total number of clones I'd have to pick up and sequence in order to have a full coverage of

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

I'm working with capscale and permutest for the first time, and having trouble getting statistical analyses for more than one constraining variable. I've read the documentation, but setting first=FALSE or using by="axis" doesn't seem to be helping. capscale seems to be fine, I receive output for more than one constrained axis. What am I doing wrong?

Here's a question for ecologists on the r-help list-- I'm addressing this to ecologists in particular because they're most likely to be familiar with the equation in question but I'll be happy to discuss the problem with anyone who's willing to take a whack at it. I'm trying to write a function to calculate the large sample variance of species richness estimates by

Hi everybody, I'm not very experienced with R software. I have used it several times for some of the population genetics analyses. I have problem with executing one of the script. The script is created by another software called Gimlet and it is aimed to calculate rarefaction curve in R software. However, when I try to execute the script, it says: *Error in library(nls) : there is no package

Dear list members, I have a data.frame so shaped: Sector Quadrants Plot Sic Time Species1 Species2 Species3 .. Species-n 1 1 1 1 5 0 0 1 . 0 2 1 1 1 12 1 1 1 . 0 3 1 1 1 34 0 1 0 . 0 4 1 1 1 23 1 1 0 . 0 5 2 1 1 22 1 1 1 . 1 6 2 1 1 10 1 1 1 . 1 7 2 1 1 2 1 0 0 . 0 8 2 1 1 2 0 0 1 . 0 9 3 1 1 12 0 0 0 . 1 . . . . . . . . . .

I've been struggling to get a loop to work. I want to create a new variable in each loop with data from some function. for example: # part of the names for the variables to be created Frags <- c("F04", "F05", "F07", "F09", "F11", "F13", "F14", "F17", "F18", "F19", "P20",

Dear all community, My issue is related to the R package (made4) that permits me to calculate the RV coefficient of co-inertia. However, it is a theoretical question. And if I am not mistaken, the list Usenet groups sci.stat.consult is not currently active. Let me explain briefly: Through different microbiota datasets, I have plotted PCoA, db-RDA and sPLS-DA using 3 different types of

Dear r-Community, Step1: I would like to calculate a NMDS (package vegan, function metaMDS) with species data. Step2: Then I want to plot environmental variables over it, using function envfit. The Problem: One of these environmental variables is cos(EXPOSURE). But 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. Therefore I