Skip to content
Prev 140546 / 398506 Next

assistance with RDAtest beta version application

Hi !!



Stanfield, Les (MNR) a ?crit :
Yon do not have to, but you can. I personally like to import the whole 
file I'm working (here "Table_11-3.txt") and go from there. If you feel 
more comfortable to do it another way there is no problem there !
OK
Here I get a bit confused

Here is what I propose:

### Load the whole file
tab11.3<-read.table("Table_11-3.txt",header=TRUE,row.names=1)

### Extract the species (Y) and the environmental variables (X) from tab11.3
Y<-tab11.3[,1:6]
X<-tab11.3[,7:10]

### Calculate the RDA
coco<-rdaTest(Y,X,testF=TRUE,nperm=999)
coco

### Draw the triplot
graph.rdaTest(coco,plot.type="F")

With these couple of lines you should get the same results as the ones 
presented in Legendre and Legendre (1998).
Have a nice day !!

Guillaume Blanchet