Scheffe: Test di Scheffe'

ScheffeR Documentation

Test di Scheffe'

Description

Esegue dei confronti multipli tra medie utilizzando la correzione di Scheffe'.

Usage

Scheffe( mx, n, MSerr, w = NULL )

Arguments

mx

vettore delle medie.

n

vettore delle numerosita' campionarie. Se le numerosita' sono tutte uguali e' possibile inserire un solo valore, se sono diverse, n e' calcolato automaticamente come media armonica delle numerosita'.

MSerr

varianza dell'errore ottenuta con l'ANOVA.

w

matrice dei pesi per i contrasti (opzionale), per default esegue tutti i confronti a coppie. Deve avere tante righe quanti sono i contrasti da eseguire e tante colonne quante sono le medie.

Value

Restistuisce una lista di tre oggetti:

$medie

medie inserite in input.

$prob

probabilita' aggiustate dei confronti effettuati.

$w

matrice dei pesi.

Author(s)

Massimiliano Pastore

References

Scheffe', H. (1953). A method for judging all contrasts in the analysis of variance. Biometrika, 40, 87-104.

Keppel, G. (1991). Design and analysis. Prentice Hall, London.

Examples

medie <- c(26.50,37.75,57.50,61.75) # medie osservate
MSerr <- 150.4583 # ricavata dall'ANOVA
# esegue tutti i confronti a coppie
Scheffe(medie,4,MSerr) 

# esempio con solo due contrasti
W <- matrix(c(1,-1,0,0,1,1,-2,0),ncol=4,byrow=TRUE)
Scheffe(medie,4,MSerr,W)

# esempio con numerosita' diverse
Scheffe(c(26.50,37.75,57.50,61.75),c(7,5,4,4),150.4583)

## 
data(relbambini,package="ADati")
mx <- tapply(relbambini$comp.ling,relbambini$cat.comp,mean)
n <- tapply(relbambini$comp.ling,relbambini$cat.comp,length)

summary(aov(comp.ling~cat.comp,data=relbambini)) # ricavo MSerr

Scheffe(mx,n,18.05)

masspastore/ADati documentation built on April 18, 2024, 5:55 p.m.