biash: Use control genes to estimate the unimodal inflation. Follow...

Description Usage Arguments

Description

Use control genes to estimate the unimodal inflation. Follow this by running ash on the non-control genes.

Usage

1
biash(betahat, ctl, sebetahat, df = NULL, stramash_args = list())

Arguments

betahat

a p vector of estimates

ctl

A vector of logicals of length length(betahat). A TRUE at position i indicates that β_i = 0, i.e. that it is a control gene. A FALSE indicates that it is not a control gene.

sebetahat

a p vector of corresponding standard errors

df

appropriate degrees of freedom for (t) distribution of betahat/sebetahat if likelihood = "t".

stramash_args

A list containing additional arguments to pass to stramash.workhorse


dcgerard/stramash documentation built on May 15, 2019, 1:24 a.m.