Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The function or variable 'getSecAbsExcMets' is not recognized #42

Open
yqq-79 opened this issue Sep 4, 2023 · 3 comments
Open

The function or variable 'getSecAbsExcMets' is not recognized #42

yqq-79 opened this issue Sep 4, 2023 · 3 comments

Comments

@yqq-79
Copy link

yqq-79 commented Sep 4, 2023

Hello!
I am following COMETS Documentation to learn the Multispecies Example under the COMETS MATLAB TOOLBOX section. When I use the getSecAbsExcMets function to obtain metabolites secreted, absorbed and exchanged by each organism, I encounter this problem: "The function or variable 'getSecAbsExcMets' is not recognized". The previous commands in the example all work smoothly.

@dukovski
Copy link
Contributor

dukovski commented Sep 5, 2023

Hi,

I'll have to take a look at this, not sure what is going on.
The MATLAB toolbox is rarely used so we focus our efforts on the cometspy python toolbox. Would you want to try this with cometspy instead?

@yqq-79
Copy link
Author

yqq-79 commented Sep 26, 2023

Hello!
I try to obtain metabolites secreted, absorbed and exchanged by each organism with cometspy. Which function do I use to achieve this to relapace the getSecAbsExcMets function?

@dukovski
Copy link
Contributor

Look at the basic run in this protocol: https://github.com/segrelab/COMETS_Protocols/blob/master/COMETS_protocols/COMETS_example_Ecoli_CoreModel_WellMixed/Python/P1_test_tube.ipynb Here first you need to make sure that you are saving the extracellular metabolites with this line: sim_params.set_param('writeMediaLog', True) in cell 6. Then you can get them with a code like this: media = experiment.media.copy()
media = media[media.conc_mmol<900]

fig, ax = plt.subplots()
media.groupby('metabolite').plot(x='cycle', ax =ax, y='conc_mmol')

I hope this helps.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants