R car package - извлечение остатков для двухстороннего анова с повторяющимися факторами
Я использую автомобильный пакет, чтобы выполнить двухстороннюю анову с повторными измерениями в R. Более конкретно, я использую функцию Anova() (обратите внимание на заглавную A). Кажется, это работает нормально, но я хочу проверить, являются ли остатки нормальными. Тем не менее, я не могу найти способ извлечь остатки или коэффициенты из вывода. Для нормальной линейной модели в R вы можете использовать остатки (мод) или график (мод), но здесь это не работает.
Любая помощь приветствуется!
Код для создания данных:
d=as.data.frame(t(matrix(c("AA","A",231,249,247,260,242,247,"BB","A",266,246,248,261,261,272,"CC","A",270,268,228,289,274,288,"DD","A",291,292,245,267,281,280,"AA","A",265,281,244,254,254,261,"BB","A",368,283,229,250,250,246,"CC","A",247,267,232,282,242,245,"DD","A",283,269,232,278,254,256,"AA","A",250,242,240,238,253,246,"BB","A",260,231,228,241,263,251,"CC","A",276,272,236,282,253,268,"DD","B",291,261,226,268,265,259,"AA","B",221,200,219,230,230,254,"BB","B",250,212,212,234,224,248,"CC","B",250,212,217,220,245,261,"DD","B",216,198,199,207,235,251,"AA","B",189,195,199,252,217,258,"BB","B",202,202,202,216,226,263,"CC","B",219,238,178,231,194,247,"DD","B",237,234,201,261,226,265,"AA","B",237,231,223,257,228,260,"BB","B",234,222,197,213,222,240,"CC","B",181,196,186,204,210,236,"DD","B",215,220,202,231,229,273),nrow = 8, ncol=24)));
colnames(d)<-c("Factor1", "Factor2", "a2011", "a2012", "a2013", "a2014", "a2015", "a2016");
d[,1]<-as.factor(d[,1]);
d[,2]<-as.factor(d[,2])
Код для ANOVA:
require(car);
options(contrasts=c('contr.sum','contr.poly'));
m=lm(cbind(a2011,a2012,a2013,a2014,a2015, a2016)~Factor1+Factor2,data=d);
niveauAnnee=factor(1:6);
niveauAnnee=as.factor(c("a2011","a2012", "a2013","a2014","a2015","a2016"));
idata=data.frame(niveauAnnee);
anov<-Anova(m,idata=idata,idesign=~niveauAnnee,type=3);
summary(anov);