J'ai réalisé une analyse de variance hiérarchique sous R avec la fonction lm selon le script:
Code : Tout sélectionner
options(contrasts=c("contr.sum", "contr.sum"))
modele1=lm(Prix~Circuit+Enseigne%in%Circuit+Marque, data=prix,subset=UB=="Cassolettes Barquette Etui")
anova(modele1,type="III")
summary(modele1)
Le probème concerne le listing. Je n'arrive pas à savoir a quels modalité correspondent les lignes (dernière partie (circuit:enseigne)):
Code : Tout sélectionner
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.84990 0.15582 31.125 2,00E-16
Circuit1 -0.18494 0.09903 -1.868 0.0647
Marque1 -0.72496 0.30311 -2.392 0.0186
Marque2 -2.04484 0.41638 -4.911 3.42e-06
Marque3 -0.35677 0.17932 -1.990 0.0493
Marque4 1.12153 0.17334 6.470 3.36e-09
Circuithyper:Enseigne1 -3.72648 0.48815 -7.634 1.20e-11
Circuitsuper:Enseigne1 0.58041 0.68938 0.842 0.4018
Circuithyper:Enseigne2 NA NA NA NA
Circuitsuper:Enseigne2 -0.00043 0.32532 -0.001 0.9989
Circuithyper:Enseigne3 NA NA NA NA
Circuitsuper:Enseigne3 -0.39193 0.27229 -1.439 0.1531
Circuithyper:Enseigne4 0.46969 0.36423 1.290 0.2001
Circuitsuper:Enseigne4 NA NA NA NA
Circuithyper:Enseigne5 0.51501 0.23042 2.235 0.0276
Circuitsuper:Enseigne5 NA NA NA NA
Circuithyper:Enseigne6 NA NA NA NA
Circuitsuper:Enseigne6 NA NA NA NA
Circuithyper:Enseigne7 NA NA NA NA
Circuitsuper:Enseigne7 NA NA NA NA
Sachant que j'ai deux circuits hypers et supers (à l'intérieur desquels ont retrouve des ensiegnes comme indiqué ci-dessous).
Code : Tout sélectionner
Supers Hypers
ATAC AUCHAN
CASINO CARREFOUR
CHAMPION CORA
INTERMARCHE GEANT
MONOPRIX LECLERC
SUPER U
MATCH
Normalement R utilise le classement par ordre alphabétique donc circuit1=hypers, circuit2=Supers et de même pour les enseignes. Le problème est que les coefficients ne "tombent pas au bon endroit": une ensiegne n'appartenant qu' un circuit, il est impossible d'avoir un coefficient pour une enseigne sur un circuit sur lequel elle n'intervient pas.
Savez-vous comment interpréter ce listing?
Désolé pour la présentation je n'arrive pas à faire mieux...
Merci d'avance,
Marion Bouet