Error in Pan-models creation

60 views
Skip to first unread message

pharitporn hianghakun

unread,
Jun 21, 2021, 6:57:50 AM6/21/21
to COBRA Toolbox
Hello, 

 

By following the tutorial instruction both live script and one provided on Creation and simulation of personalized microbiota models through metagenomic data integration — The COBRA Toolbox (opencobra.github.io) , I have problems about creating pan-models from AGORA on different taxonomical levels.  For the solver, I use glpk, LP because some installation problem of IBM. 

Related to adding unmapped species' abundance, the instruction description is not clear since I'm a beginner. Could you give me a step by step for adding the unmapped species?

The screen shown error.jpeg

Kind regards, 

Thiyada

Almut Heinken

unread,
Jun 21, 2021, 1:32:09 PM6/21/21
to cobra-...@googlegroups.com
Hi,
could you make sure that the COBRA.papers folder is up to date and then try again? It seems this is caused by an outdated version of the reaction database. You may need to update COBRA.papers separately in GitHub.
Kind regards,
Almut

--

---
You received this message because you are subscribed to the Google Groups "COBRA Toolbox" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cobra-toolbo...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cobra-toolbox/4d9f069e-3f1f-46b0-a6c0-28c65de96749n%40googlegroups.com.

Almut Heinken

unread,
Jun 21, 2021, 7:12:03 PM6/21/21
to cobra-...@googlegroups.com
Hi,
regarding your second question: I have created a function for this purpose, translateMetagenome2AGORA. It takes as input the path to a file with organism abundances with the taxa as rows and the sample names as columns. The tutorial provides an example for such a file. Unfortunately, the function may not work on every possible file since the nomenclature varies. Hence, you need to inspect the output "unmappedRows" to find the taxa that could not be mapped. You could also map the taxa in your file to the corresponding AGORA (pan-)models manually and then normalize the abundances via the function normalizeCoverage.
I hope that helps.
KInd regards,
Almut

PHARITPORN HIANGHAKUN

unread,
Jun 22, 2021, 11:26:43 AM6/22/21
to cobra-...@googlegroups.com

Hello Almut,

I have checked the COBRA.papers folder. It is up to date. After I tried it, there’s the same error.

 

Here is my screen shot in COBRA.papers.

Kind regards,

Pharitporn

 

 

 

From: Almut Heinken
Sent: 22 มิถุนายน 2564 0:32
To: cobra-...@googlegroups.com
Subject: [SOCIAL NETWORK] Re: Error in Pan-models creation

 

Hi,

could you make sure that the COBRA.papers folder is up to date and then try again? It seems this is caused by an outdated version of the reaction database. You may need to update COBRA.papers separately in GitHub.

Kind regards,

Almut

 

On Mon, Jun 21, 2021 at 11:58 AM pharitporn hianghakun <pharitp...@mail.kmutt.ac.th> wrote:

Hello, 

 

 

 

By following the tutorial instruction both live script and one provided on Creation and simulation of personalized microbiota models through metagenomic data integration — The COBRA Toolbox (opencobra.github.io) , I have problems about creating pan-models from AGORA on different taxonomical levels.  For the solver, I use glpk, LP because some installation problem of IBM. 

 

Related to adding unmapped species' abundance, the instruction description is not clear since I'm a beginner. Could you give me a step by step for adding the unmapped species?

 

 

Kind regards, 

 

Thiyada

--

---
You received this message because you are subscribed to the Google Groups "COBRA Toolbox" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cobra-toolbo...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cobra-toolbox/4d9f069e-3f1f-46b0-a6c0-28c65de96749n%40googlegroups.com.

--

---
You received this message because you are subscribed to the Google Groups "COBRA Toolbox" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cobra-toolbo...@googlegroups.com.

Almut Heinken

unread,
Jun 22, 2021, 11:34:05 AM6/22/21
to cobra-...@googlegroups.com
Hi Pharitporn,
may you check that the files cobratoolbox/papers/2021_demeter/input/ReactionDatabase.txt and MetaboliteDatabase.txt are up to date?

Kind regards,
Almut

PHARITPORN HIANGHAKUN

unread,
Jun 22, 2021, 11:56:05 AM6/22/21
to cobra-...@googlegroups.com

Hello Almut,

 

It was confirmed that the file is up to date, last modification was on 22nd of June.

Almut Heinken

unread,
Jun 22, 2021, 12:07:54 PM6/22/21
to cobra-...@googlegroups.com
Hi Pharitporn,
in that case, I am not what the issue is. Could you run the script with "Stop on errors" active to find out where exactly it crashes and take a screenshot? You probably need to change line 120 from parfor to for.
best,
Almut

PHARITPORN HIANGHAKUN

unread,
Jun 22, 2021, 12:30:03 PM6/22/21
to cobra-...@googlegroups.com

Hello Almut,

Here is the screenshot.

--

---
You received this message because you are subscribed to a topic in the Google Groups "COBRA Toolbox" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cobra-toolbox/OAEM7iZLrRo/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cobra-toolbo...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cobra-toolbox/CAKy8DryeHabmJ4iGRzoPnw-Fc1N3eXAK7FbWXMJ2s7rBaWjOmA%40mail.gmail.com.

 

Almut Heinken

unread,
Jun 22, 2021, 12:36:42 PM6/22/21
to cobra-...@googlegroups.com
Hi,
could you disable the parfor loop (by changing it to for) so it stops within rebuildModel? Sorry for the trouble.
best,
Almut

Reply all
Reply to author
Forward
0 new messages