Eta / eta-квадратные процедуры в r

Question

Помимо графической оценки линейности (метод поблепользования AT-AT-Scatterplot), который используется перед применением некоторых методики из семейства GLM, существует несколько способов сделать эту оценку арифметически (то есть без графов).

Прямо сейчас я сосредоточусь на Фишере ETA-Squared - Соотношение корреляции: арифметически, это равно в квадрате в квадрате Пирсона р (Coef. Определение: р²) Если отношения между двумя переменными являются линейными. Следовательно, вы можете сравнить значения eta. а также р и сделать оценку о типе отношений (линейных или нет). Он предоставляет информацию о процентах дисперсии в зависимой переменной, объясненной (линейно или нет) независимой переменной. Поэтому вы можете применить его, когда предположения линейности не выполняются.

Просто заявлено: есть ли рутина для ETA / ETA-квадрата в R?

Answer 1

Я все еще довольно ошеломлен, я должен признать ... нет простого и простого способа для расчета η или η² в R... так что я написал функцию в соответствии с Страница Википедии. Отказ Вот оно:

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-значения, поэтому, если кто-то знает, как это сделать ... Пожалуйста, дайте мне знать!

Answer 2

Прочитав этот вопрос и попробовать функцию в ответ, я только что нашел этоБиблиотека "Sjstats"Отказ Существует в квадрате ETA-функции. Может быть, это полезно для будущих искателей.