SimBiology

Verfolgen

When the number of groups are so high, using plot with sbiopredictionci object, it plots all the groups left to righ in a single column, is it possible to change this output?. Thanks

I don't know if I accidentally changed some plotting configuration, but in the past weeks I have been having an issue with the subplots generated by sbiopredictionci. I can edit the data normally in the subplots and I can delete subplots as I wish, but for some reason I can't move the subplots around in the Figure window. It's like they have been fixed in place somehow. Is this configurable? Thanks!

Hi, I do the simulate with my model, and compared the simulated value with observed data, and I want to adjuste the initial values of the parameters according to the comparison results. But when I compared the simulated value with observed data, I find in the Property Editor, when I select "Time" STYLE, the simulated value was very close to the observed data (figure in the left); when I selece "Percentile" STYLE, imulated value was very different with the observed data (figure in the reight). So I don't know why this is happening? This disturbs me, because when I used the model with the intial values used in the simulation step to fit my observed data, there is an error say "Optimization cannot proceed because the initial value of the negative log-likelihood is NaN. ..." In "Time" STYLE, the initial parameter looks like well, in "Percentile" STYLE, the initial parameter looks bad. How can I adjust the parameters?

I am trying to model model glycolysis, PPP and TCA using simbiology. While I get some output plots using local solvers, when I do it with global solvers options it simply does not show any results. (the solving procedure seems to be working properly) Any ideas?

I'm fitting a model to a dataset. I have two variables I'm fitting simultaneously: the therapy, and the quantity of an endogenous protein. Call the model variables Drug and Prot. The corresponding columns in the dataset woudl be DrugObs and ProtObs. Drug and Prot are of different magnitude. I fit the model: [fitcon, simdat] = sbiofit( m1, gmidata, resmap, estpars, doses, 'UseParallel',true); If I use plot(fitcon) I get a trellis plot of data and mode for each group (that is, each subject). This is great! However, since the magnitudes are different I'd like to plot the PRED and OBS of each variable in a different plot. Also, it's useful to look at these things in both log and linear form. So is there any easy way to say plot(fitcon,'PREDVarToPlot','Drug', 'OBSVarToPlot', 'DrugObs'); % and plot(fitcon,'PREDVarToPlot','Prot', 'OBSVarToPlot', 'ProtObs'); % and also plot(fitcon,'PREDVarToPlot','Drug', 'OBSVarToPlot', 'DrugObs',YScale,'log'); % and plot(fitcon,'PREDVarToPlot','Prot', 'OBSVarToPlot', 'ProtObs',YScale,'log'); % I can do this by going into the plot but its tedious: fh = gcf; fhc = fh.Children nfhc = numel(fhc) for i = 1:nfhc test the specific figure child to see if its an axes or a legend if axes get axes children fhcc = fhc(i).Childen % These are lines fhcc(1).Visible = 'off' % This is one way. Any way to delete the line completely? fhcc(3).Visible = 'off' % fhcc([1 3]) are two of the four plotted variables corresponding to etiher drug, or protein. Could be fhcc([2 4]) elseif legend castr = fhc(i).String castr = castr([2 4]) end Another approach would be to allow the second variable to be plotted using a different (right hand) yscale. I've not tried this in script, but as I understand it I'd have to loop through current axes and create a new Child axes with YAxisLocation = 'right'. Not sure then how I would transfer one set of variables (Drug and DrugObs or Prot and ProtObs) to that new axis. Again, it would be ideal if there was an ability to say: plot(fitcon, 'LeftAxis', 'Drug','RightAxis','Prot'); Are there any easy ways to do this, or must I create a script function? Thanks!