MyNixOS website logo
Description

Monotonic Association on Zero-Inflated Data.

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'.

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
Metadata

Version

0.0.2

License

Unknown

Platforms (75)

    Darwin
    FreeBSD
    Genode
    GHCJS
    Linux
    MMIXware
    NetBSD
    none
    OpenBSD
    Redox
    Solaris
    WASI
    Windows
Show all
  • aarch64-darwin
  • aarch64-genode
  • aarch64-linux
  • aarch64-netbsd
  • aarch64-none
  • aarch64_be-none
  • arm-none
  • armv5tel-linux
  • armv6l-linux
  • armv6l-netbsd
  • armv6l-none
  • armv7a-darwin
  • armv7a-linux
  • armv7a-netbsd
  • armv7l-linux
  • armv7l-netbsd
  • avr-none
  • i686-cygwin
  • i686-darwin
  • i686-freebsd
  • i686-genode
  • i686-linux
  • i686-netbsd
  • i686-none
  • i686-openbsd
  • i686-windows
  • javascript-ghcjs
  • loongarch64-linux
  • m68k-linux
  • m68k-netbsd
  • m68k-none
  • microblaze-linux
  • microblaze-none
  • microblazeel-linux
  • microblazeel-none
  • mips-linux
  • mips-none
  • mips64-linux
  • mips64-none
  • mips64el-linux
  • mipsel-linux
  • mipsel-netbsd
  • mmix-mmixware
  • msp430-none
  • or1k-none
  • powerpc-netbsd
  • powerpc-none
  • powerpc64-linux
  • powerpc64le-linux
  • powerpcle-none
  • riscv32-linux
  • riscv32-netbsd
  • riscv32-none
  • riscv64-linux
  • riscv64-netbsd
  • riscv64-none
  • rx-none
  • s390-linux
  • s390-none
  • s390x-linux
  • s390x-none
  • vc4-none
  • wasm32-wasi
  • wasm64-wasi
  • x86_64-cygwin
  • x86_64-darwin
  • x86_64-freebsd
  • x86_64-genode
  • x86_64-linux
  • x86_64-netbsd
  • x86_64-none
  • x86_64-openbsd
  • x86_64-redox
  • x86_64-solaris
  • x86_64-windows