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