Description
Monotonic Association on Zero-Inflated Data.
Description
Methods for calculating and testing the significance of pairwise monotonic association from and based on the work of Pimentel (2009) <doi:10.4135/9781412985291.n2>. Computation of association of vectors from one or multiple sets can be performed in parallel thanks to the packages 'foreach' and 'doMC'.
README.md
mazeinda
The goal of mazeinda is to provide functions to compute the monotonic association of zero-inflated vectors, according to the generalization of the tau-b formula proposed by Pimentel(2009).
Example
Mazeinda provides for the user three functions: for obtaining a matrix of association values, $test \textunderscore associations$ for obtaining the p-values of the output of associate and which gives non-zero values only if the association is significantly different from 0. The vectors to be associated can be given as a pair or as columns in matrices.
library(mazeinda)
library(gamlss.dist)
#> Loading required package: MASS
set.seed(224)
### two vectors as imput
x=abs(rBEZI(50, mu = 0.9, sigma = 1, nu = 0.7))
y=abs(rBEZI(50, mu = 0.9, sigma = 1, nu = 0.5))
associate(x,y)
#> [1] 0.07373333
test_associations(x,y)
#> [1] 0.1119453
combine(x,y)
#> [1] 0
### two matrices as imput
m1=matrix(abs(rBEZI(50, mu = 0.9, sigma = 1, nu = 0.7)),5)
m2=matrix(abs(rBEZI(30, mu = 0.9, sigma = 1, nu = 0.5)),5)
associate(m1,m2)
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> Warning in cor(x, y, method = "kendall"): the standard deviation is zero
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] -0.4444444 0.1259882 -0.0800000 -0.0800000 NA -0.0800000
#> [2,] 0.8819171 0.0800000 -0.0800000 -0.0800000 NA 0.8819171
#> [3,] 0.7777778 -0.0800000 -0.0800000 -0.0800000 NA -0.0800000
#> [4,] 0.0000000 -0.3779645 0.6666667 0.3333333 NA -0.5000000
#> [5,] 0.1054093 -0.3585686 0.9486833 0.1054093 NA -0.1054093
#> [6,] -0.5000000 -0.3779645 0.0000000 -0.5000000 NA 0.0000000
#> [7,] -0.0800000 -0.5714286 0.6299408 -0.0800000 NA -0.0800000
#> [8,] 0.2400000 -0.7559289 0.2400000 -0.4444444 NA 0.2400000
#> [9,] NA NA NA NA NA NA
#> [10,] 0.6299408 0.0800000 -0.0800000 -0.0800000 NA 0.6299408
test_associations(m1,m2)
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 0.29791112 0.77554227 0.6453059 0.6453059 0.5 0.64530594
#> [2,] 0.04594145 0.35469406 0.6453059 0.6453059 0.5 0.04594145
#> [3,] 0.06851329 0.64530594 0.6453059 0.6453059 0.5 0.64530594
#> [4,] 1.00000000 0.42919530 0.1467931 0.4681599 0.5 0.27650048
#> [5,] 0.80054211 0.40538056 0.0229774 0.8005421 0.5 0.80054211
#> [6,] 0.27650048 0.42919530 1.0000000 0.2765005 0.5 1.00000000
#> [7,] 0.64530594 0.21263450 0.1539651 0.6453059 0.5 0.64530594
#> [8,] 0.23210716 0.08711844 0.2321072 0.2979111 0.5 0.23210716
#> [9,] 0.50000000 0.50000000 0.5000000 0.5000000 0.5 0.50000000
#> [10,] 0.15396510 0.35469406 0.6453059 0.6453059 0.5 0.15396510
combine(m1,m2)
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> Warning in cor.test.default(x, y, method = "kendall"): Cannot compute exact
#> p-value with ties
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [2,] 0.8819171 0 0.0000000 0 0 0.8819171
#> [3,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [4,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [5,] 0.0000000 0 0.9486833 0 0 0.0000000
#> [6,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [7,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [8,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [9,] 0.0000000 0 0.0000000 0 0 0.0000000
#> [10,] 0.0000000 0 0.0000000 0 0 0.0000000