R中的Eta/Eta平方例程

Question

除了在應用GLM家族的某些技術之前利用的線性（凝視注視法）圖形估計之外，還有幾種方法可以用算術方法進行這種估計（即沒有圖形）。

現在，我將專注於費舍爾的埃塔平方 - 關聯性比例：（COEF測定：[R）算術，它等於平方皮爾森[R如果兩個變量之間的關係是線性的。因此，您可以比較eta和r的值，並對關係類型（線性與否）進行評估。它提供有關自變量解釋（線性或非線性）變量的變異百分比的信息。因此，您可以在線性假設不符合時應用它。

簡單地說：在R中是否有eta/eta-squared的例程？

Answer 1

我還是很愕然，我必須承認...有計算η或η 在R沒有簡單和直接的方式......所以，我根據Wikipedia page寫了一個函數。這裏有：

eta <- function(x, squared = FALSE, ...) { 
    stopifnot(is.list(x)) 
    ## unlist 
    y <- unlist(x) 
    ## group mean 
    mg <- rapply(x, mean, ...) 
    ## group size 
    ng <- rapply(x, length, ...) 
    ## total mean 
    mtot <- mean(y, ...) 
    ## SSb 
    ssb <- sum(ng * (mg - mtot)^2) 
    ## SSt 
    sst <- sum((y - mtot)^2) 
    # get eta-squared 
    if (squared) { 
     res <- ssb/sst 
    # get eta 
    } else { 
     res <- sqrt(ssb/sst) 
    } 
    return(res) 
} 

所以這產生了另一個問題，我即將發佈......你用什麼來檢查線性？但是，我無法計算p值，所以如果有人知道如何去做，請讓我知道！