wcomp.MiRKAT: Weighted MiRKAT

In UVic-omics/CoDA-wKMR:

Description

Runs MiRKAT function using the specified weights for each taxon.

Runs MiRKAT function using the specified weights for each taxon.

Usage

wcomp.MiRKAT(x, y, w = rep(1, nrow(x)), cov = NULL)

wcomp.MiRKAT(x, y, w = rep(1, nrow(x)), cov = NULL)

Arguments

x	matrix with the composition for each individual
y	response variable: dichotomous or continuous
w	vector of weights for each taxa
cov	covariates to hace into account as a possible confounding factors
x	matrix with the composition for each individual
y	response variable: dichotomous or continuous
w	vector of weights for each taxa
cov	covariates to hace into account as a possible confounding factors

Value

the p-value associated to the compositional data with the corresponding vector of weights.

the p-value associated to the compositional data with the corresponding vector of weights.

Examples

# Matrix of samples
 x <- matrix(runif(100), nrow = 10)
# Response vector
 y <- rep(c(0,1),5)
# Vector of weights
 w <- c(2,rep(1,9))

# Run wcomp.MiRKAT function
 wcomp.MiRKAT(x,y,w)

# Matrix of samples
 x <- matrix(runif(100), nrow = 10)
# Response vector
 y <- rep(c(0,1),5)
# Vector of weights
 w <- c(2,rep(1,9))

# Run wcomp.MiRKAT function
 wcomp.MiRKAT(x,y,w)

UVic-omics/CoDA-wKMR documentation built on Oct. 31, 2019, 12:56 a.m.