def decode_ideal_spec(spec): specGraph,weighGraph = decode.spectrum_graph(spec) source = min(spec) sink = max(spec) edgeTo, marked,postorder = decode.depthFirstPaths(specGraph,source) paths,peptides = decode.pathTo(specGraph,weighGraph,source,sink,edgeTo) peptide = [] for i in range(len(paths)-1): peptide.append(weighGraph[paths[i]][specGraph[paths[i]].index(paths[i+1])]) return ''.join(peptide)
def decode_ideal_spec(spec): specGraph, weighGraph = decode.spectrum_graph(spec) source = min(spec) sink = max(spec) edgeTo, marked, postorder = decode.depthFirstPaths(specGraph, source) paths, peptides = decode.pathTo(specGraph, weighGraph, source, sink, edgeTo) peptide = [] for i in range(len(paths) - 1): peptide.append(weighGraph[paths[i]][specGraph[paths[i]].index( paths[i + 1])]) return ''.join(peptide)
def ideal_spectrum(peptide): idealPrefixSpec = [] for i in range(1,len(peptide)+1): idealPrefixSpec.append(decode.get_mass(peptide[:i])) return idealPrefixSpec
def ideal_spectrum(peptide): idealPrefixSpec = [] for i in range(1, len(peptide) + 1): idealPrefixSpec.append(decode.get_mass(peptide[:i])) return idealPrefixSpec