生物医薬品デ，タの曲面近似- MATLAB & S金宝appimulink - MathWorks日本

生物医薬品デ，タの曲面近似

曲线拟合工具箱™ソフトウェアには,麻酔薬の相互作用の研究用サンプルデータがいくつか用意されています。曲線フィッタ，アプリを使用して応答曲面をこのデ，タに当てはめ，薬物相互作用の効果を分析できます。応答曲面モデルは，複合薬における薬力学的相互作用の振る舞いを理解するための優れた方法です。このデ，タは，[1]の結果に基づいています。

通常，麻酔はオピオ▪▪▪ドと鎮静催眠薬で構成される2▪▪▪以上の薬物の作用です。この例では，薬物クラスプロトタ。これらの間の相互作用は,この複合薬に対する鎮痛反応と鎮静反応の4つの異なる基準により測定されます。痛覚測定,テタニ,鎮静作用および喉頭鏡検査が,プロポフォールとレミフェンタニルのさまざまな濃度の組み合わせにおける代用薬の効果の4つの基準を構成します。

この複合薬の応答曲面を対話的に作成するには，次のようにします。

[現在のフォルダ]ブラウザを使用してmatlab工具箱\ \ curvefit \ curvefitフォルダを見けて表示します。
ファereplicationルOpioidHypnoticSynergy.txtを右クリックし，[デ、タの、ンポ、ト]を選択します。。
1. [huawei @ huawei @ huawei @ huawei @ huawei @ huawei]タブの[区切り文字]セクションで，[列区切り記号]の値は既定の[タブ]のままにします。
  ►►►►►►►►►►►►►►►异丙酚、Reminfentanil、痛觉计、强直、镇静、Laryingoscopyを確認します。
2. [huawei @ huawei @ huawei @ huawei @ huawei @ huawei]セクションで，[選択のeconp econpンポト]をクリックして[デ、タの、ンポ、ト]を選択し，用量反応デ，タをmatlab^®ワクスペスにンポトします。。
または、デ、タをプログラムによって、ンポ、トすることもできます。次のコドを入力して，ファルから用量反応デタをmatlabワクスペスに読み込みます。
```
数据= importdata(“OpioidHypnoticSynergy.txt”）;阿片催眠协同= array2table(data.data, data，.．.“VariableNames”, data.textdata);
```
応答曲面を作成するには,X入力およびY入力として2つの薬物を選択し,Z出力として4つの効果のうちの1つを選択しなければなりません。変数をワ，クスペ，スに読み込んだ後，曲線フィッタ，アプリを開いて変数を対話的に選択します。また，関数curveFitterを使用している場合は，初期の近似変数を指定できます。
この例では，曲線フィッタ，アプリを開きます。
```
curveFitter
```
曲線フィッタ，アプリ内の[曲線フィッタ]タブで，[デ，タ]セクションの[デ，タの選択]をクリックします。[近似デ，タの選択]ダ，アログボックスで，まず，[xデ.タ]、[yデ.タ]、[zデ.タ]のドロップダウンリストから(OpioidHypnoticSynergy)テ，ブルを選択します。次に，新しいドロップダウンリストから，変数の(异丙酚)、(Remifentanil)、(痛觉计)を選択します。

アプリによって，痛觉计デ，タの新しい応答曲面が作成されます。既定の近似はデ，タ点を通過する内挿曲面です。
次のいずれかの方法で現在の曲面近似のコピ，を作成します。
1. [曲線フィッタ]タブの[ファ电子邮箱ル]セクションで[複製]ボタンを選択する。
2. [近似テ，ブル]ペ电子邮箱ンで近似を右クリックし，["untitled fit 1"を複製]を選択する。
デ，タを近似する独自の方程式を定義します。[曲線フィッタ]タブの[近似タ电子邮箱プ]セクションで，矢印をクリックしてギャラリ，を開きます。ギャラリ，で，[カスタム]グル，プの[カスタム式]をクリックします。
[近似オプション]ペesc escンで，エディットボックスに例として入力されているカスタム式テキストを選択して削除します。
カスタム式のエディットボックスを使用して，matlabコ，ドを入力しモデルを定義できます。モデルを定義する方程式は，デ，タ変数xおよびyと,一連の固定パラメーターまたは推定可能なパラメーターあるいはその両方のパラメーターに依存しなければなりません。
この論文のモデルは次のとおりです。

