compute.corr: Compute the actual correlation

In MPBA/r-sparcc: R-sparcc

Usage

compute.corr(x, iter = 10, th = 0.1)

Arguments

x
iter
th

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, iter = 10, th = 0.1) 
{
    fracs <- counts2frac(x)
    V <- variation.mat(fracs)
    ll1 <- basis.var(fracs, V)
    ll2 <- cor.from.basis(V, ll1[["Vbase"]])
    excluded <- NULL
    for (i in 1:iter) {
        ll3 <- exclude.pairs(ll2[["Cor.mat"]], ll1[["M"]], th = th, 
            excluded = excluded)
        excluded <- ll3[["excluded"]]
        if (!ll3[["flag"]]) {
            ll1 <- basis.var(fracs, V, M = ll3[["M"]], excluded = excluded)
            ll2 <- cor.from.basis(V, ll1[["Vbase"]])
        }
    }
    return(list(Vbase = ll1[["Vbase"]], Cor.mat = ll2[["Cor.mat"]], 
        Cov.mat = ll2[["Cov.mat"]]))
  }

MPBA/r-sparcc documentation built on May 8, 2019, 3:22 p.m.