Skip to content

Commit 7324bf1

Browse files
committed
feat: ✨ New MPS refinement model
1 parent 8bcdcab commit 7324bf1

File tree

2 files changed

+38
-2
lines changed

2 files changed

+38
-2
lines changed

pangebin/mps.tex

+29-1
Original file line numberDiff line numberDiff line change
@@ -34,4 +34,32 @@ \subsection{The Maximum Plasmidness Score problem \MPS{}}\label{sec:method:mps}
3434
\objlabel{mps:obj:max_plasmidness_score} % chktex 25 % chktex 35
3535
\end{equation}
3636
where \(\displaystyle PlasmidnessScore = \sum_{i \in \Fragments{}} \plm{i}\frag{i}\).
37-
\end{definition}
37+
\end{definition}
38+
39+
\subsubsection{Find a solution of \MPS{} maximizing the coverage score and then the GC score}\label{sec:method:mps:refinement}
40+
41+
The solution set \(S_\MPS{}\) for \MPS{} contains solution with the same value of \(PlasmidnessScore\) and different values for \(CoverageFlow\) and \(GCScore\) because of~\Cref{mgc:cst:coverage_score_near_optimal,mps:cst:gc_probability_score_near_optimal}.
42+
The final refinement stage consists in fixing the value of \(PlasmidnessScore\) and finding the subset of \(S_\MPS{}\) maximizing first \(CoverageFlow\) and then \(GCScore\).
43+
It results in a last MILP formulation \(\MPS{}'\):
44+
45+
\begin{definition}{\MPS{}' MILP variables}{milp:mpsprime:variables}
46+
The varibales of \MPS{} are kept.
47+
\end{definition}
48+
49+
\begin{definition}{\MPS{}' MILP constraints}{milp:mpsprime:constraints}
50+
To the constraints in \Cref{definition:milp:mps:constraints} we add the following:
51+
52+
The plasmidness score must be equal to the plasmidness score \(P\) found in \MPS{}:
53+
\begin{equation}
54+
PlasmidnessScore = P
55+
\cstlabel{mpsprime:cst:fixed_plasmidness_score} % chktex 25
56+
\end{equation}
57+
\end{definition}
58+
59+
\begin{definition}{\MPS{}' MILP objective function}{milp:mpsprime:objective}
60+
The objective function aims to maximize the coverage flow prior to the GC score:
61+
\begin{equation}
62+
\max ~ |\Fragments{}| CoverageFlow + GCScore
63+
\objlabel{mpsprime:obj:max_coverage_flow_then_gc_probability_score} % chktex 25 % chktex 35
64+
\end{equation}
65+
\end{definition}

prelude.tex

+9-1
Original file line numberDiff line numberDiff line change
@@ -18,7 +18,15 @@ \subsection*{Aims}\label{sec:prelude:aims}
1818
\end{itemize}
1919

2020
\begin{newfeatbox}
21-
Done since last times:
21+
\textbf{2025--03--24} Done since last times:
22+
23+
\begin{CheckList}{Task}
24+
\Task{done}{Propose a refinement of the \MPS{} problem (\Cref{sec:method:mps:refinement})}
25+
\end{CheckList}
26+
\end{newfeatbox}
27+
28+
\begin{newfeatbox}
29+
\textbf{2025--03--17} Done since last times:
2230

2331
\begin{CheckList}{Task}
2432
\Task{done}{Model a connected-flow in MILP to avoid iteratively resolving the problem until the solution subgraph has only one connected component. This answers one point of the PlasBin-flow discussion (\Cref{sec:method:mcf})}

0 commit comments

Comments
 (0)