E' disponibile una funzione eta nel pacchetto labstatR, legato al testo di Iacus e Massarotto, Laboratorio di statica con R (McGraw-Hill, 2007).

eta(x,y)

Dove x è la variabile categoriale (o che deve suddividere i gruppi), e t la variabile metrica. Tale funzione produce, oltre che al risultato di eta quadro, anche il boxplot delle variabili (vedi anche: http://cran.r-project.org/web/packages/labstatR/labstatR.pdf).

Se si fosse interessati ad avere la funzione senza il grafico, scaricare o copiare questo script:

Eta2.R

eta2 <- function(x,y){
 my <- by(y,x,mean)
 ny <- by(y,x,length)
 sm <-sum( (my - mean(y))^2*ny )/length(y)
 sm/mean((y-mean(y))^2)
}