Comportement étrange pour var ?

[Christophe Genolini]

Bonjour,

Les bras m'en sont tombés :

Code : Tout sélectionner

a <- as.factor(c("H","S","S","H","S"))
var(a)
[1] 0.3


R calcule la variance d'une variable nominale...
Une raison a cela ? Ou bien faut-il ajouter cela à la liste des bizarreries R-iennes?

Christophe

[Nicolas Péru]

Tu peux récupérer tes bras :) :

Code : Tout sélectionner

>contrasts(a)
  S
H 0
S 1

Et si on fait :

Code : Tout sélectionner

> var(c(1,0,0,1,0))
[1] 0.3

R calcule la variance de l'indicatrice.

C'est peut être encore discutable sur le fond...

[Samir Messad]

ou bien la fonction var force la variable en numerique :

Code : Tout sélectionner

a <- as.factor(c("H","S","S","H","S"))
var(as.numeric(a))
0.3


ce qui peut se discuter aussi...

[jean lobry]

Bonjour,

Ici, la documentation de R (2.10) ne me semble pas très "clean". Si nous voulons la changer c'est sur R-Dev qu'il faut en discuter (avec arguments à l'appui, of course).

Amicalement,

Jean

[Matthieu Stigler]

Le plus troublant c'est que si var(factor) donne un résultat numérique, cov(factor, factor) n'en donne pas...

Code : Tout sélectionner

a <- as.factor(c("H","S","S","H","S"))
b <- as.factor(c("S","S","S","H","S"))
var(a)
var(b)

cov(a,b)
cov(as.numeric(a), as.numeric(b))

S'agit-t-il d'une incohérence de R? Ou peut-on argumenter que cela fait sens de calculer la variance d'une variable qualitative mais pas la covariance?

J'ai en effet entendu parler de corrélation poliychorique (voir pkg polychor) pour des variables qualitatives?