icmm (version 1.2)

get.zeta: Local posterior probability estimation

Description

This function estimates the local posterior probability when assuming no prior on structured predictors.

Usage

get.zeta(SS, w, alpha)

Arguments

SS

a scalar value of sufficient statistic for regression coefficient.

w

a scalar value of mixing weight.

alpha

a scalar value of hyperparameter alpha.

Value

Return a scalar value of local posterior probability.

Details

Given all other parameters, this function estimates the local posterior probability or the probability that a regression coefficient is not zero conditional on other parameters. This function is called when assuming no prior on structured predictors.

Examples

# NOT RUN {
data(simGaussian)
Y<-as.matrix(simGaussian[,1])
X<-as.matrix(simGaussian[,-1])
n<-dim(X)[1]
# Obtain initial values from lasso
data(initbetaGaussian)
initbeta<-as.matrix(initbetaGaussian)
# Obtain the final output from ebvs
# }
# NOT RUN {
output<-icmm(Y, X, b0.start=0, b.start=initbeta, family = "gaussian", 
        ising.prior = FALSE, estalpha = FALSE, alpha = 0.5, maxiter = 100)
# }
# NOT RUN {
b0<-output$coef[1]
# }
# NOT RUN {
beta<-matrix(output$coef[-1], ncol=1)
# }
# NOT RUN {
# Get all parameters for function arguments
# }
# NOT RUN {
w<-get.wprior(beta)
# }
# NOT RUN {
alpha<-0.5
# }
# NOT RUN {
sigma<-get.sigma(Y,X,beta,alpha)
# }
# NOT RUN {
# Estimate local posterior probability
# }
# NOT RUN {
j<-1
# }
# NOT RUN {
Yres<-Y-b0-X%*%beta+X[,j]*beta[j,1]
# }
# NOT RUN {
sxy<-t(Yres)%*%X[,j]
# }
# NOT RUN {
ssx<-sum(X[,j]^2)
# }
# NOT RUN {
SS<-sqrt(n-1)*sxy/(sigma*ssx)
# }
# NOT RUN {
zeta<-get.zeta(SS=SS, w=w, alpha=alpha)
# }