sum(dudi.pca(kasc2df(map), row.w = dataenfa1$pr/sum(dataenfa1$pr), scan=FALSE)$eig) which gives you 5 only. The weights should sum to 1 (i.e. they are proportions). But then, how would you interpret this? This is the same as for the "global specialization".
Yes, but be careful: kasc2df returns a list with one component "tab" (the component of interest) and one component "index" (the component allowing to rebuild the original kasc), so that the correct code is:
sum(dudi.pca(kasc2df(map)$tab, row.w = dataenfa1$pr/sum(dataenfa1$pr), scan=FALSE)$eig)
Best, Clément Calenge _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo