2010-06-09 611 views
0

除了在应用GLM家族的某些技术之前利用的线性(凝视注视法)图形估计之外,还有几种方法可以用算术方法进行这种估计(即没有图形)。R中的Eta/Eta平方例程

现在,我将专注于费舍尔的埃塔平方 - 关联性比例:(COEF测定:[R)算术,它等于平方皮尔森[R如果两个变量之间的关系是线性的。因此,您可以比较etar的值,并对关系类型(线性与否)进行评估。它提供有关自变量解释(线性或非线性)变量的变异百分比的信息。因此,您可以在线性假设不符合时应用它。

简单地说:在R中是否有eta/eta-squared的例程?

+0

如果r是你的皮尔逊r r 2 – John 2010-06-09 05:21:38

+0

哦......谢谢!深夜登录... = D – aL3xa 2010-06-09 05:42:49

回答

3

我还是很愕然,我必须承认...有计算η或η 在R没有简单和直接的方式......所以,我根据Wikipedia page写了一个函数。这里有:

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值,所以如果有人知道如何去做,请让我知道!