Eta/Eta carré de routines dans la R

Question

En dehors de graphique de l'estimation de la linéarité (regard-à-nuage de points de la méthode), qui est utilisé avant d'appliquer une technique de GLM de la famille, il y a plusieurs façons de faire cette estimation arithmétique (c'est à diresans les graphiques).

Maintenant, je vais me concentrer sur de Fisher eta au carré - rapport de corrélation:arithmétiquement, il est égal au carré de Pearson r (coef.de détermination: r²) si la relation entre les deux variables est linéaire.Par conséquent, vous pouvez comparer les valeurs de eta et r et de faire une évaluation sur le type de relation (linéaire ou non).Il fournit des informations à propos de pour cent de la variance de la variable dépendante expliquée (linéaire ou non) par la variable indépendante.Par conséquent, vous pouvez l'appliquer lors de la linéarité hypothèses ne sont pas respectées.

Simplement déclaré:est-il une routine pour les eta/eta carré dans R?

Answer 1

Je suis toujours très étonné, je dois l'avouer...il n'y a pas de facile et simple pour le calcul de η ou η² dans R...J'ai donc écrit une fonction Page Wikipedia.Va ici:

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)
}

Donc cela donne une autre question, je suis sur le point de publier peu de temps...qu'utilisez-vous pour vérifier la linéarité?Cependant, je ne peux pas calculer les valeurs de p, donc si quelqu'un sait comment faire...s'il vous plaît, laissez-moi savoir!

Answer 2

Après la lecture de cette Question, et en essayant de la fonction dans la réponse, je viens de trouver cette la bibliothèque "sjstats".Il y a un Eta Carré de fonction inclus.C'est peut-être utile pour les futurs demandeurs.