$E ＝ \frac{E_{马克斯} ． {（ \frac{C_{一个}}{我 C 50_{一个}} + \frac{C_{B}}{我 C 50_{B}} + α ． \frac{C_{一个}}{我 C 50_{一个}} ． \frac{C_{B}}{我 C 50_{B}} ）}^{n}}{1 + {（ \frac{C_{一个}}{我 C 50_{一个}} + \frac{C_{B}}{我 C 50_{B}} + α ． \frac{C_{一个}}{我 C 50_{一个}} ． \frac{C_{B}}{我 C 50_{B}} ）}^{n}}$

ここで，c_一个およびc_Bは薬物の濃度，IC50A, IC50B， αおよびnは推定する係数です。
これをmatlabコ，ドで次のように定義できます。
```
效果= Emax*(CA/IC50A + CB/IC50B +.．.α* (CA / IC50A)。* (CB / IC50B))。^ n。/.．.(1 + (ca / ic50a + cb / ic50b +.．.α* (CA / IC50A)。* (CB / IC50B)) ^ n);
```
近似する変数と推定するパラメ，タ，をアプリに指示するために，変数名CAおよびCBをxおよびyに書き換える必要があります。エディットボックスにカスタム式を入力するときに，xとyを含めなければなりません。効果の出力は正規化されるため，Emax = 1と仮定します。
カスタム式のエディットボックスに次のテキストを入力します。
```
(x/IC50A + y/IC50B + alpha*(x/IC50A).*(y/IC50B))。^ n。/（1+(x/IC50A + y/IC50B + alpha*(x/IC50A).*(y/IC50B))。^n)
```
曲線フィッタ，アプリが，カスタム式モデルを使用して曲面をデ，タに当てはめます。
[近似オプション]ペesc escンで，近似オプションの一部を変更します。[詳細オプション]をクリックしてセクションを展開します。
1. [ロバスト]の値を(政治)に設定します。
2. [係数の制約]テ，ブルで，(α)係数にいて，[開始点]の値を1に,[下限]を5に設定します。
  
  新しいオプションを使用して近似が更新されます。
[結果]ペ@ @ンを見てみます。次の結果のいずれも表示できます。
- モデルの方程式
- 推定係数の値
- 適合度の統計量
残差のプロットを表示し，曲面に対する点の分布を確認します。[曲線フィッタ]タブの[可視化]セクションで[残差プロット]をクリックします。
曲線フィッターアプリのセッションで現在選択されている近似と,開かれているプロットについて,コードを生成します。[曲線フィッタ]タブの[エクスポ，ト]セクションで[エクスポ，ト]をクリックし，[コ，ド生成]を選択します。
曲線フィッターアプリはセッションからコードを生成し,MATLABエディターにファイルを表示します。このファesc escルには現在のセッションで選択した近似と，開いていたプロットが含まれます。
ファ@ @ルをcreateFit.mという既定の名前で保存します。
近似とプロットを再作成するには,コマンドラインから元のデータまたは新しいデータを入力引数に指定してそのファイルを呼び出します。
この場合,元のデータがワークスペースに引き続き表示されるので,元のデータの変数で関数を実行できます。
```
[fitresult,gof] = createFit(阿片催眠协同。异丙酚,.．.OpioidHypnoticSynergy。Remifentanil,.．.OpioidHypnoticSynergy.Algometry)
```
この関数は、セッション内で選択した近似の图ウィンドウを作成します。曲線フィッターアプリで対話的に作成した曲面と残差の両方のプロットがカスタム近似の图に表示されます。
痛觉计の代わりに强直応答の新しい近似を作成します。
```
[fitresult,gof] = createFit(阿片催眠协同。异丙酚,.．.OpioidHypnoticSynergy。Remifentanil,.．.OpioidHypnoticSynergy.Tetany)
```
プロットに新しい応答ラベルが必要な場合は，ファ。生成したコ，ドを開始点として，ニ，ズに合うように曲面近似およびプロットを変更できます。使用できるメソッドの一覧にいては，sfitを参照してください。

同じ問題例をプログラムによって曲面で近似する方法にいては，カスタム式による生物医薬品デ，タの曲面近似を参照してください。

参照

[1]克恩，史蒂文·E，谢国明，茱莉亚·l·怀特，塔玛奇·d·伊根。丙泊酚-瑞芬太尼在志愿者体内药效相互作用的响应面分析麻醉学100年,没有。6(2004年6月1日):1373-81。https://doi.org/10.1097/00000542 - 200406000 - 00007。