J'ai un petit échantillon de 10 palmiers (4 sains, 3 moyennement malades et 3 fortement malades). J'ai fait une ACP (avec ade4) et j'aimerais faire maintenant une AFD :
Code : Tout sélectionner
T<-read.table("Couvert.txt",header=TRUE)
T.num<-T[,sapply(T,is.vector)]
l<-T$Palmier
afc<-dudi.coa(T.num,scannf=FALSE,nf=3)
afcd<-discrimin(afc,fac=l,scannf=FALSE,nf=3)
sumvp<-sum(afcd$eig)
inertie<-afcd$eig/sumvp*100
inertie.c<-cumsum(inertie)
et jobtiens un résultat étrange puisque j'ai 9 axes qui ont tous une inertie de 11.11%... est-ce possible ou bien ais-je commis une erreur que je n'arrive pas à trouver dans mon code ?
Merci pour votre aide
SImon