Scheffe | R Documentation |
Esegue dei confronti multipli tra medie utilizzando la correzione di Scheffe'.
Scheffe( mx, n, MSerr, w = NULL )
mx |
vettore delle medie. |
n |
vettore delle numerosita' campionarie. Se le numerosita' sono tutte uguali e' possibile
inserire un solo valore, se sono diverse, |
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. |
Restistuisce una lista di tre oggetti:
$medie |
medie inserite in input. |
$prob |
probabilita' aggiustate dei confronti effettuati. |
$w |
matrice dei pesi. |
Massimiliano Pastore
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.
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.