At 09:46 13/06/2003 +0200, Sylvain Dolédec wrote:
>sous ADE-4
>setosa 0.12176 0.14082 0.029556 0.010884
>versicolor 0.2611 0.0965 0.2164 0.038324
>virginica 0.39626 0.10192 0.2985 0.073924
>
>sous R
> > setosa versicolor virginica
> > 0.1242490 0.2664327 0.4043429
> > setosa versicolor virginica
> > 0.14368980 0.09846939 0.10400408
> > setosa versicolor virginica
> > 0.03015918 0.22081633 0.30458776
> > setosa versicolor virginica
> > 0.01110612 0.03910612 0.07543265
>
>Tiens cela ne donne pas tout à fait le même résultat. Ai-je fait une erreur
>de calcul où y-a-t-il des considérations théoriques sous-jacentes ?
Il y a des considération théoriques sous jacentes. ADE-4 calcule la
variance observée:
(1/n)*Somme(xi - moyenne(x))
Alors que R calcule l'estimation de la variance par
(1/n-1)*Somme(xi - moyenne(x))
Donc, pour retrouver les résultats d'ADE avec R:
> library(mass)
> data(iris)
> lapply(1:4,function(x) tapply(iris[,x],iris$Species, function(x)
(length(x)-1)*var(x)/length(x)))
[[1]]
setosa versicolor virginica
0.121764 0.261104 0.396256
[[2]]
setosa versicolor virginica
0.140816 0.096500 0.101924
[[3]]
setosa versicolor virginica
0.029556 0.216400 0.298496
[[4]]
setosa versicolor virginica
0.010884 0.038324 0.073924
Cordialement,
Clément.
======================================
Clément CALENGE
Doctorant / PhD.
UMR CNRS 5558 - Equipe "Ecologie Statistique"
Laboratoire de biométrie
Université Claude Bernard Lyon 1
43, Boulevard du 11 novembre 1918
69622 Villeurbanne Cedex
tel. (+33) 04.72.43.29.35
fax. (+33) 04.72.43.13.88
This archive was generated by hypermail 2b30 : Tue Sep 07 2004 - 13:45:22 MEST