ashr (version 0.9)

ashn: Multi-Mode Adaptive Shrinkage function

Description

This is a wrapper function that takes a grid value of $alpha$ and then consider the model $betahat_j ~ g()$,and eqnbeta_j ~ N(0,sebetahat^2) or student t distribution. $alpha$ should be a grid of mu,this wrapper function would select the best $alpha$ and reports the ash item based on that $alpha$.

Usage

ashn(betahat, sebetahat, mixcompdist = c("uniform", "halfuniform", "normal"), df = NULL, alpha = 4, ncores = FALSE, ...)

Arguments

betahat

a p vector of estimates

sebetahat

a p vector of corresponding standard errors

mixcompdist

distribution of components in mixture ( "uniform","halfuniform" or "normal"), the default value would be "uniform"

df

appropriate degrees of freedom for (t) distribution of betahat/sebetahat, default is NULL(Gaussian)

alpha

Could be a vector of grid values for mu. that this wrapper would select based on likelihood principle. Could also be a positive integer greater or equal to 4, then alpha number of grid values would be generated from [-abs(mean(betahat)),2*abs(mean(betahat)], equally spaced.

ncores

Whether to use parallel computing, defaults to FALSE, user could specify number of cores they would like to use. Further, if user does not specify and length(betahat)>50000, then the function would perform parallel computation using number of CPU cores on the current host.

Value

ashm returns a list of objects ashm returns a list of objects

Details

All other inputs are exactly the same as the main function ash, and would pass to the main function to evaluate the likelihood.

See Also

ash the main function that this wrapper function is calling

Examples

beta = c(rep(0,100),rnorm(100))+0.2
sebetahat = abs(rnorm(200,0,1))
betahat = rnorm(200,beta,sebetahat)
beta.ashn = ashn(betahat, sebetahat,alpha=6)
#beta.ashn4 = ashn(betahat, sebetahat,alpha=6,ncores=4)
print(beta.ashn[[1]])  #best ash object
print(beta.ashn[[2]])  #corresponding mode (0 or some other values)
print(beta.ashn[[3]])  #log-likelihood for all models