Eta/Eta-quadrado rotinas em R
-
25-09-2019 - |
Pergunta
Além do gráfico de estimativa de linearidade (olhar-em-gráfico de dispersão método), que é utilizado antes de aplicar alguma técnica de GLM família, existem várias maneiras de se fazer esta estimativa aritmeticamente (i.e.sem gráficos).
Agora, vou me concentrar em Fisher eta ao quadrado - correlação taxa:aritmeticamente, é igual ao quadrado de Pearson r (coeficiente.de determinação: r2) se a relação entre duas variáveis é linear.Assim, você pode comparar os valores de eta e r e fazer uma avaliação sobre o tipo de relação (linear ou não).Ele fornece uma informação sobre a percentagem de variância na variável dependente explicada (linear ou não) pela variável independente.Portanto, você pode aplicá-lo quando linearidade premissas não forem atendidas.
Simplesmente afirmou:existe uma rotina para eta/eta-quadrado em R?
Solução
Eu ainda sou muito atordoado, devo admitir...não há fácil e rapidamente para o cálculo de η ou η2 no R
...Então eu escrevi uma função de acordo com A página da wikipédia.Aqui vai:
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)
}
Então isso gera uma outra pergunta, que eu vou postar em breve...o que você usa para verificar a linearidade?No entanto, eu não posso calcular os valores de p, então se alguém sabe como fazer isso...por favor, me avise!
Outras dicas
Depois de ler esta Pergunta, e tentando a função na resposta, eu só encontrei este o biblioteca "sjstats".Há uma Eta-Quadrado-função incluída.Talvez seja útil para o futuro de asilo.