这是蒙特卡洛分析的示例。我们将模拟 500 个随机变化的数据集,并对每个数据集进行剂量反应曲线拟合。该脚本将创建一个表格,其中包含所有 500 个拟合数据的 logEC50 和 HillSlope。在运行此脚本之前,您必须创建一个文件,该文件的第一个数据表格格式为单分组变异性数据集(无 X、无 Y 子列)。第一个结果表必须是模拟误差数据的结果表。第二个结果表必须是对数据拟合剂量反应曲线的非线性回归。
脚本命令 |
注释 |
表Prism 1 清除 |
指定文件中的第一个数据表包含本脚本的结果。现在删除该表的所有内容(清除)。请注意,此脚本要求第一个数据表的设置数据表格式为一个分组变异性数据,不含 X 列。 |
ForEach 500 |
循环 500 次。 |
转到 R 1 |
转到第一个结果表,即模拟结果表。 |
重新生成 |
使用新的随机数重新生成模拟结果。 |
转到 R 2 |
转到第二个结果表。由于该工作表包含非线性回归结果,因此Prism 会在每个循环中自动重新编辑数据。 |
Wtable "logEC50", 5,1 |
将第 5 行第一列 Y 的值写入输出表的第一列,并将该列标为 "logEC50"。 |
Wtable "Hill 斜率", 6,1 |
将第 6 行第一列 Y 的值写入输出表的下一列(第二列),并在该列上标注 "Hill 斜坡率"。 |
下一步 |
循环回到 ForEach 语句。 |
脚本完成后,表 1 将有两列 500 行。每一行都将有一个模拟实验的 logEC50 和 Hill 斜率。