R: Likelihood object for normal mixture error distribution

lik_normalmix {ashr} | R Documentation

## Likelihood object for normal mixture error distribution

### Description

Creates a likelihood object for ash for use with normal mixture error distribution

### Usage

    lik_normalmix(pilik, sdlik)

### Arguments

`pilik` | a k vector of mixture proportions (k is the number of mixture components), or an n*k matrix that the j'th row the is mixture proportions for betahat_j
`sdlik` | a k vector of component-wise standard deviations, or an n*k matrix that the j'th row the is component-wise standard deviations for betahat_j

### Examples

       e = rnorm(100,0,0.8)
       e[seq(1,100,by=2)] = rnorm(50,0,1.5) # generate e~0.5*N(0,0.8^2)+0.5*N(0,1.5^2)
       betahat = rnorm(100)+e
       ash(betahat, 1, lik=lik_normalmix(c(0.5,0.5),c(0.8,1.5)))

[Package ashr version 2.2-54 Index]