runcor:

In musto101/wilcox_R: What the package does (short line)

Usage

runcor(x, y, z, fr = 1, corflag = FALSE, corfun = pbcor, plotit = TRUE, rhat = FALSE)

Arguments

x
y
z
fr
corflag
corfun
plotit
rhat

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, y, z, fr = 1, corflag = FALSE, corfun = pbcor, plotit = TRUE, 
    rhat = FALSE) 
{
    temp <- cbind(x, y, z)
    temp <- elimna(temp)
    x <- temp[, 1]
    y <- temp[, 2]
    z <- temp[, 3]
    plotit <- as.logical(plotit)
    rmd <- NA
    if (!corflag) {
        for (i in 1:length(x)) {
            flag <- near(z, z[i], fr)
            if (sum(flag) > 2) 
                rmd[i] <- corfun(x[flag], y[flag])$cor
        }
    }
    if (corflag) {
        for (i in 1:length(x)) {
            flag <- near(z, z[i], fr)
            if (sum(flag) > 2) 
                rmd[i] <- cor(x[flag], y[flag])
        }
    }
    if (plotit) {
        plot(c(max(z), min(z), z), c(1, -1, rmd), xlab = "Modifier", 
            ylab = "Correlation", type = "n")
        sz <- sort(z)
        zorder <- order(z)
        sysm <- rmd[zorder]
        lines(sz, sysm)
    }
    if (!rhat) 
        rmd <- "Done"
    rmd
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.