diff --git a/4-qtl-integration.R b/4-qtl-integration.R
index f7b73e37d28cb5541ebe0f3aa6cae5547c612c5b..8a906c7050595cf782ca7965b2f1c84e61d23b9a 100644
--- a/4-qtl-integration.R
+++ b/4-qtl-integration.R
@@ -91,7 +91,8 @@
   qtl.hist <- ggplot(all.peak, aes(x=qtl_bp, fill = stage)) +
     geom_histogram(binwidth = 2000000, right = T, origin = 0, alpha = 1) +
     #geom_density(alpha = 0.5) +
-    facet_grid(factor(qtl_level, level = c('phQTL', 'mQTL', 'eQTL')) ~ qtl_chromosome, scales="free_y") +
+    facet_grid(factor(qtl_level, level = c('phQTL', 'mQTL', 'eQTL')) ~ qtl_chromosome, scales="free") +
+    
     presentation +
     scale_fill_manual(values = c('#000000', '#ccbb44', '#228833', '#4477aa', '#cc3311')) +
     theme(legend.position = "top") +
diff --git a/5-create-community-network.R b/5-create-community-network.R
index e0e1ebcec9c5a19ddf03c6b97e95dfb0d554c60d..277d04086204b4417cc908c7e7005e5223cac73f 100644
--- a/5-create-community-network.R
+++ b/5-create-community-network.R
@@ -314,6 +314,7 @@
                          mean = round(apply(degree.table[2:ncol(degree.table)], 1, mean), 2),
                          sd = round(apply(degree.table[2:ncol(degree.table)], 1, sd), 2))
   all.rank <- all.rank[order(all.rank$mean), ]
+  write.csv(degree.table, paste0('networks/', dev.stage, chr, '-ranks-for-all-threshold.csv'))
   write.csv(all.rank, paste0('networks/', dev.stage, chr, '-ranks.csv'))
   
   # visualize the network