这是一个蒙特卡洛分析的示例。我们将模拟 500 组具有随机变异性的数据集，并为每组拟合一条剂量反应曲线。该脚本会生成一个表格，其中包含所有 500 次拟合结果的 logEC50 和 HillSlope 值。在运行此脚本之前，您必须创建一个文件，其中第一个数据表的格式需符合单分组变量数据集的要求（不包含 X、Y 子列）。 第一个结果表必须是模拟带误差的数据。第二个结果表必须是针对数据拟合剂量反应曲线的非线性回归。
脚本命令 |
注释 |
Table Prism 1 Clear |
指定文件中的第一个数据表用于存储本脚本的结果。现在清空该表的所有内容（clear）。请注意，本脚本要求第一个数据表格式为单分组变量数据集，且不包含 X列。 |
ForEach 500 |
循环 500 次。 |
GoTo R 1 |
转到第一个结果工作表，该工作表为模拟结果。 |
Regenerate |
使用新的随机数重新生成模拟。 |
GoTo R 2 |
转到第二个结果表。由于该表包含非线性回归结果，Prism 将在每次循环时自动对数据进行重新拟合。 |
Wtable “logEC50”, 5,1 |
将第 5 行 Y 列的第一个值写入输出表的第一列，并将该列命名为“logEC50”。 |
Wtable “Hill Slope”， 6,1 |
将第 6 行第一个 Y 列的值写入输出表的下一列（即第二列），并将该列命名为“Hill斜率”。 |
Next |
循环回到 ForEach 语句。 |
脚本执行完毕后，表 1 将包含两列和 500 行。每行将包含来自某次模拟实验的 logEC50 和 Hill斜率值。
