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