sim_stats_mo {DBpower}

Simulate power starting from individual-level data for multiple outcomes setting.

Description

Simulate power starting from individual-level data for multiple outcomes setting.

Usage

sim_stats_mo(B, covY, xMat, gVec, alphaMat, gammaVec, checkpoint = FALSE)

Arguments

B

Number of simulations.

covY

Covariance matrix of outcomes.

xMat

Design matrix of non-genetic covariates, n*p.

gVec

n*1 vector of genotypes.

alphaMat

p*K vector of regression coefficients for xMat.

gammaVec

K*1 vector of regression coefficients for each outcome.

checkpoint

Boolean, if true then print message every 50 simulations.

Value

A list with the elements:

zMat

Matrix of test statistics Z.

zVecGBJ

Check on Z statistics, vector should match first row of zMat.

iMat

Innovated statistics using correlation matrix under the null.

Examples

## Not run: 
covY <- matrix(data=0.3, nrow=10, ncol=10); diag(covY) <- 1
xMat <- cbind(1, rnorm(n = 1000), rbinom(n = 1000, size=1, prob=0.5))
gVec <- rbinom(n= 1000, size = 2, prob=0.3)
alphaMat <-matrix(data = 1, nrow=3, ncol=10)
gammaVec <- rep(0, 10)
sim_stats_mo(B=10000, covY = covY, xMat = xMat, gVec = gVec,
alphaMat = alphaMat, gammaVec = gammaVec)

## End(Not run)