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