Вычислите неразделенный на отсеки анализ (NCA) параметры (требует Statistics and Machine Learning Toolbox),
ncaparameters = sbionca(data,opt)data. Объект 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);