Exponenta Banner
exponenta event 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');

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

Создайте модель с двумя отсеками.

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

Поддержка