similar to: Customized LDA (MASS) object plot

Hi list, I was searching how to properly write a command line for a mixed ANOVA. Well honestly, there are so many material on the older post of the list that just confused me. I have five factors. Season (fixed) Beach (fixed) Line (fixed) Block (random) Strata (random) nested in Block And for each of the tree strata per block I got 3 replicates. I saw lots of things about

Hi list members. I’ll try to plot the abundance means of nine transects as lines, with five points on each transect (A to I). I will also need to plot for each point, it’s standard deviation (once each point will have tree replicates) as whiskers. Another problem will be that all the points should exist for each transect, but we know that six of the nine transects will have blanks (A1, A2, A3, -,

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 list, I would like to know if there is a way in R to execute a post-hoc test (factor levels comparison, like Tukey for ANOVA) of a non-parametric analysis of variance (kruskal.test() function object). I was looking for something like the Dunn Test. I’ve searched the internet but nothing seemed to define a way to it. Do You have any suggestion? Other equivalent test maybe! Thanks in

Hi again List, Well this time I’m writing for a friend (really J). He needs to create a distance matrix based on an abundance matrix using the 1-Pearson’s R index. Well I told him to look at the proxy package, but there is only Pearson Index. He needs it to perform a clustering. Well, as soon as he told me there proxy only had the Pearson index I thought: “He could just do something like

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 list. I would like to know if someone had already successfully loaded and Arc/Info Ungenerated (.dat – extracted from http://rimmer.ngdc.noaa.gov/mgg/coast/getcoast.html) file on R. I’m able to do so with shape files using maptools functions, but these ungenerated data are making me crazy. I don’t have any of ArcGIS tools (one of them probably can do that), and I can’t find any free tool that

Hi List, here I go again. Well I need to save plotted objects as .eps using the postscript() function, well I can do that but all resulting object are perfect squared dimensions (x = y). I need a rectangular output something like x = 2y dimension. Is it possible? The I’m new to R and postscript image format! I’m doing this…

Hi List, I would appreciate any suggestion on how can I make a text I’ve inserted in a plot show some contrast? With this I mean that I have a white text on a plot and I would like to make a tiny border around it in black, so even being small sized and the entire graphic being small in the text page I’ll insert it. It will be visible. Well I tried inserting a box() function inside the text()

Hi everybody, Well here is another doubt! I’m elaborating some maps, and I pretend to plot over it the abundance of species. I could do that using points(), but I’d like to use bubble() from sp package. When trying to use bubble() I face the fallowing issue: When I call the bubble function it overwrite everything I previously plotted. It doesn’t happen with points (obviously), only with

Hi everybody, well I'm trying to use barplot for the first time. I get a suitable graphic but I'm having a weird problem. My graphic is generated from a matrix with several species (rowNames) and two abundance columns. So each species will have two bars, and the species name should appear at left. Until now everything is Ok, but when a I plot it the species names (my labels) are cut by

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

Hello, I have multiple data sets I would like to plot on the same XY scatterplot. The data sets have in common the same Y values. Could somebody tell me how to do this? I tried par(T=new) (I think this was it) but it literally overlays plots on each other. Thanks in advance. George

Hi list, recently, for didactic reasons, I had the need to change a cluster topology (just the order of my groups, defined in an height of choice), not the results, obviously. Well I’m looking at hclust or agnes (cluster package) objects for something that defines the line breakings of the graphics, so I could try something to change this plot ordering. But, I’m unable to figure that out. Is

Hi list, I’m trying to plot lines in the same graphic. There are lots of similar topic on the list history. I’ve read about xyplot (lattice) and other ways to do that but I only got complex error messages, and wasn’t able even to start a decent R scrip. Well I have some paired information (each H W ), were H is the y axis and W is the x axis, how could I plot all the four lines at the same image?

Hello, i have a xls file. I will read it in r, what library-command i use for this?? any ideas?? Thanks Michele [[alternative HTML version deleted]]

Try the function sqlFetch to import the data in the spreadsheet. jo<-odbcConnectExcel("i:\\all\\sorkinjohn\\stats\\silvermannatalie\\NEMOcomb ined06-24-09.xls",readOnly = TRUE) mo<-sqlFetch(jo,'Your Sheet Name or Number',colnames=F,rownames=F) mo Hope it helps. ------------------------------------------------------------- MSc. Rodrigo Aluizio Centro de Estudos do

Hi list, Well, this time I’ve a doubt with mapping generation. I was already able to read and plot shapefiles, plot point on this map. All this without any problems, but now I want to do something that I think, only Golden Software Surfer is capable of. I would like to plot a kriging result on the existing map (map script below). Well I looked for on the net, but I was not able to be sure of

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