AovSum: Analysis of variance with the contrasts sum (the sum of the...

Description Usage Arguments Value Author(s) See Also Examples

Description

Analysis of variance with the contrasts sum (the sum of the coefficients is 0)
Test for all the coefficients
Handle missing values

Usage

1

Arguments

formula

the formula for the model 'y~x1+x2+x1:x2'

data

a data-frame

na.action

(where relevant) information returned by model.frame on the special handling of NAs.

...

other arguments, cf the function lm

Value

Retourne des objets

Ftest

a table with the F-tests

Ttest

a table with the t-tests

Author(s)

Francois Husson husson@agrocampus-ouest.fr

See Also

aov, lm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## Example two-way anova
data(senso)
res <- AovSum(Score~ Product + Day , data=senso)
res

## Example two-way anova with interaction
data(senso)
res2 <- AovSum(Score~ Product + Day + Product : Day, data=senso)
res2

## Example ancova
data(footsize)
res3 <- AovSum(footsize ~ size + sex + size : sex, data=footsize)
res3

Example output

Ftest
               SS df      MS F value   Pr(>F)    
Product    2.0278  2  1.0139  3.4428  0.04148 *  
Day       25.4725  1 25.4725 86.4959 1.18e-11 ***
Residuals 12.0742 41  0.2945                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Ttest
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  4.590736   0.081464 56.3528  < 2e-16 ***
Product - 1  0.237057   0.120675  1.9644  0.05628 .  
Product - 2 -0.284745   0.114663 -2.4833  0.01720 *  
Product - 3  0.047688   0.118456  0.4026  0.68935    
Day - 1      0.803378   0.086382  9.3003  < 2e-16 ***
Day - 2     -0.803378   0.086382 -9.3003  < 2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Ftest
                 SS df      MS F value    Pr(>F)    
Product      1.7773  2  0.8886  2.9876   0.06206 .  
Day         25.8418  1 25.8418 86.8787 1.804e-11 ***
Product:Day  0.4738  2  0.2369  0.7964   0.45813    
Residuals   11.6004 39  0.2974                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Ttest
                        Estimate Std. Error t value Pr(>|t|)    
(Intercept)            4.5842929  0.0871239 52.6181  < 2e-16 ***
Product - 1            0.2557071  0.1225836  2.0860  0.04357 *  
Product - 2           -0.2509596  0.1203159 -2.0858  0.04358 *  
Product - 3           -0.0047475  0.1266521 -0.0375  0.97029    
Day - 1                0.8120707  0.0871239  9.3209  < 2e-16 ***
Day - 2               -0.8120707  0.0871239 -9.3209  < 2e-16 ***
Product - 1 : Day - 1  0.0279293  0.1225836  0.2278  0.82096    
Product - 2 : Day - 1 -0.1454040  0.1203159 -1.2085  0.23412    
Product - 3 : Day - 1  0.1174747  0.1266521  0.9275  0.35935    
Product - 1 : Day - 2 -0.0279293  0.1225836 -0.2278  0.82096    
Product - 2 : Day - 2  0.1454040  0.1203159  1.2085  0.23412    
Product - 3 : Day - 2 -0.1174747  0.1266521 -0.9275  0.35935    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Ftest
              SS df     MS F value    Pr(>F)    
size      72.322  1 72.322 70.4901 1.353e-12 ***
sex        3.114  1  3.114  3.0353   0.08531 .  
size:sex   3.789  1  3.789  3.6930   0.05820 .  
Residuals 82.079 80  1.026                      
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Ttest
                Estimate Std. Error t value Pr(>|t|)    
(Intercept)     4.973278   4.214939  1.1799  0.24153    
size            0.204172   0.024318  8.3958  < 2e-16 ***
sex - F         7.343305   4.214939  1.7422  0.08531 .  
sex - M        -7.343305   4.214939 -1.7422  0.08531 .  
sex - F : size -0.046733   0.024318 -1.9217  0.05820 .  
sex - M : size  0.046733   0.024318  1.9217  0.05820 .  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

FactoMineR documentation built on Jan. 8, 2021, 2:18 a.m.