¿Cómo se hace un ANOVA SS de tipo III en R con códigos de contraste?

Question

Por favor, proporcione el código de R que permite realizar un ANOVA entre sujetos con contrastes de -3, -1, 1, 3. Entiendo que hay un debate sobre el tipo de suma de cuadrados (SS) apropiado para dicho análisis. Sin embargo, como el tipo de SS por defecto utilizado en SAS y SPSS (Tipo III) se considera el estándar en mi área. Por lo tanto, me gustaría que los resultados de este análisis coincidieran perfectamente con los generados por esos programas estadísticos. Para ser aceptada una respuesta debe llamar directamente a aov(), pero otras respuestas pueden ser votadas (especialmente si son fáciles de entender/usar).

sample.data <- data.frame(IV=rep(1:4,each=20),DV=rep(c(-3,-3,1,3),each=20)+rnorm(80))

Editar: Tenga en cuenta que el contraste que solicito no es un simple contraste lineal o polinómico, sino que se trata de un contraste derivado de una predicción teórica, es decir, el tipo de contrastes analizados por Rosenthal y Rosnow.

Answer 1

Prueba el comando Anova en la biblioteca de coches. Utilice el argumento type="III", ya que por defecto es del tipo II. Por ejemplo:

library(car)
mod <- lm(conformity ~ fcategory*partner.status, data=Moore, contrasts=list(fcategory=contr.sum, partner.status=contr.sum))
Anova(mod, type="III")

Answer 2

También me autopromociono, escribí una función para exactamente esto: https://github.com/samuelfranssens/type3anova

Instálelo de la siguiente manera:

library(devtools)
install_github(samuelfranssens/type3anova)
library(type3anova)

sample.data <- data.frame(IV=rep(1:4,each=20),DV=rep(c(-3,-3,1,3),each=20)+rnorm(80))

type3anova(lm(DV ~ IV, data = sample.data))

También necesitará tener el car paquete instalado.