View source: R/EBelasticNet.GaussianCV.R
EBelasticNet.GaussianCV | R Documentation |
Hyperparameter controls degree of shrinkage, and is obtained via Cross Validation (CV). This program calculates the maximum lambda that allows one non-zero basis; and performs a search down to 0.0001*lambda_max at even steps. (20 steps)
EBelasticNet.GaussianCV(BASIS, Target, nFolds,foldId, Epis = FALSE, verbose = 0)
BASIS |
sample matrix; rows correspond to samples, columns correspond to features |
Target |
Response each individual |
nFolds |
number of n-fold cv |
Epis |
TRUE or FALSE for including two-way interactions |
foldId |
random assign samples to different folds |
verbose |
from 0 to 5; larger verbose displays more messages |
If Epis=TRUE, the program adds two-way interaction K*(K-1)/2 more columns to BASIS
CrossValidation |
col1: hyperparameter; col2: loglikelihood mean; standard ERROR of nfold mean log likelihood |
Lmabda_optimal |
the optimal hyperparameter as computed |
Alpha_optimal |
the optimal hyperparameter as computed |
Anhui Huang; Dept of Electrical and Computer Engineering, Univ of Miami, Coral Gables, FL
Huang, A., Xu, S., and Cai, X. (2013). Empirical Bayesian elastic net for multiple quantitative trait locus mapping. submitted.
library(EBEN)
data(BASIS)
data(y)
#reduce sample size to speed up the running time
n = 50;
k = 100;
BASIS = BASIS[1:n,1:k];
y = y[1:n];
## Not run:
CV = EBelasticNet.GaussianCV(BASIS, y, nFolds = 3,Epis = FALSE)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.