Omega kuadrat untuk mengukur efek dalam R?

Buku statistik yang saya baca merekomendasikan omega kuadrat untuk mengukur efek percobaan saya. Saya telah membuktikan menggunakan desain plot terpisah (campuran antara subyek dan desain antara subyek) bahwa faktor dalam subyek saya signifikan secara statistik dengan p <0,001 dan F = 17.

Sekarang saya ingin melihat seberapa besar perbedaannya ... apakah ada implementasi omega kuadrat di suatu tempat untuk R (atau python? Saya tahu ... orang dapat bermimpi;) Mencari di internet untuk hal-hal yang berkaitan dengan R adalah sakit * , saya tidak tahu bagaimana saya bisa menemukan barang dengan C.

Terima kasih!

Fungsi untuk menghitung omega kuadrat mudah untuk ditulis. Fungsi ini mengambil objek yang dikembalikan oleh tes aov, dan menghitung serta mengembalikan dan omega kuadrat:

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)
}

edit: fungsi yang diperbarui untuk model aov n-way:

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)
}

Saya baru-baru ini harus melaporkan .$\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
}

Ini adalah fungsi yang berantakan yang dapat dengan mudah dibersihkan. Ini menghitung parsial , dan mungkin seharusnya hanya digunakan pada desain faktorial antar subyek.$\omega^2$

Saya akan menyarankan bahwa eta square umum dianggap ( ref , ref ) ukuran ukuran efek yang lebih tepat. Ini termasuk dalam output ANOVA dalam paket ez untuk R.

Saya menemukan fungsi kuadrat omega di profil seseorang. Profil yang mereka sediakan online:

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

Paket Daniel "strengejacke" Lüdecke sjstatstidak dapat melakukan omega-squared, partial-omega-squared dll untuk model ANOVA. Coba lihat.

Berikut ini adalah sketsa yang menunjukkan bahwa:

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)