R/localYule.R

localYule<-
function (X) 
{
    sumrow <- apply(X, 1, sum)
    sumcol <- apply(X, 2, sum)
    n <- sum(X)

    result <- X
    for (i in (1:nrow(X))) {
    for (j in (1:ncol(X))) {
        a <- X[i, j]
        b <- sumrow[i] - a
        c <- sumcol[j] - a
        d <- n - a - b - c
        Q <- (a * d - b * c)/(a * d + b * c)
        result[i,j] <- Q
    }
}
    result
}

Try the RcmdrPlugin.pointG package in your browser

Any scripts or data that you put into this service are public.

RcmdrPlugin.pointG documentation built on May 2, 2019, 3:26 p.m.