1 | twodcor8(x, y)
|
x |
|
y |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ##---- 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)
{
X <- elimna(cbind(x, y))
Z1 <- (X[, 1] - mean(X[, 1]))/sqrt(var(X[, 1]))
Z2 <- (X[, 2] - mean(X[, 2]))/sqrt(var(X[, 2]))
temp <- cor.test(Z1 - Z2, X[, 3])
temp <- temp[3]$p.value
list(p.value = temp)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.