Learn R Programming

ALLSPICER (version 0.1.9)

ALLSPICE_simulation: ALLSPICE_simulation

Description

Simulate data and run ALLSPICE

Usage

ALLSPICE_simulation(n_ind, n_var, c, r, pi, sigma, mle = TRUE, null = TRUE)

Value

A list of two pieces of results: 1) ALLSPICE test results 2) effect size table: true effect size simulated, effect size estimate from linear model, effect size estimated from MLE

Arguments

n_ind

total number of individuals

n_var

total number of variants

c

slope between the two sets of variant effect sizes, only applicable when `null` == TRUE

r

phenotypic correlation between the two phenotypes

pi

probability of variant of having no effect on the phenotype

sigma

variance of the two sets of effect sizes

mle

whether to use MLE of c to compute the test statistic, use true c value if FALSE

null

whether to simulate data under the null hypothesis (no linear relationship) or the alternative hypothesis

Examples

Run this code
ALLSPICE_simulation(n_ind=10000, n_var=100, c=0.6, r=0.5, pi=0.5, sigma=1, mle = TRUE, null=TRUE)

Run the code above in your browser using DataLab