similar to: panel.arrows problem in custom panel function

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,

Hi I have just (finally) started to poke around in R and wanted to analyse a stream fish dataset with 28 sites and 18 species. When trying to follow the Vegan manual to run nmds from distance measures calculated by the vegdist function it turns out that I have two sites (streams) with the exactly the same four species (I have used pres-abs data in this case). When I try to run isoMDS I get an

Hi, I have been trying to create a function to generate a Procrustes plot, generated from package "vegan" in lattice. standard vegan code as follows library(vegan) pro=protest(P1, P8, permutations=4999,choices=1:4) plot(pro) Now, here is the code for the function that I have failed to get to work properly. panel.procrustes=function(x,y) }Pro=protest(x,y,permutations=4999,choices=1:4)

Hi, I am trying to create 3 mds plots: one with axis 1 vs axis 2, one with axis 2 vs axis 3, and one with axis 1 vs axis 3. When inputting my code, I only end up with one diagram and when inputting mdsg.mds$dims, the program returns 2 for 2 dimensions. How can I create the other two plots? Any help would be greatly appreciated, Calla Carbone The table I use is number of each taxa on each

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)

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)

Dear All, I would like to ask how to customize the graph corresponding to a procrustes analysis. I have to distance matrices, that I transform to two set of coordinates by means of muti dimensional scaling: library(mva) c1<-cmdscale(mat.dist1) c2<-cmdscale(mat.dist2) I vant to rotate c2 on c1, and I use the "procrustes" analysis from the {vegan} library. library(vegan)

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

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

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

Let's say I have a program that returns variables whose names may be any string within the vector NAMES=c("varA","varB","varC","varD","varE","varF"..."varZ"), but I do not ever know which ones have actually been created. So in one example output, "varA", "varC", and "varD" could exist, but

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

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

Dear R-list, I have a data set (in the following example called "a") which have: one "subject indicator" variable (called "id") three dependent variables (varD, varE, var F) three independent variables (varA, varB, varC) I want to fit 9 lme models, one per posible combination (DA, DB, DC, EA, EB, EC, FA, FB, FC). In stead of writting the 9 lme models, I want to

HI Guys, I know that this forum is not for homework but I am trying to interpret R output code. I was just wondering if someone might be able to help. I have been given the following. For (X1,X2) distributed bivariate normal with parameters mu1 = 5.8 mu2 = 5.3 sd1 = sd2 = 0.2 and p = 0.6 The r-code and inpit/output are as follows input m <- 5.3 + 0.6*(6.3 - 5.8) s <-

I don?t see how neither a) or b) applies to this question nor the technical merit of the remark about mixture models. Do you have a suggestion for a more appropriate forum for this issue/question? (stackoverflow basically sent me here). Kind regards > On 29. Jun 2017, at 16:58, Bert Gunter <bgunter.4567 at gmail.com> wrote: > > Offlist, because this is (a) an opinion and (b)

Dear All, I am using a Procrustes analysis to compare two NMDS ordinations for the same set of sites. One ordination is based on fish data, the other is based on invertebrate data. Ordinations were derived using metaMDS() from the {vegan} library as follows: fish.mds<-metaMDS(fish.data, distance="bray", k=3, trymax=100, wascores=TRUE, trace=TRUE, zero="add")

Would package "teigen" help? Ranjan On Thu, 29 Jun 2017 14:41:34 +0200 vare vare via R-help <r-help at r-project.org> wrote: > Hello! > > I am new to R (before used python exclusively and would actually call the R solution for this issue inside a python notebook, hope that doesn?t disqualify me right of the batch). > > Right now I am looking for a piece of

Offlist, because this is (a) an opinion and (b) about statistics and therefore offtopic. I don't know whether any such package exists, but I would predict that this is likely to be overdetermined (too many parameters) and therefore unlikely to be a successful strategy. Fitting a mixture of Gaussians is already difficult enough. Feel free to ignore, of course, and no need to reply. Cheers,