Beispiel #1
0
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)
Beispiel #2
0
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)
Beispiel #3
0
def ideal_spectrum(peptide):
	idealPrefixSpec = []
	for i in range(1,len(peptide)+1):
		idealPrefixSpec.append(decode.get_mass(peptide[:i]))
	return idealPrefixSpec
Beispiel #4
0
def ideal_spectrum(peptide):
    idealPrefixSpec = []
    for i in range(1, len(peptide) + 1):
        idealPrefixSpec.append(decode.get_mass(peptide[:i]))
    return idealPrefixSpec