taba.test: Robust Correlation Test

In Taba: Taba Robust Correlations

Description

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

Usage

taba.test(x, y, method = c("taba", "tabarank", "tabwil", "tabwilrank"),
          alternative = c("less", "greater", "two.sided"),
          omega, alpha = 0.05)

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", "tabarank", "tabwil", or "tabwilrank" determining if one wants to calculate Taba linear, Taba rank (monotonic), TabWil, or TabWil rank 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 for Taba and Taba rank, and 0.1 for TabWil and TabWil rank. Range is between 0 and 1.
alpha	Type I error rate. Numeric must be between 0 and 1. Default set to 0.05.

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 tuning constant omega.

Value

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

References

Tabatabai, M., Bailey, S., Bursac, Z. et al. An introduction to new robust linear and monotonic correlation coefficients. BMC Bioinformatics 22, 170 (2021). https://doi.org/10.1186/s12859-021-04098-4

doi: 10.1186/s12859-021-04098-4

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

x = rnorm(10)
y = rnorm(10)
taba.test(x, y)
taba.test(x, y, method = "tabarank", alternative = "less")$p.value
taba.test(x, y, method = "tabwil", omega = .1)

Example output

$correlation
[1] 0.5507491

$t.statistic
[1] 1.866309

$p.value
[1] 0.09896679

[1] 0.9948266
$correlation
[1] 0.3819498

$t.statistic
[1] 1.168943

$p.value
[1] 0.2760759

Taba documentation built on April 1, 2021, 1:06 a.m.