Description Usage Arguments Details Value Author(s) References Examples
A function to fit regularized multivariate regression model using the MAP penalty.
1 |
X.m |
numeric matrix (n by p): columns correspond to predictor variables and rows correspond to samples. Missing values are not allowed. |
Y.m |
numeric matrix (n by q): columns correspond to response variables and rows correspond to samples. Missing values are not allowed. |
lamL1 |
numeric value: l_1 norm penalty parameter. |
lamL2 |
numeric value: l_2 norm penalty parameter. |
phi0 |
numeric matrix (p by q): an initial estimate of the coefficient matrix of the multivariate regression model; default(=NULL): univariate estimates are used as initial estimates. |
C.m |
numeric matrix (p by q): C_m[i,j]=0 means the corresponding coefficient beta[i,j] is set to be zero in the model; C_m[i,j]=1 means the corresponding beta[i,j] is included in the MAP penalty; C_m[i,j]=2 means the corresponding beta[i,j] is not included in the MAP penalty; default(=NULL): C_m[i,j] are all set to be 1. |
remMap
uses a computationally efficient approach for performing multivariate regression
under the high-dimension-low-sample-size setting (Peng and et. al., 2008).
A list with two components
phi |
the estimated coefficient matrix (p by q) of the regularized multivariate regression model. |
rss.v |
a vector of length q recording the RSS values of the q regressions. |
Jie Peng, Pei Wang, Ji Zhu
J. Peng, J. Zhu, A. Bergamaschi, W. Han, D.-Y. Noh, J. R. Pollack, P. Wang, Regularized Multivariate Regression for Identifying Master Predictors with Application to Integrative Genomics Study of Breast Cancer. (http://arxiv.org/abs/0812.3671)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 | ############################################
############# Generate an example data set
############################################
n=50
p=30
q=30
set.seed(1)
## generate X matrix
X.m<-matrix(rnorm(n*p),n,p)
## generate coefficient
coef.m<-matrix(0,p,q)
hub.n=10
hub.index=sample(1:p, hub.n)
for(i in 1:q){
cur=sample(1:3,1)
temp=sample(hub.index, cur)
coef.m[temp,i]<-runif(length(temp), min=2, max=3)
}
## generate responses
E.m<-matrix(rnorm(n*q),n,q)
Y.m<-X.m
##############################################
# 1. ## fit model for one pair of (lamL1, lamL2)
##############################################
try1=remMap(X.m, Y.m,lamL1=10, lamL2=5, phi0=NULL, C.m=NULL)
################################################################################################
# 2. ## Select tuning parameters with BIC:
## ## computationally easy; but the BIC procedure assumes orthogonality of the design matrix
## ## to estimate the degrees of freedom;
## ## thus it tends to select too small models when the actual design matrix (X.m) is far
## ## from orthogonal
################################################################################################
lamL1.v=exp(seq(log(10),log(20), length=3))
lamL2.v=seq(0,5, length=3)
df.m=remMap.df(X.m, Y.m, lamL1.v, lamL2.v, C.m=NULL)
## The estimated degrees of freedom can be used to select the ranges of tuning parameters.
try2=remMap.BIC(X.m, Y.m,lamL1.v, lamL2.v, C.m=NULL)
pick=which.min(as.vector(t(try2$BIC)))
result=try2$phi[[pick]]
FP=sum(result$phi!=0 & coef.m==0) ## number of false positives
FN=sum(result$phi==0 & coef.m!=0) ## number of false negatives
#BIC selected tuning parameters
print(paste("lamL1=", round(result$lam1,3), "; lamL2=", round(result$lam2,3), sep=""))
print(paste("FP=", FP, "; FN=", FN, sep=""))
#################################################################################################
# 3. ## Select tuning parameters with v-fold cross-validation;
## ## computationally demanding;
## ## but cross-validation assumes less assumptions than BIC and thus is recommended unless
## ## computation is a concern;
# ## alos cv based on unshrinked estimator (ols.cv) is recommended over cv based on shrinked
## ## estimator (rss.cv);
## ## the latter tends to select too large models.
################################################################################################
lamL1.v=exp(seq(log(10),log(20), length=3))
lamL2.v=seq(0,5, length=3)
try3=remMap.CV(X=X.m, Y=Y.m,lamL1.v, lamL2.v, C.m=NULL, fold=5, seed=1)
############ use CV based on unshrinked estimator (ols.cv)
pick=which.min(as.vector(try3$ols.cv))
#pick=which.min(as.vector(try3$rss.cv))
lamL1.pick=try3$l.index[1,pick] ##find the optimal (LamL1,LamL2) based on the cv score
lamL2.pick=try3$l.index[2,pick]
##fit the remMap model under the optimal (LamL1,LamL2).
result=remMap(X.m, Y.m,lamL1=lamL1.pick, lamL2=lamL2.pick, phi0=NULL, C.m=NULL)
FP=sum(result$phi!=0 & coef.m==0) ## number of false positives
FN=sum(result$phi==0 & coef.m!=0) ## number of false negatives
##CV (unshrinked) selected tuning parameters
print(paste("lamL1=", round(lamL1.pick,3), "; lamL2=", round(lamL2.pick,3), sep=""))
print(paste("FP=", FP, "; FN=", FN, sep=""))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.