这是蒙特卡洛分析的一个示例。我们将模拟500个随机变化的数据集,并为每个数据集拟合一条剂量反应曲线。该脚本将创建一个包含源自所有500条拟合的logEC50和HillSlope的表格。运行该脚本前,必须创建一份文件,该文件的第一张数据表格式为一个分组变量数据集(无X,无Y子列)。第一张结果工作表必须是模拟有错误数据的结果表。第二张结果工作表必须是将剂量响应曲线拟合到数据的非线性回归。
“脚本”命令 |
备注 |
Table Prism 1 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 Slope。