Description Usage Arguments Details Value Author(s) Examples
With the Ising prior on structured predictors, this function gets the posterior probability of mixing weight.
1 |
SS |
a scalar value of sufficient statistic for regression coefficient. |
beta |
a (p*1) matrix of regression coefficients. |
alpha |
a scalar value of hyperparameter |
hyperparam |
a two-dimensional vector of hyperparameters |
structure |
a data frame stores the information of structure among predictors. |
edgeind |
a vector stores primary keys of |
j |
an index ranges from 1 to p. This function estimates a posterior probability of a mixing weight corresponding to predictor |
With the Ising prior on structured predictors, the problem is transformed into the realm of empirical Bayes thresholding with Laplace prior by estimating the posterior probability of mixing weight. The posterior probability is used to find the posterior median of a regression coefficient.
Return a scalar value of a posterior probability of mixing weight for predictor.
Vitara Pungpapong, Min Zhang, Dabao Zhang
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | data(simGaussian)
Y<-as.matrix(simGaussian[,1])
X<-as.matrix(simGaussian[,-1])
n<-dim(X)[1]
data(linearrelation)
edgeind<-sort(unique(linearrelation[,1]))
# Obtain initial values from lasso
data(initbetaGaussian)
beta<-as.matrix(initbetaGaussian)
# Initiate all other parameters
alpha<-0.5
sigma<-get.sigma(Y=Y, X=X, beta=beta, alpha=alpha)
hyperparam<-get.ab(beta, linearrelation, edgeind)
# Estimate the posterior probability of first predictor
j<-1
Yres<-Y-X%*%beta+X[,j]*beta[j,1]
sxy<-t(Yres)%*%X[,j]
ssx<-sum(X[,j]^2)
SS<-sqrt(n-1)*sxy/(sigma*ssx)
wpost<-get.wpost(SS=SS, beta=beta, alpha=alpha, hyperparam=hyperparam,
structure=linearrelation, edgeind=edgeind, j=j)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.