Exponenta Banner
exponenta event banner

график

График результатов доверительного интервала для прогнозов модели

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

Синтаксис

fh = график (predCI)

Описание

пример

fh = plot(predCI) строит график доверительного интервала в результате predCI, a 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

Поддержка