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