Exponenta Banner
exponenta event banner

график

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

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

Синтаксис

fh = plot(predCI)

Описание

пример

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

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

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

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

Примеры

свернуть все

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

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

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

clear all
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 должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

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

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

Совет

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

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

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

свернуть все

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

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

| | |

Введенный в R2017b

Документация SimBiology
Поддержка