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

Problems with code reproduction in tutorials #186

Open
lishijia13579 opened this issue Jul 16, 2023 · 1 comment
Open

Problems with code reproduction in tutorials #186

lishijia13579 opened this issue Jul 16, 2023 · 1 comment

Comments

@lishijia13579
Copy link

lishijia13579 commented Jul 16, 2023

Hi! I am trying to repeat the code from the BacArena tutorial, but I am getting results that don't match those in the tutorial after I set upt the enviroment. For example, after setting up the arena (with the code unchanged), the output information is displayed:
`
[1] "substances with 6562.5 fmol per gridcell:"
[1] "EX_glc(e)"

[1] "substances with 625000 fmol per gridcell:"
[1] "EX_co2(e)" "EX_h(e)" "EX_h2o(e)" "EX_nh4(e)" "EX_o2(e)" "EX_pi(e)"
`
And the example result in the tutorial is as follows:

`
[1] "substances with 322.5 fmol per gridcell:"
[1] "EX_glc(e)"

[1] "substances with 1000 fmol per gridcell:"
[1] "EX_co2(e)" "EX_h(e)" "EX_h2o(e)" "EX_nh4(e)" "EX_o2(e)" "EX_pi(e)"
`
Obviously the amount of substrates in each gridcell is very different from the example result, what could be the reason for this? Is it possibly due to the version of GLPK (mine is GLPK ver. 5.0)?

@lishijia13579
Copy link
Author

Another question.
When I performed futher operations in the tutorial, I found that I got the growth curve of ecoli_core_model in line with the results in the tutorial. However, the curves for the changes of substrates were completely different, I got a plot that with all horizontal straight lines, and I further observed the result of the change in each substrate using the getSubHist() function, I found the the initial amount of products were 0 and gradually increased during the simulation, while the initial amount of the substrates except glucose was 250000000, so even if the quantity decreased, the change could not be observed due to the order of magnitude difference with products. In addition, the initial amount of glucose was different from all other substrates, which is 2625000, and there is also an order of magnitude difference. What is the reason for these significant differences from the results in the tutorial?
Thank you!
image

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

1 participant