Comments (7)
I just fixed GEMs issues. I am using Carveme 1.5.2 so I need to change something in the snakefile when I ran the modelVis #52
But I still have a problem with CompositionVis.
from metagem.
Ok thank you so much! It seems that all issues have been fixed.
from metagem.
Hi @SoleilSu,
Regarding GEMs, were you able to install CarveMe and the academic version of the CPLEX solver? Please check your log files for the CarveMe jobs to see what the errror message there was.
Regarding the plot, I suspect that perhaps the abundance estimation output folders are empy, did you run the abundance calculation jobs?
Best,
Francisco
from metagem.
Thank you for the prompt response! I think Carveme works fine and I am able to get the XML files, and when I ran FBA on the XML files, it looks like there are reactions. Here is my modelVis pdf:
I ran the abundance job and I got a nonempty abundance.stats files:
from metagem.
Great! You have your GEMs, taxonomy assignments, and abundance mappings 👍
Seems like the problem then is coming from the plotting script itself compositionVis.R
metaGEM/workflow/scripts/compositionVis.R
Lines 1 to 24 in 69e0629
Could you please share with me your GTDBTk.stats
and abundance.stats
so that I can try to reproduce your error?
I suspect maybe the two files are not being joined properly for some reason and therefore you get a blank plot.
Best,
Francisco
from metagem.
Yea, I just noticed that the files were not joined as expected because the sampleID
was not appended to the user_genome
in the GTDBTk.stats
like this issue here: #110
But it seems that my bug occurs because the same variable is not accessible in the context of my Snakefile. So I made a change in the following code in CompositionRule
in the Snakefile.
for folder in */;do
samp=$(echo $folder|sed 's|^.*/||');
cat $folder/classify/*summary.tsv| sed 's/orig/o/g' | sed 's/permissive/p/g' | sed 's/strict/s/g' | sed "s/^/$samp./g";
done > GTDBTk.stats
I changed the $samp
to $(basename $folder)
like the following:
for folder in */;do
cat $folder/classify/*summary.tsv| sed 's/orig/o/g' | sed 's/permissive/p/g' | sed 's/strict/s/g' | sed "s/^/$(basename $folder)./g";
done > GTDBTk.stats
In this way, my GTDBTk.stats
has the sampleID
, and my plot looks like this now:
However, I am expecting that the different sampleID
have different y-axis of species similar to the one in the tutorial. So I am not sure if I should edit the compositionVis.R
Thank you!
from metagem.
In line 19 of the script shown in my last comment, play around with the scales
value from free
to free_x
or free_y
, or simply remove this parameter from the facet_wrap()
arguments
from metagem.
Related Issues (20)
- feature: add dead-end counter for assembly evaluation
- Installation error HOT 5
- feature: test & enable bioconda installation on macbook HOT 1
- improvement: use `wc` to count bp instead of summing contig header info for abundance HOT 1
- low number of species after binning for metaGEM model reconstruction HOT 2
- automate the response to y/n HOT 6
- [Question]: Install on mac M2? HOT 5
- Job submission with qsub HOT 2
- [Bug]: 'BiGG_gene' is both an index level and a column label, which is ambiguous. HOT 4
- [Usage]: running workflow on workstation with local flag HOT 17
- Questions about media, gapfilling, and predicting interactions HOT 11
- Getting the following error while running the bash metaGEM.sh -t check
- Getting the following error while running the bash metaGEM.sh -t check HOT 8
- refined_bins output remains empty after successful binRefine step HOT 2
- [Question]: How to define and construct a custom culture medium component that can be recognized by CarveMe? HOT 2
- [Question]:Why, when I use CarveMe for gap-filling, does it show that my custom medium does not exist in the database? HOT 1
- [Question]: I meet some errors when I use CarveMe for gap-filling? HOT 3
- [Question]: How to use the GEM output of CarveMe to generate these two files? HOT 1
- [Bug]: Metawrap Installation failure HOT 5
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from metagem.