主要内容
このページの翻訳は最新ではありません。ここをクリックして,英語の最新版を参照してください。
カスタム式による生物医薬品データの曲面近似
この例では,曲线拟合工具箱™を使用して応答曲面を麻酔薬のデータに当てはめ,薬物相互作用の効果を分析する方法を示します。応答曲面モデルは,複合薬における薬力学的相互作用の振る舞いを理解するための優れた方法です。
このデータは,次の論文の結果に基づいています。谢国平，谢国忠，王志强。阿片-催眠协同作用:志愿者丙泊酚-瑞芬太尼药效学相互作用的响应面分析。麻醉学2004;100: 1373 - 81。
通常,麻酔はオピオイドと鎮静催眠薬で構成される2つ以上の薬物の作用です。この例では,薬物クラスプロトタイプとしてプロポフォールとレミフェンタニルを使用します。これらの間の相互作用は,この複合薬に対する鎮痛反応と鎮静反応の4つの異なる基準により測定されます。痛覚測定,テタニ,鎮静作用および喉頭鏡検査が,プロポフォールとレミフェンタニルのさまざまな濃度の組み合わせにおける代用薬の効果の4つの基準を構成します。
次のコードは,曲线拟合工具箱のメソッドを使用して,曲線近似ツールで対話的に曲面作成を行います(生物医薬品データの曲面近似を参照)。
データの読み込み
ファイルからデータを読み込みます。
数据= importdata (“OpioidHypnoticSynergy.txt”）;异丙酚= data.data (: 1);Remifentanil = data.data (:, 2);痛觉计= data.data (: 3);强直= data.data (: 4);镇静= data.data (: 5);Laryingoscopy = data.data (: 6);
モデル近似タイプの作成
関数fittypeを使用して,論文のモデルを定義できます。ここで,CAおよびCBは薬物の濃度,IC50A、IC50B、αおよびnは推定する係数です。モデル近似タイプを作成します。
英国《金融时报》= fittype ('Emax*(CA/IC50A + CB/IC50B + alpha*(CA/IC50A) *(CB/IC50B))^n /(CA/IC50A + CB/IC50B) *(CA/IC50A))^n + 1)'，．..“独立”,{“CA”，“CB”}，“依赖”，“z”，“问题”，“Emax”）
ft =一般模型:ft(IC50A,IC50B,alpha,n,Emax,CA,CB) = Emax*(CA/IC50A + CB/IC50B + alpha*(CA/IC50A) *(CB/IC50B))^n /(CA/IC50A + CB/IC50B + alpha*(CA/IC50A) *(CB/IC50B))^n + 1)
効果の出力は正規化されるため,Emax = 1と仮定します。
Emax = 1;
近似オプションの設定
ロバスト近似,範囲および開始点の近似オプションを設定します。
Opts = fitoptions(ft);选择。[0, 0， -5， -0];选择。健壮的=“守护神”；选择。StartPoint = [0.0089, 0.706, 1.0, 0.746];
痛覚測定の曲面の近似とプロット
[f, gof] = fit([异丙酚，瑞芬太尼]，alogometry, ft，．..选择,“问题”Emax)
成功，但适合停止了，因为变化在残差小于容忍(TolFun)。
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 4.148 (4.123, 4.173) IC50B = 9.043(8.97, 9.117)α= 8.499 (8.314,8.685)n = 8.294(8.136, 8.452)问题参数:Emax = 1
gof =结构体字段:adrsquare: 0.9991 rmse: 0.0146
plot(f，[异丙酚，瑞芬太尼]，痛觉测定法);
テタニへの曲面近似
同じfittypeを再利用して,テタニの応答曲面を作成します。
[f, gof] = fit([异丙酚，瑞芬太尼]，Tetany, ft, opts，“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 4.544 (4.522, 4.567) IC50B = 21.22(21.04, 21.4)α= 14.94 (14.67,15.21)n = 6.132(6.055, 6.209)问题参数:Emax = 1
gof =结构体字段:adrsquare: 0.9993 rmse: 0.0117
plot(f，[异丙酚，瑞芬太尼]，手足搐动);
鎮静作用への曲面近似
[f, gof] = fit([异丙酚，瑞芬太尼]，镇静，ft, opts，“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 1.843 (1.838, 1.847) IC50B = 13.7(13.67, 13.74)α= 1.986 (1.957,2.015)n = 44.27(42.56, 45.98)问题参数:Emax = 1
gof =结构体字段:adrsquare: 0.9994 rmse: 0.0121
plot(f，[异丙酚，瑞芬太尼]，镇静);
喉頭鏡検査への曲面近似
[f, gof] = fit([异丙酚，瑞芬太尼]，喉镜检查，ft, opts，“问题”Emax)
通用模型:f (CA, CB) = Emax * (CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n / ((CA / IC50A + CB / IC50B +α* (CA / IC50A) * (CB / IC50B)) ^ n + 1)系数(95%置信界限):IC50A = 5.192 (5.177, 5.207) IC50B = 37.77(37.58, 37.97)α= 19.67 (19.48,19.86)n = 37(35.12, 38.87)问题参数:Emax = 1
gof =结构体字段:Sse: 0.1555 rsquare: 0.9982 dfe: 393 adjrsquare: 0.9982 rmse: 0.0199
plot(f，[异丙酚，瑞芬太尼]，喉镜检查);
你也可以从以下列表中选择一个网站:
