Eta/Eta-quadrado rotinas em R

Question

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: r²) 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?

Answer 1

Eu ainda sou muito atordoado, devo admitir...não há fácil e rapidamente para o cálculo de η ou η² 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!

Answer 2

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.