Description Usage Arguments Value Author(s) References Examples
A function to achieve covariance selection by glasso with EM, missing data is allowed
1 | glasso.miss(Y, emIter, rho, ...)
|
Y |
numeric matrix. Columns are for variables and rows are for samples. |
emIter |
numeric value. The maximum number of iteration for EM |
rho |
The l1 norm tuning parameter. |
... |
additional control parameter passed to glasso |
a list of following comopnents.
Y.imputed |
The sample matrix Y with missing data imputed by EM |
bic |
the BIC(Bayesian Information Criterion) of current fit. |
... |
other values returned by glasso |
Shiyuan He
He, S.Y., Wang, X., and Yuan, W.(2012), Discovering Co-movement Structure of Chinese Stock Market by SPACE method with EM
Jerome Friedman, Trevor Hastie and Robert Tibshirani (2007). Sparse inverse covariance estimation with the lasso. Biostatistics 2007.
Meinshausen, N. and Buhlmann, P.(2006) High dimensional graphs and variable selection with the lasso. Annals of Statistics,34, p1436-1462.
Daniela Witten, Jerome Friedman, and Noah Simon (2011). New insights and faster computations for the graphical lasso. To appear in Journal of Computational and Graphical Statistics.
1 2 3 4 5 6 7 | data(finStocksCn) ##data of finance sector of Chinese Stock Market, from January 4th, 2011 to July 26th, 2011
finStocksCn$names ##stock names
y.m<-scale(finStocksCn$returns)
n=nrow(y.m)
p=ncol(y.m)
res=glasso.miss(y.m,rho=0.030,emIter=25,penalize.diagonal=FALSE)
res$bic ##bic returned
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.