library(ade4)
Bscda <- discrimin(dudi.coa(Bs, scann = FALSE, nf = 1), Bsfac, scann = FALSE, nf = 1)
Eccda <- discrimin(dudi.coa(Ec, scann = FALSE, nf = 1), Ecfac, scann = FALSE, nf = 1)
Hicda <- discrimin(dudi.coa(Hi, scann = FALSE, nf = 1), Hifac, scann = FALSE, nf = 1)
Mgcda <- discrimin(dudi.coa(Mg, scann = FALSE, nf = 1), Mgfac, scann = FALSE, nf = 1)