žPrognozuj Cohena dla efektu wielkości

Question

podanych dwóch wektorów:

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

Jak obliczyć d Cohena do wielkości efekt?

Na przykład chcę użyć pwr package do oszacowania mocy testu t z nierównymi wariancjami i wymaga to d. Cohena.

Answer 1

Following this link i wikipedia, d Cohena dla testu t wydaje się być:

Gdzie sigma (mianownik) jest:

Tak, z twoimi danymi:

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

Istnieje kilka pakietów zapewniających funkcję obliczania współczynnika Cohena. Można na przykład użyć funkcji cohensD tworzą pakiet lsr:

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