Simbiology createSimFunction data access

7 Ansichten (letzte 30 Tage)
Mohammad Jafarnejad
Mohammad Jafarnejad am 21 Apr. 2021
Beantwortet: Priya Moorthy am 5 Nov. 2021
I am using createSimFunction to run many scenarios. The results come out in an array of simData. How can I efficiently access one of the output observables (e.g. species of interest) from each of the simData results to plot the ranges or confidence intervals from these many scenarios?

Akzeptierte Antwort

Jeremy Huard
Jeremy Huard am 21 Apr. 2021
Bearbeitet: Jeremy Huard am 21 Apr. 2021
Hi Mohammad,
the method selectbyname allows you to extract the data of a specific state.
Alternatively, you can use the attached function to generate a percentile plot for one of more states of your simData array:
sbioplotscan(simdataObj,"Name",["Receptor","Complex"],"Alpha",0.05)
or starting R2021a
sbioplotscan(simdataObj,Name=["Receptor","Complex"],Alpha=0.05)
For Alpha=0.05, the plot will display the median as well as a an area between the 2.5 and 97.5 percentiles.
Disclaimer: This function is a helper function I wrote for my personal use. It has not been thoroughly tested and is provided as is.
I hope this helps.
Best,
Jérémy
  8 Kommentare
Mohammad Jafarnejad
Mohammad Jafarnejad am 22 Apr. 2021
Thanks Jeremy.
The problem with simData.DataNames is the following: you can have species with the same name in different compartments that will show up as duplicates in the DataNames. Assume you have comp1 and comp2. you can define the drug (D) in compa1 (comp1.D) and in comp2 (comp2.D). Now in the DataNames it shows up as 'D' and 'D'. not able to distinguish unless tracking down the parents. selectbyname works well as it can take 'comp1.D' as the input and distinguish which D we are looking for. I think sbioscanplot.m could also benefit from such a feature. I am getting around that by doing a selectbyname before feeding the data to the sbioscanplot os as not to have duplicates.
Also it would be nice to be able to specifiy DataName for the median line to facilitate legends later.
Thanks again for sharing the function. I have made some modifications to it for my personal use.
Cheers
Mohammad
Jeremy Huard
Jeremy Huard am 22 Apr. 2021
Hi Mohammad,
I see what you mean now. You are right, DataNames are not necessarily unique. However, the DataInfo property contains information about the compartment that contains about the species. In my original code, I was only looking in DataNames. I have now updated the function and uploaded the updated version in my original answer.
Best,
Jérémy

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (1)

Priya Moorthy
Priya Moorthy am 5 Nov. 2021
Mohammad,
The SimBiology Model Analyzer app introduced the percentile plot type in R2021b. You can use this plot type to visualize how summary statistics (percentiles, mean, standard deviation, min/max) for a set of timecourses (simulations from a scan or experimental data from imported datasets) change over time.
This new plot type is documented here:
Best,
Priya

Communitys

Weitere Antworten in  SimBiology Community

Kategorien

Mehr zu Scan Parameter Ranges finden Sie in Help Center und File Exchange

Produkte

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by