anova3nest: Anàlisi de la Variància amb models de tres factors...

In alexsanchezpla/AnovaTools: Taules ANOVA amb factors aleatoris

Description

Aquesta funció recalcula els elements de la taula ANOVA en el cas d'un model de tres factors quan algun d'ells és aleatori.

Usage

anova3nest(aovObj, random = NULL)

Arguments

aovObj	Un objecte de la classe aov o lm.
random	Un vector amb el nom del factor o factors aleatoris.

Details

Per tal de contrastar correctament la significació del factor aleatori, una vegada tenim un objecte aov o lm del model amb factors fixos, aquesta funció recalcula els estadístics F i els p-valors corresponents de la taula ANOVA amb les mitjanes de les sumes de quadrats apropiades quan un factor és aleatori.

Value

Un objecte de la classe c("anova", "data.frame").

Note

Les estimacions dels paràmetres associats a factors aleatoris que proporciona el summary no tenen sentit.

Author(s)

Jordi Ocaña [e-mail]

Examples

# Exemple dels "bessons"
CAnaerobia <- c(20.6, 20.62, 20.3, 20.7, 22.1, 22.4, 22.3, 23, 21.4, 23.4, 
              24, 23.2, 24, 23.1, 23.5, 22, 21.6, 22.62, 20.5, 21.7, 22.2, 
              24.4, 23.3, 25, 24.4, 22.4, 21, 21.2, 20, 21.1, 22.1, 22)
zigosi <- gl(2,16,labels=c("mono","di"))
parella <- gl(4,4,length=32)
individu <- gl(2,2,length=32)
# La forma més explícita d'indicar el model lineal en el qual
# "parella" és jeràrquicament dins "zigosi" i "individu" és dins les 
# combinacions de "zigosi" i "parella" és:
bessons.aov <- aov(CAnaerobia ~ zigosi + parella %in% zigosi 
                      + individu %in% (parella %in% zigosi))

# Forma equivalent, més concisa
# bessons.aov <- aov(CAnaerobia ~ zigosi/parella/individu)

anova(bessons.aov)
# La taula ANOVA anterior no és correcta.

# La taula ANOVA correcta donat que el primer factor és fix i els dos darrers
# factors són aleatoris és:
anova3nest(bessons.aov, random = c("parella", "individu"))

# Equivalent a:
anova3nest(bessons.aov, random = c("individu", "parella"))
# l'ordre en el qual enumerem els factors aleatoris no importa

# També equivalent a:
anova3nest(bessons.aov, random = "parella")
# tots els jeràrquicament inferiors d'un aleatori els considera aleatoris

alexsanchezpla/AnovaTools documentation built on May 10, 2019, 8:57 a.m.