the technical problem of matrix inversion, but not the
problem of the meaning of the results, which is an open issue...
HTH,
Best regards,
Clément Calenge
--
Clément CALENGE
Cellule d'appui à l'analyse de données
Direction des Etudes et de la Recherche
Office national de la chasse et de la faune
each segment then use the overlay method to
obtain info in each of these polygon.
Arnaud
--
Clément CALENGE
Cellule d'appui à l'analyse de données
Direction des Etudes et de la Recherche
Office national de la chasse et de la faune sauvage
Saint Benoist - 78610 Auffargis
tel. (33
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
On 06/09/2010 09:05 AM, Clément Calenge wrote:
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
:
it is provided by scatter.enfa (see the paper cited above). So I do not
understand how the result of scatter.enfa could be inconsistent with the
biplot, since the result of scatter.enfa *is* the biplot.
Best,
Clément Calenge
--
Clément CALENGE
Cellule d'appui à l'analyse de données
Office national de la
is the dataframe containing the values of environmental
variables (columns) in each pixel (rows), and that wei is as above the
vector describing the utilization weight of each pixel, you can
calculate the tolerance with:
sum(dudi.pca(tab, row.w=wei, scan=FALSE)$eig)
HTH,
Clément Calenge
here...
Best regards,
Clément Calenge
--
Clément CALENGE
Cellule d'appui à l'analyse de données
Office national de la chasse et de la faune sauvage
Saint Benoist - 78610 Auffargis
tel. (33) 01.30.46.54.14
___
R-sig-Geo mailing list
R-sig-Geo
of 1/eigenvalues just before the eigenvalue P-3
(where P is the total number of eigenvalues), then it would be a good
idea to keep the last three axes. Then factorial maps and other tools
described on the help page and in the paper would help to interpret the
results.
Hope this helps,
Clément
,
Clément Calenge
On 05/04/2010 01:41 PM, Consuelo Hermosilla wrote:
Hello guys,
I have a couple of questions regarding enfa (adehabitat package).
The first is related to this message: In predict.enfa(enfa1,
octopus.hab$index, octopus.hab$attr) :
the enfa *is not mathematically optimal
- as.data.frame(attr([EMAIL PROTECTED], att))
should be changed to
ka - [EMAIL PROTECTED]
Many thanks for this comment. This will be included in the next version
which will be uploaded on CRAN before the end of the week.
Best wishes,
Clément CALENGE
10 matches
Mail list logo