这是一个蒙特卡洛分析的示例。我们将模拟500组具有随机变异的数据集，并为每组拟合一条剂量反应曲线。该脚本会生成一个表格，其中包含所有500次拟合结果的logEC50和HillSlope值。 在运行此脚本之前，您必须创建一个文件，其中第一个数据表的格式应为单分组变量数据集（无 X、Y 子列）。第一个结果表必须是模拟带误差的数据。第二个结果表必须是非线性回归，用于对数据拟合剂量反应曲线。
脚本命令 |
注释 |
Table Prism 1 Clear |
指定文件中的第一个数据表用于存储本脚本的结果。现在清空该表的所有内容（clear）。请注意，本脚本要求第一个数据表格式为单分组变量数据集，且不包含 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斜率值。
