How to plot categorical data to compare with SimBiology simulation result?
2 views (last 30 days)
Show older comments
Fearass Zieneddin on 8 Mar 2022
Answered: Arthur Goldsipe on 11 Mar 2022
Is there a way to extract a simbiology simulation result for one time point and then overaly the result to compare with actual data? I'd like to create something like this plot below. The x-axis in my case would represent different sections of tissue samples such as stomach or small intestine and all measuremets would be compared with simulation result from day= 5.
Arthur Goldsipe on 11 Mar 2022
If I understand your question, you're primarily asking how to extract the simulation results at a specific time (for example, 5 days). It sounds like you already have the simulation results, which you obtained by calling sbiosimulate. I'm assuming your simulation time units are day. However, it's not clear to me whether you have those results as a SimData object or as numeric arrays. If you have the results in SimData you can estimate the results at 5 days using the resample method. If you have the results in numeric arrays, you can estimate the resutls at 5 days using the interp1 function. In both cases, I default to the pchip interpolation method. Here's some sample code.
t = 5;
% Case 1: SimData
simdata = sbiosimulate(model);
simdata1 = resample(simdata, t, 'pchip');
[~,x1] = getdata(simdata1);
% Case 2: Numeric arrays
[tManyTimes,xManyTimes] = sbiosimulate(model);
x2 = interp1(tManyTimes, xManyTimes, 5, 'pchip');
% Now plot the data
plot(t, x1, 'o')
plot(t, x2, 'o')
Find more on Import Data in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!