genSignal: Generate the Vector of SNP Effect Sizes

In andreyshabalin/simPRS: Fast Simulation of Polygenic Risk Scores

Description

Generate a vector of SNP effect sizes. Supports various signal strength distributions.

Usage

genSignal(NSignalSnps, NTotalSNPs, heritability, signalDistr = "Same")

Arguments

NSignalSnps	Number of causal SNPs.
signalDistr	Distribution of effect sizes. Has to be "Same", "Normal", or "Uniform".
NTotalSNPs	Total number of SNPs.
heritability	SNP heritability, a value between 0 and 1.

Details

The function generates a signal effect size vector that can be used in gwasFast.

Value

Return a vector of length NTotalSNPs.

Author(s)

Andrey A Shabalin [e-mail]

Examples

signal = genSignal(
            NSignalSnps = 10,
            NTotalSNPs = 20,
            heritability = 0.1,
            signalDistr = "Same")
            
gwas = gwasFast(signal = signal, N = 1000)

cbind(beta = gwas$beta, pv = gwas$pv)

andreyshabalin/simPRS documentation built on May 21, 2019, 2 p.m.