# How to plot the parallelized example in "Parallelized Morris and Sobol Sensitivity Analysis of an ODE"

**URL:** <https://discourse.julialang.org/t/how-to-plot-the-parallelized-example-in-parallelized-morris-and-sobol-sensitivity-analysis-of-an-ode/135034>\
**Category:** Modelling & Simulations\
**Tags:** globalsensitivityjl\
**Created:** [January 13, 2026, 5:03pm UTC](https://discourse.julialang.org/t/how-to-plot-the-parallelized-example-in-parallelized-morris-and-sobol-sensitivity-analysis-of-an-ode/135034 "2026-01-13T17:03:07Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![kraqbi](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/kraqbi/32/206682_2.png) [@kraqbi](https://discourse.julialang.org/u/kraqbi)\
**Post date:** [January 13, 2026, 5:03pm UTC](https://discourse.julialang.org/t/how-to-plot-the-parallelized-example-in-parallelized-morris-and-sobol-sensitivity-analysis-of-an-ode/135034/1 "2026-01-13T17:03:07Z")

</div>

Hello!  
I am trying GolbalSensitivity.jl for the first time, so please forgive me if this is a rookie mistake.

I have trouble with the example in the doc :

> **[Parallelizing the Global Sensitivity Analysis - Parallelized Morris and Sobol...](https://docs.sciml.ai/GlobalSensitivity/stable/tutorials/parallelized_gsa/#Parallelizing-the-Global-Sensitivity-Analysis)**
>
> Let's run GSA on the Lotka-Volterra model to study the sensitivity of the maximum of predator population and the average prey population. | Documentation for GlobalSensitivity.jl.

I reproduce the example until there :  
sobol\_result = gsa(f1, Sobol(), A, B, batch = true)  
but I do not understand how I can plot the result obtained.

In this case, the sobol\_result.ST returns a 400x4 matrix instead of a 2x4 matrix for the non-parallized case.  
Am I doing something wrong? Is something deprecated? How can I extract the result I want?

Thank you for your help!

---

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [January 13, 2026, 5:54pm UTC](https://discourse.julialang.org/t/how-to-plot-the-parallelized-example-in-parallelized-morris-and-sobol-sensitivity-analysis-of-an-ode/135034/2 "2026-01-13T17:54:12Z")

</div>

The doc looks like it has a bug because it doesn’t pass the output\_func so it doesn’t reduce the solution, it gives back the full 200 times timeseries for 2 species, hence the 400. I’ll fix that.