Kidney

c(“#bdbdbd”,“#a7a7a7”, “#919191”, “#7b7b7b”,“#656565”, “#4f4f4f”) c(“#7974b4”,“#6c82bf”, “#5e8fc5”, “#539cc8”,“#4fa7c6”, “#53b2c3”,“#5fbbbd”,“#71c4b8”,“#86ccb3”,“#9cd3af”) c(“#9cd3af”, “#86ccb3”, “#71c4b8”, “#5fbbbd”, “#53b2c3”, “#4fa7c6”,“#539cc8”,“#5e8fc5”, “#6c82bf”, “#7974b4”)

#tiff("kidney_vp.tif",height = 12, width = 10, units = 'in', res = 200)

plot(x=kidney_vp$n_students_t_test_difference_ko_wt, y=kidney_vp$n_log_students_t_test_p_value_ko_wt, cex=1.9, xlab=expression(italic(log[2])~Fold~change), ylab=expression(italic(-log[10])~P-value), xlim=c(-10,10), col=dr$col, pch=20,  main = "Kidney")
  

 points(labels_vpk$n_students_t_test_difference_ko_wt, labels_vpk$n_log_students_t_test_p_value_ko_wt, col='black', type='p')
 points(kidney_vpsm$n_students_t_test_difference_ko_wt, kidney_vpsm$n_log_students_t_test_p_value_ko_wt, col="#F9AA75",type= "p", pch=19, cex=1.2)
 points(kidney_vpsp$n_students_t_test_difference_ko_wt, kidney_vpsp$n_log_students_t_test_p_value_ko_wt, col="#7974b4",type= "p", pch=19, cex=1.2)
  addTextLabels(labels_vpk$n_students_t_test_difference_ko_wt, labels_vpk$n_log_students_t_test_p_value_ko_wt, labels_vpk$t_gene_names, cex.label=1.2, col.label='black', col.background=NA, avoidPoints=TRUE)
 lines(kidney_vpc$x, kidney_vpc$y)

 #dev.off()

liver

#tiff("liver_vp.tif",height = 12, width = 10, units = 'in', res = 200)

plot(x=liver_vp$n_students_t_test_difference_ko_wt, y=liver_vp$n_log_students_t_test_p_value_ko_wt, cex=1.9, xlab=expression(italic(log[2])~Fold~change), ylab=expression(italic(-log[10])~P-value), xlim=c(-10,10), col=dr$col, pch=20,  main = "Liver")
  
  
 lines(liver_vpc$x, liver_vpc$y)
 points(labels_vpl$n_students_t_test_difference_ko_wt, labels_vpl$n_log_students_t_test_p_value_ko_wt, col='black', type='p')
 
 points(liver_vpsm$n_students_t_test_difference_ko_wt, liver_vpsm$n_log_students_t_test_p_value_ko_wt, col="#F9AA75",type= "p", pch=19, cex=1.2)
 points(liver_vpsp$n_students_t_test_difference_ko_wt, liver_vpsp$n_log_students_t_test_p_value_ko_wt, col="#7974b4",type= "p", pch=19, cex=1.2)
  
 addTextLabels(labels_vpl$n_students_t_test_difference_ko_wt, labels_vpl$n_log_students_t_test_p_value_ko_wt, labels_vpl$t_gene_names, cex.label=1.2,
               col.label='black', col.background=NULL, avoidPoints=TRUE)

 #dev.off()

spleen

#tiff("spleen_vp.tif",height = 12, width = 10, units = 'in', res = 200)

plot(x=spleen_vp$n_students_t_test_difference_ko_wt, y=spleen_vp$n_log_students_t_test_p_value_ko_wt, cex=1.9, xlab=expression(italic(log[2])~Fold~change), ylab=expression(italic(-log[10])~P-value), xlim=c(-10,10), col=dr$col, pch=20, main = "Spleen")
  

addTextLabels(labels_vps$n_students_t_test_difference_ko_wt, labels_vps$n_log_students_t_test_p_value_ko_wt, labels_vps$t_gene_names, cex.label=1.2,
               col.label='black', col.background=NULL, avoidPoints=TRUE)
 lines(spleen_vpc$x, spleen_vpc$y)
 points(labels_vps$n_students_t_test_difference_ko_wt, labels_vps$n_log_students_t_test_p_value_ko_wt, col='black', type='p')
  

 points(spleen_vpsm$n_students_t_test_difference_ko_wt, spleen_vpsm$n_log_students_t_test_p_value_ko_wt, col="#F9AA75",type= "p", pch=19, cex=1.2)
 points(spleen_vpsp$n_students_t_test_difference_ko_wt, spleen_vpsp$n_log_students_t_test_p_value_ko_wt, col="#7974b4",type= "p", pch=19, cex=1.2)
  
 lines(spleen_vpc$x, spleen_vpc$y)

 #dev.off()