Eta / eta-квадратные процедуры в r
-
25-09-2019 - |
Вопрос
Помимо графической оценки линейности (метод поблепользования AT-AT-Scatterplot), который используется перед применением некоторых методики из семейства GLM, существует несколько способов сделать эту оценку арифметически (то есть без графов).
Прямо сейчас я сосредоточусь на Фишере ETA-Squared - Соотношение корреляции: арифметически, это равно в квадрате в квадрате Пирсона р (Coef. Определение: р2) Если отношения между двумя переменными являются линейными. Следовательно, вы можете сравнить значения eta. а также р и сделать оценку о типе отношений (линейных или нет). Он предоставляет информацию о процентах дисперсии в зависимой переменной, объясненной (линейно или нет) независимой переменной. Поэтому вы можете применить его, когда предположения линейности не выполняются.
Просто заявлено: есть ли рутина для ETA / ETA-квадрата в R?
Решение
Я все еще довольно ошеломлен, я должен признать ... нет простого и простого способа для расчета η или η2 в 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-значения, поэтому, если кто-то знает, как это сделать ... Пожалуйста, дайте мне знать!
Другие советы
Прочитав этот вопрос и попробовать функцию в ответ, я только что нашел этоБиблиотека "Sjstats"Отказ Существует в квадрате ETA-функции. Может быть, это полезно для будущих искателей.