Exponenta Banner

plot

Постройте результаты доверительного интервала для предсказаний модели

свернуть все на странице

Синтаксис

fh = plot(predCI)

Описание

пример

fh = plot(predCI) доверительный интервал графиков следует из predCI, PredictionConfidenceInterval объект или вектор из объектов.

Графики функций точки данных наблюдения как черные знаки "плюс" и предсказания модели как сплошные линии.

  • Если состояние доверительного интервала (predCI.Status), constrained или not estimable, функция использует второй цвет по умолчанию (красный), чтобы построить доверительные интервалы.

  • В противном случае функция использует первый (синий) цвет по умолчанию и строит доверительные интервалы как заштрихованные области.

Примеры

свернуть все

Открыть скрипт

Загрузка данных

Загрузите выборочные данные, чтобы соответствовать.

load data10_32R.mat
gData = groupedData(data);
gData.Properties.VariableUnits = {'','hour','milligram/liter','milligram/liter'};
sbiotrellis(gData,'ID','Time',{'CentralConc','PeripheralConc'},'Marker','+',...
            'LineStyle','none');

Создайте модель

Создайте модель 2D отсека.

pkmd                 = PKModelDesign;
pkc1                 = addCompartment(pkmd,'Central');
pkc1.DosingType      = 'Infusion';
pkc1.EliminationType = 'linear-clearance';
pkc1.HasResponseVariable = true;
pkc2                 = addCompartment(pkmd,'Peripheral');
model                = construct(pkmd);
configset            = getconfigset(model);
configset.CompileOptions.UnitConversion = true;

Задайте дозирование

Задайте капельное внутривенное введение.

dose             = sbiodose('dose','TargetName','Drug_Central');
dose.StartTime   = 0;
dose.Amount      = 100;
dose.Rate        = 50;
dose.AmountUnits = 'milligram';
dose.TimeUnits   = 'hour';
dose.RateUnits   = 'milligram/hour';

Задайте параметры

Задайте параметры, чтобы оценить.

responseMap = {'Drug_Central = CentralConc','Drug_Peripheral = PeripheralConc'};
paramsToEstimate   = {'log(Central)','log(Peripheral)','Q12','Cl_Central'};
estimatedParam     = estimatedInfo(paramsToEstimate,...
                                   'InitialValue',[1 1 1 1],...
                                   'Bounds',[0.1 3;0.1 10;0 10;0.1 2]);

Подбирайте модель

Выполните необъединенную подгонку, то есть, один набор предполагаемых параметров для каждого пациента.

unpooledFit = sbiofit(model,gData,responseMap,estimatedParam,dose,'Pooled',false);

Вычислите доверительные интервалы для предсказаний модели

Вычислите 95% доверительных интервалов для предсказанных ответов модели в необъединенной подгонке с помощью Гауссова приближения.

ciPredUnpooled = sbiopredictionci(unpooledFit);

Постройте доверительные интервалы. Если состоянием оценки доверительного интервала является constrained или not estimable, функция использует второй (красный) цвет по умолчанию. В противном случае функция использует первый (синий) цвет по умолчанию. Чтобы видеть последовательность цветов, введите get(groot,'defaultAxesColorOrder').

Каждая группа отображена в каждом столбце, слева направо, в том же порядке, что они появляются в GroupNames свойство объекта, который используется, чтобы пометить каждый столбец.

plot(ciPredUnpooled)

График с помощью одного (черного) цвета.

plot(ciPredUnpooled,'Color',[0 0 0])

Входные параметры

свернуть все

Доверительный интервал заканчивается для предсказаний модели в виде PredictionConfidenceInterval возразите или вектор из объектов.

Если существует несколько групп, каждая группа отображена в каждом столбце, слева направо, в том же порядке, что они появляются в GroupNames свойство объекта. Каждая строка представляет ответ модели.

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'Color',[0 0 0] задает, чтобы построить использование только черный цвет.

Красный Зеленый Синий цветной триплет в виде разделенной запятой пары, состоящей из 'Color' и трехэлементный вектор-строка. По умолчанию доверительные интервалы, которые не ограничиваются границами параметра, заданными в исходной подгонке, построены с помощью первого (синего) цвета по умолчанию, и те, которые ограничиваются границами, построены с помощью второго (красного) цвета по умолчанию. Если доверительный интервал не является допускающим оценку, он также построен в красном. Чтобы видеть порядок цвета по умолчанию, введите get(groot,'defaultAxesColorOrder') или смотрите свойство ColorOrder.

Совет

Используйте эту пару "имя-значение", когда это необходимо, чтобы создать графики с одним цветом, например, в целях публикации.

Пример: 'Color',[0 0 0]

Выходные аргументы

свернуть все

Изобразите указатель графика, возвращенного как указатель фигуры.

Смотрите также

| | |

Введенный в R2017b

Документация SimBiology

Поддержка