Skip to content
Snippets Groups Projects
Unverified Commit da87177b authored by St. Elmo's avatar St. Elmo
Browse files

format

parent 9ec83298
No related branches found
No related tags found
No related merge requests found
......@@ -55,9 +55,15 @@ productions = screen_variants(
using CairoMakie
disp_rxns = rand(1:95, 10) # only display 10 random fluxes to save space
barplot(productions[disp_rxns], direction = :x, axis = (yticks=(1:10, reactions(model)[disp_rxns]),
xlabel="Flux [mmol/gDW/h]",
ylabel="Reaction ID"))
barplot(
productions[disp_rxns],
direction = :x,
axis = (
yticks = (1:10, reactions(model)[disp_rxns]),
xlabel = "Flux [mmol/gDW/h]",
ylabel = "Reaction ID",
),
)
# ## Knocking out reaction combinations
#
......@@ -82,8 +88,14 @@ productions = screen_variants(
)
f = Figure()
ax = Axis(f[1,1], ylabel="Reaction KO", xlabel="Reaction KO", yticks=(1:10, reactions(model)[disp_rxns]),
xticks=(1:10, reactions(model)[disp_rxns]), xticklabelrotation =-pi/4)
ax = Axis(
f[1, 1],
ylabel = "Reaction KO",
xlabel = "Reaction KO",
yticks = (1:10, reactions(model)[disp_rxns]),
xticks = (1:10, reactions(model)[disp_rxns]),
xticklabelrotation = -pi / 4,
)
hm = heatmap!(ax, productions[disp_rxns, disp_rxns])
Colorbar(f[1,2], hm, label="Flux [mmol/gDW/h]")
Colorbar(f[1, 2], hm, label = "Flux [mmol/gDW/h]")
f
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment