similar to: ordisymbol - changing symbols used in plotting factor levels

Hi! Is there a way to get biplot and identify to work togheter. Having the output of prcomp I would like to draw a biplot that. Instead of plotting the sample (row-names) names plots some "pch" symbols. (thats easy with xlabs) But now I would like to add using identify the names to only some of the points. I have noticed that both biplot.prcomp and biplot.default does a lot of

Dear all, I'm performing a detrended correspondence analysis on vascular plant community data (296 species), and I have a question on the species scores projected in the ordination diagram. When I run a ordiplot all species are projected in the output graph, but I'd like to restrict the number of species plotted in the final graph. Some species are so rare in the data, that no relevant

I've tried a few things both with prcomp(), and rda() and its friends in vegan (including biplot.rda and ordiplot), but can't find a solution. I'd like to associate subsets of the points in a resulting biplot ("sites" in the rda object) with different plotting colors/text styles to emphasize certain sets of points. I can't figure out how to keep the arrows (for

Hi, First I should note I am relatively new to R so I would appreciate answers that take this into account. I am trying to perform an MDS ordination using the function ?metaMDS? of the ?vegan? package. I want to ordinate species according to a set of functional traits. ?Species? here refers to ?sites? in traditional vegetation analyses while ?traits? here correspond to ?species? in such

Dear all, I have 4 sites and want to determine how different they are from each other. For this I have decided to use R though it seems a bit daunting to learn. I have read data in from a CSV the structure is : Species1 Species2 Species3 Site1 4 4 7 Site2 3 1 0 Site3 0 99 6 Site4 75 3 33 There are many more species than shown above this is just an example. Here are the

I have looked all over the internet for being able to color sites differently in a plot of an MDS (metaMDS)- I would like to color the different sites in the ordination plot (plot or ordiplot). I have set the matrix up so that my site code is 1 .... 296 (first column). is there a way to make 1-23 blue, 24-40 red etc. thanks Stephen -- Let's not spend our time and resources thinking

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 R experts, I'm looking for some help with plotting vectors from envfit in vegan, onto a 3d plot using ordiplot3d. So far I have data.mds <- metaMDS(data, k=3,trace = FALSE) vect_data<-envfit(data.mds,vegdata[,3:21],choices=1:3,permu=9999) ordiplot3d(data.mds,envfit=vect_data) ordixyplot(data.mds,pch=pts,envfit=vect_data) (my data's not really called data, I thought it might be

Hi, I am working with a species-by-trait .csv file (columns=traits, rows=species) and get the following warning message when trying to plot results of both metaMDS and pcoa:  "Warning message: In ordiplot(x, choices = choices, type = type, display = display,  :   Species scores not available" I am using a Gower's transformation in both procedures within the metaMDS or pcoa

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

Hi all (especially R-core) I suppose, With the introduction of the new functional programming functions into base I thought I'd ask for a Curry() function. I use a simple one that looks this: Curry = function(FUN,...) { .orig = list(...);function(...) do.call(FUN,c(.orig,list(...))) } This comes in really handy when using say, heatmap():

Hi! Thanks for providing great help in R-related statistics. Now, however I'm stuck. I'm not a statistics person but I was recommended to use R to perform a nmds plot and PerMANOVA of my dataset. Sample(treatment) in the columns and species (OTU) in the rows. I have 4 treatments (Ambient Temperature, Ambient temperature+Low pH, High temperature, High temperature+low pH), and I have 16

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

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

Dear all, I sometimes use the following function: Curry <- function(FUN,...) { # by Byron Ellis, https://stat.ethz.ch/pipermail/r-devel/2007-November/047318.html .orig <- list(...) function(...) do.call(FUN,c(.orig,list(...))) } ... and have thought it might be convenient to have a method for [ doing this. As a simple example, > apply(M, 1, mean[trim=0.1]) # hypothetical

Hi, I successfully created depending sliders with prototype and slider.js. See the result here: http://www.soa-world.de/dev/slider/slider.html (here''s the code: http://pastebin.com/ma1e82eb) There is one problem with my code. What I do is I have a sliders class "function sliders()" inside I create the 5 sliders in a for loop and store them inside an array which is a member of

Hi is there a way to plot only the loadings in a biplot (with the nice arrows), and to skip the scores? thanks christoph

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 there Icecast, Hope all is well. I have a macOS application that takes audio in from an external mixer, and uses Apple core audio to process. I want to add a live streaming feature to an Icecast server. Is there any Swift API examples. I just want a basic API that will take the core audio and stream it out to an Icecast service. In my endeavors to find this, I have found the dark ice