Вычислить параметры некомпартментального анализа (NCA) (требуется панель инструментов статистики и машинного обучения)
ncaparameters = sbionca(data,opt)data. Объект options opt определяет столбцы данных и другие параметры расчета. ncaparameters - таблица вычисленных значений параметров NCA для каждой группы.
Загрузите синтетический набор данных, который содержит измерения концентрации лекарств у четырех человек после внутривенной болюсной дозы.
load data1.matУстанавливают дозу NaN в моменты времени, когда доза не вводилась.
data1.Dose(data1.Dose(:) == 0) = NaN;
Просмотрите данные.
sbiotrellis(data1,'ID','Time','DrugConc','Marker','o','LineStyle','--');
Классифицируйте столбцы данных с помощью объекта параметров NCA.
opt = sbioncaoptions; opt.groupColumnName = 'ID'; opt.concentrationColumnName = 'DrugConc'; opt.timeColumnName = 'Time'; opt.IVDoseColumnName = 'Dose';
Вычислите параметры NCA для каждого отдельного пользователя.
ncaparameters = sbionca(data1,opt);
Отображение первых нескольких столбцов таблицы. Каждая строка ncaparameters таблица представляет отдельный объект (или группу), и в каждом столбце отображается соответствующее значение параметра NCA.
ncaparameters(:,1:15)
ans=4×15 table
ID doseSchedule administrationRoute Lambda_Z R2 adjusted_R2 Num_points AUC_0_last Tlast C_max C_max_Dose T_max MRT T_half AUC_infinity
__ ____________ ___________________ ________ _______ ___________ __________ __________ _____ ______ __________ _____ ______ ______ ____________
1 {'Single'} {'IVBolus'} 0.57893 0.99991 0.9999 11 143.61 48 74.412 1488.2 0 1.5408 1.1973 143.61
2 {'Single'} {'IVBolus'} 0.66798 0.99998 0.99998 11 299.37 48 191.96 1919.6 0 1.3352 1.0377 299.37
3 {'Single'} {'IVBolus'} 0.62124 0.99999 0.99999 11 766.5 48 411.06 1644.2 0 1.4476 1.1157 766.5
4 {'Single'} {'IVBolus'} 0.58011 0.99995 0.99995 11 1301.8 48 648.33 1296.7 0 1.5721 1.1949 1301.8
Можно также указать пользовательский временной диапазон для вычисления T_max и C_max в пределах этого временного диапазона, скажем, от времени = 0 до 20. Это можно сделать, установив C_max_ranges свойство как массив ячеек двухэлементного вектора строки.
opt.C_max_ranges = {[5.5 20]};
ncaparameters2 = sbionca(data1,opt);Функция сообщает значения T_max и C_max в пределах диапазона, добавляя два новых столбца: T_max_5_5__20 и C_max_5_5__20. Обратите внимание, что в именах этих двух столбцов последнему моменту времени предшествуют два последовательных подчеркивания (__).
ncaparameters2.T_max_5_5__20(:)
ans = 4×1
6
6
6
6
ncaparameters2.C_max_5_5__20(:)
ans = 4×1
2.2719
3.0213
10.0233
19.9006
Аналогично, можно задать пользовательский временной диапазон для вычисления частичного значения AUC для каждой группы.
opt.PartialAreas = {[0 20]};
ncaparameters3 = sbionca(data1,opt);
ncaparameters3.AUC_0__20(:)ans = 4×1
103 ×
0.1436
0.2994
0.7665
1.3017
Можно также указать несколько временных диапазонов для C_max_ranges и PartialAreas.
opt.C_max_ranges = {[0 20],[0 10],[0 15]};
opt.PartialAreas = {[0 12],[0 30]};
ncaparameters4 = sbionca(data1,opt);