def postproc(pcResultFolder, pcimFolder, lgnFilePath) : code = DEFAULT_OUTPUT_FILENAME genes_in_lgn, edges_in_lgn = pl.read_lgn(lgnFilePath) list_of_expanded_network_file = pen.merge_different_pc_run(code, pcResultFolder) genes_in_tiles, blocks = pen.read_expanded_network(list_of_expanded_network_file) # 1st internal validation IV = InternalValidation(genes_in_tiles, blocks, genes_in_lgn, edges_in_lgn) expansion_list = IV.expansion_list()
xml_filename = os.path.basename(xml_filepath)
obs_filename, lgn_filename = px.read_single_xml(xml_filepath)
lgn_filepath = helper.generate_filepath(LGN_FOLDER, lgn_filename)
genes_in_lgn, edges_in_lgn = pl.read_lgn(lgn_filepath)
list_of_expanded_network_file = pen.merge_different_pc_run(xml_filename,
PC_RESULT_FOLDER)
# no expanded network found for this xml file
# maybe the result will be run at a later time
if len(list_of_expanded_network_file) == 0:
continue
genes_in_tiles, blocks = pen.read_expanded_network(list_of_expanded_network_file)
print "Internal Validation for %s" % xml_filename
# 1st internal validation
IV = InternalValidation(genes_in_tiles, blocks, genes_in_lgn, edges_in_lgn)
expansion_list = IV.expansion_list()
expansion_filepath = helper.generate_filepath(PCIM_FOLDER_1, xml_filename)
we.write_expansion_list(expansion_filepath, expansion_list)
# 2nd internal validation
# IV2 = InternalValidationRls(genes_in_tiles, blocks, genes_in_lgn, edges_in_lgn)
# expansion_list_2 = IV2.expansion_list()
# expansion_filepath = helper.generate_filepath(PCIM_FOLDER_2, xml_filename)
