similar to: plot(my.procrustes.model) from library {vegan}

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)

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,

Check out textplot in the gregmisc package ... Federico Calboli <f.calboli@ucl.ac.uk> Sent by: r-help-bounces@stat.math.ethz.ch 07/13/2004 12:06 PM Please respond to f.calboli To: r-help <r-help@stat.math.ethz.ch> cc: Subject: [R] plotting a table together with graphs Dear All, I would like to ask how to add a table to a

[ This was posted to the R-packages list (which I moderate), but definitely doesn't belong there. Martin Maechler ] ------- start of forwarded message ------- From: Federico Calboli <f.calboli at ucl.ac.uk> To: r-packages at stat.math.ethz.ch Subject: how to upload Date: 15 Jul 2004 18:30:26 +0100 Dear All, I just finished a add on lib called Malmig, to calculate Malecot

Dear All, I am in the somewhat unfortunate position of having to reproduce the results previously obtained from (non-metric?) MDS on a "kinship" matrix using Statistica. A kinship matrix measures affinity between groups, and has its maximum values on the diagonal. Apparently, starting with a nxn kinship matrix, all it was needed to do was to feed it to Statistica flagging that the

Dear All, how do I compute the left eigenvector of a matrix? I gather that "eigen" computes the right eigenvectors... Regards, Federico Calboli -- ================================= Federico C. F. Calboli PLEASE NOTE NEW ADDRESS Dipartimento di Biologia Via Selmi 3 40126 Bologna Italy tel (+39) 051 209 4187 fax (+39) 051 251 208 f.calboli at ucl.ac.uk

Dear All, I would like to ask how to get the Sum of Squares from fitted lme model. I appreciate that lme maximises the likelihood (or REML) and uses likelihood ratio tests, but I just fail why I could not get the SS if I want them. I could use lm to calculate them, but it looks quite pointless to fit a second (and wrong in a way) model for something so trivial. I know that the issue has been

Dear All, I apologize if the following has been already asked, but I could not find anything in the archives. I am running Debian Sid (on an Athlon XP), with R 1.8.1.cvs.20040307-1 (it describes itself as "R 1.9.0 experimental" when I fire R up) and I am unable to load the library "nlme" (I apt-getted r-recommended, so the library is on the system) >library(nlme) Error in

Dear List, I am writing a custom panel function and xyplot method to plot the results of a procrustes analysis from the vegan package. I am having trouble getting the call to panel.arrows to work as I wish when conditioning. The attached file contains the function definitions for the xyplot method and the custom panel and prepanel functions I am using. This example, using data and functions from

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

> -----Original Message----- > From: Federico Calboli > Sent: Tuesday, January 20, 2004 5:40 PM > To: r-help > Subject: [R] matrix exponential: M^0 > > I would like to ask why the zeroeth power of a matrix gives me a matrix > of ones rather than the identity matrix: > > > D<-rbind(c(0,0,0),c(0,0,0),c(0,0,0)) > > D<-as.matrix(D) > > D >

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

How about this: > expm function(x,pow=2) { xd <- diag((eigen(x)$values^pow)) xm <- eigen(x)$vector %*% xd %*% t(eigen(x)$vector) return(xm) } > xa [,1] [,2] [1,] 2 1 [2,] 1 3 > expm(xa,pow=3) [,1] [,2] [1,] 15 20 [2,] 20 35 > xa %*% xa %*% xa [,1] [,2] [1,] 15 20 [2,] 20 35 > From: Federico Calboli <f.calboli at ucl.ac.uk>

Hi all, I wondered how I can get the translated, rotated and scaledmatrix in a procrustes analysis. This does not help: ############################## library (vegan) #defining the target matrix: a= cbind (c (1,2,3,4), c( 10,12,14,16)) #defining the matrix to be rotated: b= a; b[,2]= c(9,11,15,16) c=procrustes (a,b); g=fitted (c) # trying to get to the procrustes matrices: d= t

Dear All, I would like to ask if there is any package that would calculate Malecot migration matrix (or the Imaizumi or Pena implementations of the MMM). I did a site-wide search, and went trough the table of contents of a number of packages that seemed promising, but with no luck. Regards, Federico Calboli -- ================================= Federico C. F. Calboli PLEASE NOTE NEW

Dear all, I would like to ask you the following: I have data about different manuscripts, together with data about the presence/absence of copying errors, in the days when manuscript were really manuscripts. I would ideally like to use the data to draw a phylogenetic tree, so that I can infer which ms was copied from which. The error presence/absence is coded in binary format. The plan is to use

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20111207/0d0b9fcc/attachment.pl>

Dear all, I am putting together a business case for R and I would like to know if there exist statistics on the distribution of R in academia, research and commercial institutions. The same information about S-plus would be welcome. Please reply only if you have consistent and certain information. Thank you very much, giovanni

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

Can anyone point me in the right direction of figuring out what downweight() is doing? I am using vegan to perform CCA on diatom assemblage data. I have a lot of rare species, so I want to reduce the influence of rare species in my CCA. I have read that some authors reduce rare species by only including species with an abundance of at least 1% in at least one sample (other authors use 5% as a