清空筛选条件
清空筛选条件

A 2 compartment PK model

5 次查看（过去 30 天）
Skye Cameron
Skye Cameron 2021-3-18
编辑： Sheida V 2022-7-13
I have the following equations and I am not 100% sure how to start my PK model in the matlab command window? Thanks in advance for any help or advice.

请先登录，再回答此问题。

回答（1 个）

SaiDileep Kola
SaiDileep Kola 2021-3-25
you can start with the commands described here and example models executed.
  2 个评论
Skye Cameron
Skye Cameron 2021-3-25
Hello, you may be able to help me - I have attached my code below, I am using a matrix in a separate file to use my different data for the code. I have attached firstly my code, then the matrix code in the separate file and then the error message. Are you able to help me solve the error please:)
clear all
close all
%Known parameters
alpha = -4.10e-6; %(1/s)
beta = 1.22e-5; %(1/s)
gamma = -8.1e-6; %(1/s)
k_AK = 5e-7; %(1/s)
%transvascular permeability values of cytotoxic drugs
p_DOX = 3e-6; %transvascular permeability in tumour for DOXIRUBICIN
p_PTX = 7e-9; %transvascular permeability in tumour for PACLITAXEL
p_TMZ = 8e-8; %%transvascular permeability in tumour for TEMOZOLOMIDE
p_BCNU = 7e-7; %transvascular permeability in tumour for CARMUSTINE
%elimination rate values for cytotoxic drugs
ke_DOX = 5.8e-4; %elimination rate for DOXIRUBICIN
ke_PTX = 6.8e-7; %elimination rate for PACLITAXEL
ke_TMZ = 1.06e-4; %elimination rate for TEMOZOLOMIDE
ke_BCNU = 1.1e-4; %elimination rate for CARMUSTINE
%elimination rate for Bevacizumab
ke_BEV = 1.2e-5; %elimination rate for BEVACIZUMAB
%half life
thalf_BEV = 57750; %(s)
thalf_DOX = 1194.8; %(s)
thalf_PTX = 1019917.6; %(s)
thalf_TMZ = 9000; %(s)
thalf_BCNU = 6300; %(s)
%plasma clearance
kclr_BEV = (log(2))/thalf_BEV;
kclr_DOX = (log(2))/thalf_DOX;
kclr_PTX = (log(2))/thalf_PTX;
kclr_TMZ = (log(2))/thalf_TMZ;
kclr_BCNU = (log(2))/thalf_BCNU;
%initial ratio of vessel surface area over tissue volume
sv_bt = 20000; %in brain tumour
sv_nt = 7000; %in normal tissue
%PK model
pkm = PKModelDesign;
pkc1 = addCompartment (pkm, 'Central');
pkc1.DosingType = 'Infusion';
pkc1.EliminationType = 'linear-clearance';
pkc1.HasResponseVariable = true;
pkc2 = addCompartment (pkm, 'Peripheral');
pkc1.DosingType = 'Infusion';
pkc2.EliminationType = 'linear-clearance';
pkc2.HasResponseVariable = true;
model = construct(pkm);
configset = getconfigset(model);
configset.CompileOptions.UnitConversion = true;
%Dose for model
dose = sbiodose('dose', 'TargetName', 'Drug_Central');
dose.StartTime = 0;
dose.Amount = 525;
dose.Rate = 175;
dose.AmountUnits = 'milligram';
dose.TimeUnits = 'hour';
dose.RateUnits = 'milligram/hour';
ResponseMap = {'Drug_Central = CentralConc', 'Drug_Peripheral = PeripheralConc'};
paramsToEstimate = {'log(Central)', 'log(Peripheral)', 'Q12', 'Clearance'};
estimatedParam = estimatedInfo(paramsToEstimate, 'InitialValue', [1 1 1 1])';
pooledFit = sbiofit(model,DrugData,responseMap, estimatedParam,dose, 'Pooled', true);
plot(pooledFit)
%-----------------------------------
%Drug, CentralV, PeripheralV, k12, ClearanceCentral
m = [1.90 4.19 1.2e-5 1.2003e-5; %BEV
1.90 4.19 5.8e-4 5.8014e-4; %DOX
1.90 4.19 1.1e-4 1.1002e-4; %BCNU
1.90 4.19 6.8e-7 6.7961e-7; %PTX
1.90 4.19 7.7e-5 7.7016e-5] %TMZ
>> ModelForDrugDelivery
Attempt to execute SCRIPT DrugData as a function:
/Users/skye/Desktop/MatLab/DrugData.m
Error in ModelForDrugDelivery (line 74)
pooledFit = sbiofit(model,DrugData,responseMap, estimatedParam,dose, 'Pooled', true);
Sheida V
Sheida V 2022-7-13
编辑：Sheida V 2022-7-13
One thing I noticed is that you assigned pkc1.DosingType = 'Infusion' twice. I beleive the second time you meant to define pkc2.DosingType.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Simulate Responses to Biological Variability and Doses 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by