Iterative reconstruction of high-dimensional Gaussian graphical models based on a new method to estimate partial correlations under constraints


R code and datasets necessary to reproduce the results presented in the paper
Vincent Guillemot and Anne-Laure Boulesteix

IMPORTANT: please do not forget to change the working directory when asked.

The R code fo the reproduction of the results presented in this article is structured as follows: 
   - three main R files for the reproduction of figures 2 to 7, namely figures2and3.R, figures4and5.R, figures6and7.R; 
   - one R file to reproduce the results on the real dataset, realData.R;
   - the data associated to the latter, "segal_et_al_2001_property.dat";
   - a folder "functions" containing intermediary R functions;
   - the R package "pacose" available on the CRAN (http://cran.r-project.org).



