library(RODBC)
library(ade4)
Env <-read.table("sang_en.txt")
Spec <-read.table("sang_sp.txt")
Tra <- read.table("sang_tr.txt")

coaspec <- dudi.coa(Spec, scannf = FALSE, nf = 10)
dudienv <- dudi.pca(Env, scannf = FALSE, nf = 8, row.w = coaspec$lw)
duditra <- dudi.acm(Tra, scannf = FALSE, nf = 6, row.w = coaspec$cw)
rlqout <- rlq(dudienv, coaspec, duditra, scannf = FALSE, nf = 4)
plot(rlqout)
summary(rlqout)
randtest.rlq(rlqout)
Qqual=Tra[,unlist(lapply(Tra,is.factor))]
sco.boxplot(rlqout$mQ[,1],Qqual)