Taba (version 0.1.0)

taba.test: Robust Correlation Test

Description

Tests the association between two numeric vectors using Taba robust linear or Taba rank (monotonic) correlation.

Usage

taba.test(x, y, method = c("taba", "tabarank"),
          alternative = c("less", "greater", "two.sided"),
          omega = 0.45)

Arguments

x

A numeric vector of length greater than 2 must be same length as y

y

A numeric vector of length greater than 2 must be same length as x

method

A character string of "taba" or "tabarank" determining if one wants to calculate Taba linear or Taba rank (monotonic) correlation, respectively. If no method is specified, the function will output Taba linear correlation.

alternative

Character string specifying the alternative hypothesis must be one of "less" for negative association, "greater" for positive association, or "two.sided" for difference in association. If the alternative is not specified, the function will default to a two sided test.

omega

Numeric allowing the user to alter the tuning constant. If one is not specified, the function will default to 0.45. Range is between 0 and 1.

Value

This function returns the robust linear or monotonic association between two numeric vectors, along with it's respective test statistic, and p-value.

Details

This function tests the association of two non-empty numeric vectors of length greater than two, or two columns of a data frame or matrix composed of more than two numeric elements. Covariates are combined colomn-wise and can be numeric vectors, matricies, or data frames with numeric cells. Each column in the matrix or data frame will be treated as a different covariate, and must have different names. Missing values in either x or y are deleted row-wise. The two sided test with the null hypothesis correlation is equal to zero. The default is a two sided test using Taba Linear correlation, with the tuning constant omega equal to 0.45.

References

The paper is under review for possible publication.

See Also

taba for calculating Taba linear or Taba rank (monotonic) correlations taba.partial for partial and semipartial correlations taba.gpartial for generalized partial correlations taba.matrix for calculating correlation, p-value, and distance matricies

Examples

# NOT RUN {
x = rnorm(10)
y = rnorm(10)
taba.test(x,y)
taba.test(x,y,method = "tabarank", alternative = "less")$p.value
# }