getMeans: get means

In AcceptABEL: A package for controlling false positives rate in GWAS by estimating local inflation factor

Description

get means

Usage

getMeans(data, snp.name, trait.name)

Arguments

data
snp.name
trait.name

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (data, snp.name, trait.name) 
{
    snp.num <- which(snpnames(data) == snp.name)
    gtypes <- as.numeric(data@gtdata[, snp.num])
    m1 <- mean(na.omit(data@phdata[which(gtypes == 0), trait.name]))
    m2 <- mean(na.omit(data@phdata[which(gtypes == 1), trait.name]))
    m3 <- mean(na.omit(data@phdata[which(gtypes == 2), trait.name]))
    c(m1, m2, m3)
  }

AcceptABEL documentation built on May 2, 2019, 5:24 p.m.