Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Yoann DUFRESNE
linked reads molecule ordering
Commits
122475c5
Commit
122475c5
authored
Aug 08, 2019
by
Yoann Dufresne
Browse files
Merge branch 'master' of gitlab.pasteur.fr:ydufresne/10x-deconvolve
parents
4c84978d
80c59a67
Changes
1
Hide whitespace changes
Inline
Side-by-side
deconvolution/evaluate.py
View file @
122475c5
...
...
@@ -405,8 +405,8 @@ def verify_graph_edges(d2_component):
# Construct the molecule(s) that this udg really 'reflects'
# i.e. the udg has a central node and two cliques
# that central node
s
is the result of merging of several molecules
# ideally, only one of th
em
is connected to the cliques
# that central node is the result of merging of several molecules
# ideally, only one of th
ose molecules
is connected to the
molecules of the
cliques
# (there could be more than one though; in that case the udg is 'ambiguous')
# udg_molecules aims to reflect the molecule(s) underlying this udg
udg_molecules
=
set
()
...
...
@@ -431,7 +431,7 @@ def verify_graph_edges(d2_component):
nexts
.
sort
(
key
=
lambda
x
:
x
)
quality
=
sum
([
1.0
/
x
if
mol_idx
+
x
in
nexts
else
0
for
x
in
range
(
1
,
6
)])
/
sum
([
1.0
/
x
for
x
in
range
(
1
,
6
)])
if
quality
>
0.6
:
if
quality
>
0.6
:
# eyeballed but still arbitrary
udg_molecules
.
add
(
mol_idx
)
#print("mol",mol_idx,molecule_idxs,"quality",quality,"nexts",nexts)
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment