get.wprior {icmm}

Mixing weight estimation.

Description

Given other parameters, this function estimates a mixing weight from the mode of its full conditional distribution function.

Usage

get.wprior(beta)

Arguments

beta

a (p*1) matrix of regression coefficients.

Details

Given other parameters, this function estimates a mixing weight from the mode of its full conditional distribution function. This function is called when use the independent prior of predictors (no prior on structured predictors).

Value

Return a scalar value of a mixing weight.

Author(s)

Vitara Pungpapong, Min Zhang, Dabao Zhang

Examples

data(simGaussian)
Y<-as.matrix(simGaussian[,1])
X<-as.matrix(simGaussian[,-1])
# Obtain initial values from lasso
data(initbetaGaussian)
beta<-as.matrix(initbetaGaussian)
# Estimate the mixing weight
w<-get.wprior(beta)