[R-sig-Geo] randtest.enfa {adehabitat}

Clément Calenge clement.calenge at gmail.com
Wed Jun 9 09:05:54 CEST 2010


>
> 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



More information about the R-sig-Geo mailing list