Commits

Casey Dunn committed 50a0850

changed plot to greyscale

Comments (0)

Files changed (1)

 
 contrastcor <- cor.table(contrasts)$r
 
+p <- ncol( contrasts )
+n <- nrow( contrasts )
 
 # Regularization by thresholding
 # Bickel, P. J. & Levina, E. Covariance regularization by thresholding. 
 # Ann. Statist. 36, 2577–2604 (2008). http://dx.doi.org/10.1214/08-AOS600
 
-p <- ncol( contrasts )
-n <- nrow( contrasts )
-
 bickel <- contrastcor * ( abs(contrastcor) > sqrt(log(p)/n) )
 
 # Regularization by convex Minimization
 
 
 # Set up plotting parameters
-mycols <- rainbow(100)
-mycols <- heat.colors(100)
-mycols <- tim.colors(100)
+# Color plots
+# mycols <- tim.colors(100)
+# Greyscale plots
+mycols <- gray(1:100 / 100)
+
 mybreaks <- seq(-1, 1, length.out=101)