오메가 R의 효과를 측정하기 위해 제곱?

내가 읽고있는 통계 책은 실험의 효과를 측정하기 위해 오메가 제곱을 권장합니다. 필자는 개체 내 요인이 p <0.001 및 F = 17에서 통계적으로 유의하다는 분할 플롯 설계 (피사체와 개체 간 설계의 혼합)를 사용하여 이미 입증했습니다.

이제 차이가 얼마나 큰지보고 싶습니다 ... R에 대한 어딘가에 오메가 제곱이 있습니까 (또는 파이썬? 알고 있습니다 ... 하나는 꿈을 꿀 수 있습니다.) 인터넷에서 R 관련 물건을 검색하는 것은 고통 * , 나는 C로 물건을 찾는 방법을 모른다.

감사!

오메가 제곱을 계산하는 기능은 간단합니다. 이 함수는 aov 테스트에서 반환 된 객체를 가져 와서 계산하고 반환합니다.

omega_sq <- function(aovm){
    sum_stats <- summary(aovm)[[1]]
    SSm <- sum_stats[["Sum Sq"]][1]
    SSr <- sum_stats[["Sum Sq"]][2]
    DFm <- sum_stats[["Df"]][1]
    MSr <- sum_stats[["Mean Sq"]][2]
    W2 <- (SSm-DFm*MSr)/(SSm+SSr+MSr)
    return(W2)
}

편집 : n-way aov 모델의 기능 업데이트 :

omega_sq <- function(aov_in, neg2zero=T){
    aovtab <- summary(aov_in)[[1]]
    n_terms <- length(aovtab[["Sum Sq"]]) - 1
    output <- rep(-1, n_terms)
    SSr <- aovtab[["Sum Sq"]][n_terms + 1]
    MSr <- aovtab[["Mean Sq"]][n_terms + 1]
    SSt <- sum(aovtab[["Sum Sq"]])
    for(i in 1:n_terms){
        SSm <- aovtab[["Sum Sq"]][i]
        DFm <- aovtab[["Df"]][i]
        output[i] <- (SSm-DFm*MSr)/(SSt+MSr)
        if(neg2zero & output[i] < 0){output[i] <- 0}
    }
    names(output) <- rownames(aovtab)[1:n_terms]

    return(output)
}

최근에 $\omega^2$.

partialOmegas <- function(mod){
    aovMod <- mod
    if(!any(class(aovMod) %in% 'aov')) aovMod <- aov(mod)
    sumAov     <- summary(aovMod)[[1]]
    residRow   <- nrow(sumAov)
    dfError    <- sumAov[residRow,1]
    msError    <- sumAov[residRow,3]
    nTotal     <- nrow(model.frame(aovMod))
    dfEffects  <- sumAov[1:{residRow-1},1]
    ssEffects  <- sumAov[1:{residRow-1},2]
    msEffects  <- sumAov[1:{residRow-1},3]
    partOmegas <- abs((dfEffects*(msEffects-msError)) /
                  (ssEffects + (nTotal -dfEffects)*msError))
    names(partOmegas) <- rownames(sumAov)[1:{residRow-1}]
    partOmegas
}

쉽게 정리할 수있는 지저분한 기능입니다. 그것은 부분을 계산$\omega^2$, 아마도 개체 간 요인 설계에만 사용해야합니다.

일반화 된 에타 스퀘어는 ( ref , ref )보다 적절한 효과 크기 측정치로 간주 됩니다 R 용 ez 패키지 의 ANOVA 출력에 포함됩니다 .

누군가의 .Rprofile에서 오메가 제곱 함수를 찾았습니다.

http://www.estudiosfonicos.cchs.csic.es/metodolo/1/.Rprofile

Daniel "strengejacke"Lüdecke의 패키지 sjstats는 ANOVA 모델에 대해 오메가 제곱, 부분 오메가 제곱 등을 수행 할 수 없습니다. 확인 해봐.

다음은이를 보여주는 비 네트입니다.

https://cran.r-project.org/web/packages/sjstats/vignettes/anova-statistics.html

install.packages("sjstats")
library(sjstats)

mod1 <- aov(y~x, data= d.frame)

anova_stats(mod1)