J'ai un modèle de régression linéaire avec une variable catégorique (mâle et femelle) et une variable continue .
J'ai mis en place des codes de contrastes en R avec options(contrasts=c("contr.sum","contr.poly"))
. Et maintenant, j'ai des sommes de carrés de type III pour , et leur interaction (A: B) en utilisant .drop1(model, .~., test="F")
Ce que je suis coincé est de savoir comment la somme des carrés est calculé pour . Je pense que oui sum((predicted y of the full model - predicted y of the reduced model)^2)
. Le modèle réduit ressemblerait y~A+A:B
. Mais lorsque j'utilise predict(y~A+A:B)
, R renvoie des valeurs prédites qui sont les mêmes que les valeurs prédites du modèle complet. Par conséquent, les sommes des carrés seraient égales à 0.
(Pour les sommes de carrés de , j'ai utilisé un modèle réduit de , qui est le même que .)y~B+A:B
y~A:B
Voici un exemple de code pour des données générées aléatoirement:
A<-as.factor(rep(c("male","female"), each=5))
set.seed(1)
B<-runif(10)
set.seed(5)
y<-runif(10)
model<-lm(y~A+B+A:B)
options(contrasts = c("contr.sum","contr.poly"))
#type3 sums of squares
drop1(model, .~., test="F")
#or same result:
library(car)
Anova(lm(y~A+B+A:B),type="III")
#full model
predFull<-predict(model)
#Calculate sum of squares
#SS(A|B,AB)
predA<-predict(lm(y~B+A:B))
sum((predFull-predA)^2)
#SS(B|A,AB) (???)
predB<-predict(lm(y~A+A:B))
sum((predFull-predB)^2)
#Sums of squares should be 0.15075 (according to anova table)
#but calculated to be 2.5e-31
#SS(AB|A,B)
predAB<-predict(lm(y~A+B))
sum((predFull-predAB)^2)
#Anova Table (Type III tests)
#Response: y
# Sum Sq Df F value Pr(>F)
#(Intercept) 0.16074 1 1.3598 0.2878
#A 0.00148 1 0.0125 0.9145
#B 0.15075 1 1.2753 0.3019
#A:B 0.01628 1 0.1377 0.7233
#Residuals 0.70926 6