SimBiology

关注

I have some questions regarding how sensitivity analysis is performed in SimBiology. What is the specific mathematical equation used for computing the sensitivity coefficient of a parameter? Normally what value range for the parameter (e.g. ± 5% of original value) is used for the computation? Is there any way to tweak this? I cannot seem to find any reference for #2 specifically. Thank you!

Hi, I want to do the sensitivity analysis for the PBPK model. I add a program "calculate sensitivities" , but in "SENSIVITITIES TO COMPUTE" ,I can't select my output. My output is the observation in the model, and the observation is the sum of several terms. How can I use the observation as my output to run the sensitivity analysis? Thank you very much.

Is there any built-in matlab function for global sensitivity analysis (moment based like AMA) ?

How to make the most out of the sensitivity analysis? is it enough to just look at the bargraphs and specifying the one(s) with highest magnitude? Why is it unitless? How to interpret the numbers on the y-axis and relate them to the model? Which type of visualization do I choose: simple plots or bargraphs? is it possible to visualize the data as heatmap?

Hello, The title pretty much covers my question: I'm trying to use the 'Calculate sensitivities' Task in SimBiology desktop on a dummy variable whose value is enforced through a repeatedAssignment rule in my model, but I get the error: 'An input or output cannot be the left-handside variable of a repeated assignment rule' I understand what the error is, but I'm wondering if there are any tricks/strategies around this restriction in the sensitivity calculator? Thanks in advance!

I have a PBPK model with the partition coefficients calculated using a mechanistic model. I have a MATLAB function that defines the partition coefficient as an initial assignment in the SimBiology model. This causes a problem when running a sensitivity analysis: "Complex values cannot be converted to logicals" Is there a way to overcome this or do I have to do a manual sensitivity analysis?

Hi I am trying to replicate the results of Sensitivity Analysis I get from Simbiology GUI in my script file. The code is written below: % set the outputs set(csObj.SensitivityAnalysisOptions, 'Outputs',[s1 s2 s3]); % set the inputs set(csObj.SensitivityAnalysisOptions,'Inputs', [k1 k2 k3]); % Set Normalization set(csObj.SensitivityAnalysisOptions,'Normalization','Full'); % Enable Sensitivity set(csObj.SolverOptions, 'SensitivityAnalysis',true); % accelerate sbioaccelerate(m); % Simulate simDataObj = sbiosimulate(m); % Get the data [T, R, snames, ifacs] = getsensmatrix(simDataObj); The Sensitivity matrix I get ('R') has NaNs and Infs in it, so I convert them to zero and sum across the time axis. R(isnan(R)) = 0; R(isinf(R)) = 0; R = abs(R); R1 = sum(R,1); % Sum across time R1 = squeeze(R1); % size of R1 now is states * parameters And compare the results with the plot obtained from GUI [Value, Index] = sort(R1(1,:),'descend') I see that my results are different from those in the GUI, both in terms of values of Sensitivity and order of sensitive parameter. Can anyone help as to what I am doing wrong here in the normalization of 'R' matrix? Thanks S N