Mouches

Données

On s'intéresse à la relation entre l'activité sexuelle et la longévité
L'étude se base sur une expérience faite sur 5 groupes de 25 mouches mâles.
Aux mouches du premier groupe, on a fourni une femelle vierge par jour, et à celle du deuxième groupe une tous les huit jours.
Les groupes 3 et 4 correspondent à la même expérience mais avec des femelles récemment inséminées (ce qui rend le rapport impossible); les mâles du groupe 5 sont seuls.

Variables

Longev : longévité en jours (réponse)
Npart : nombre de partenaires (0, 1 ou 8)
Virg : vierge (v), inséminée (i), aucun (a, si nombre=0)
Thorax : longueur du thorax du mâle en m
Sleep : Pourcentage de chaque journée passé à dormir

Analyse

Il s'agit d'expliquer Longev en fonction de Npart, Virg, et Thorax.
Sur les 9 combinaisons possibles pour Npart et Virg, seules 5 ont un sens.
Pour éviter ce problème on peut retirer dans un premier temps le 5ième groupe de l'étude, ce qui fait un plan plus simple à 4 possibilités.
On observe que cela ne change quasiment pas les résultats.

1/ Lire des données.

Pour exclure le 5ième groupe, on pourra faire

flyn=fly[fly$virg!='a',];

2/ Faire l'analyse de variance. Quel modèle garde-t-on ?

3/ Refaire alors l'analyse avec les cinq groupes; on créera un nouveau facteur par aggrégation des facteurs virg et npart :

gr=paste(fly$virg,fly$npart,sep='');

fly$groupe=factor(gr);

Utiliser les résultats de l'analyse pour comparer les durées de vie moyennes dans les 5 groupes :
plutôt que d'utiliser directement les coefficients du modèle, on va faire un autre ensemble de 5 individus où toutes les mouches ont la même valeur de Thorax (p.ex. la valeur moyenne), mais (Npart,Virg) prend les 5 valeurs considérées:

m5=rep(mean(flyn$thorax),5);

cinq=data.frame(thorax=m5);

cinq$groupe=factor(c('i1','i8','v1','v8','a0'));

et utiliser "predict".

4/ Vérifier que "la variable Thorax explique 1/3 de la variance de Longev".

5/ Faire de même pour la variable Sleep au lieu de Longev.

Source

"Sexual Activity and the Lifespan of Male Fruitflies" par Linda Partridge and Marion Farquhar, Nature, 294, 580-581, 1981).