Stima Cohen per effetto di scala

Question

dati due vettori:

x <- rnorm(10, 10, 1) 
y <- rnorm(10, 5, 5) 

Come calcolare la d di Cohen per grandezza dell'effetto?

Ad esempio, desidero utilizzare pwr package per stimare la potenza di un t-test con varianze disuguali e richiede Cohen's d.

Answer 1

seguito this link e wikipedia, d di Cohen per un t-test sembra essere:

Dove sigma (denominatore) is:

Quindi, con i tuoi dati:

set.seed(45)      ## be reproducible 
x <- rnorm(10, 10, 1)     
y <- rnorm(10, 5, 5) 

cohens_d <- function(x, y) { 
    lx <- length(x)- 1 
    ly <- length(y)- 1 
    md <- abs(mean(x) - mean(y))  ## mean difference (numerator) 
    csd <- lx * var(x) + ly * var(y) 
    csd <- csd/(lx + ly) 
    csd <- sqrt(csd)      ## common sd computation 

    cd <- md/csd      ## cohen's d 
} 
> res <- cohens_d(x, y) 
> res 
# [1] 0.5199662 

Answer 2

Ci sono diversi pacchetti che forniscono una funzione per il calcolo di d di Cohen. È possibile ad esempio utilizzare la funzione cohensD formano il pacchetto lsr:

library(lsr) 
set.seed(45) 
x <- rnorm(10, 10, 1) 
y <- rnorm(10, 5, 5) 
cohensD(x,y) 
# [1] 0.5199662