Advertisement
Guest User

Untitled

a guest
Apr 20th, 2019
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.00 KB | None | 0 0
  1. import sys
  2. import os
  3.  
  4. #path = "/Users/rmadupuri/Desktop/case_lists_test"
  5. path = "/Users/rmadupuri/GitHub/datahub/public"
  6.  
  7. #dictionary of case_list and the corresponding name
  8. case_names = {
  9. "cases_all.txt": "All samples",
  10. "cases_sequenced.txt": "Samples with mutation data",
  11. "cases_cna.txt": "Samples with CNA data",
  12. "cases_cnaseq.txt": "Samples with mutation and CNA data",
  13. "cases_log2CNA.txt": "Samples with log2 copy-number data",
  14. "cases_methylation_all.txt": "Samples with methylation data",
  15. "cases_methylation_hm27.txt": "Samples with methylation data (HM27)",
  16. "cases_methylation_hm450.txt": "Samples with methylation data (HM450)",
  17. "cases_miRNA.txt": "Samples with microRNA data (microRNA-Seq)",
  18. "cases_mRNA_U133.txt": "Samples with mRNA data (U133 microarray)",
  19. "cases_mRNA.txt": "Samples with mRNA data (Agilent microarray)",
  20. "cases_protein_quantification.txt": "Samples with protein data (Mass Spec)",
  21. "cases_RNA_Seq_mRNA.txt": "Samples with mRNA data (RNA Seq)",
  22. "cases_RNA_Seq_v2_mRNA.txt": "Samples with mRNA data (RNA Seq V2)",
  23. "cases_rppa.txt": "Samples with protein data (RPPA)",
  24. "cases_complete.txt": "Complete samples"
  25. }
  26.  
  27. #dictionary of case_list and their corresponding descriptions
  28. case_desc = {
  29. "cases_all.txt": "All samples",
  30. "cases_sequenced.txt": "Samples with mutation data",
  31. "cases_cna.txt": "Samples with CNA data",
  32. "cases_cnaseq.txt": "Samples with mutation and CNA data",
  33. "cases_rppa.txt": "Samples protein data (RPPA)",
  34. "cases_log2CNA.txt": "Samples with log2 copy-number data",
  35. "cases_methylation_all.txt": "Samples with methylation data",
  36. "cases_methylation_hm27.txt": "Samples with methylation data (HM27)",
  37. "cases_methylation_hm450.txt": "Samples with methylation data (HM450)",
  38. "cases_miRNA.txt": "Samples with microRNA data",
  39. "cases_mRNA_U133.txt": "Samples with mRNA expression data",
  40. "cases_mRNA.txt": "Samples with mRNA expression data",
  41. "cases_protein_quantification.txt": "Samples with protein data (Mass Spec)",
  42. "cases_RNA_Seq_mRNA.txt": "Samples with mRNA expression data",
  43. "cases_RNA_Seq_v2_mRNA.txt": "Samples with mRNA expression data",
  44. "cases_complete.txt": "Samples with mutation, CNA and expression data"
  45. }
  46.  
  47.  
  48. all_studies = os.listdir(path)
  49.  
  50. for st in all_studies:
  51. if st != ".DS_Store":
  52. files = os.listdir(path+"/"+st)
  53. if "case_lists" in files:
  54. case_list_files = os.listdir(path+"/"+st+"/case_lists")
  55. for val in case_names:
  56. if val in case_list_files:
  57. seq_data = str()
  58. with open(path+"/"+st+"/case_lists"+"/"+val,'r') as seqf:
  59. for line in seqf:
  60. if line.startswith("case_list_name:"):
  61. seq_data += "case_list_name: "+ case_names[val] +"\n"
  62. elif line.startswith("case_list_description:"):
  63. description = line.rstrip('\n').split(' ')
  64. seq_data += "case_list_description: "+case_desc[val]+" "+description[len(description)-2]+" "+description[len(description)-1]+"\n"
  65. else:
  66. seq_data += line
  67. os.remove(path+"/"+st+"/case_lists"+"/"+val)
  68. f= open(path+"/"+st+"/case_lists"+"/"+val,"w+")
  69. f.write(seq_data)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement