Scripts underlying the publication: The Impact Of Flocculation on In Situ and Ex Situ Particle Size Measurements by Laser Diffraction

Scripts corresponding to the Monte Carlo bootstrap analysis of the referred paper. The scripts use particle size distributions (PSDs) measured with a LISST-200X laser diffraction instrument to determine the optimum measurement time. The optimal measurement time is a trade of between accuracy and time investment. In the script, the PSD is read from raw LISST output data making use of the provided functions. Next, it randomly drews a subset of measurements and calculated its D50. The size of the subset ranged from one measurement to all measurements in the entire set. Next, a Monte Carlo bootstrap analysis is performed 1000 times for each subset size to determine the deviation of the subset from the data set mean D50. The minimum and maximum values were taken from each run. The measurement frequency (which could be more than 1 measurement per second) can be used to convert the number of measurements, as calculated by the Monte Carlo bootstrap analysis, to measurement time. By studying the change in maximum deviation from the data set mean when adding more measurement readings (when measurement time increases), one can give an estimate on how many readings (and hence measurement time) were needed to give a representative estimate of the D50 of the sample.<br><br><br><br><br><br>

本脚本集对应所引论文的蒙特卡洛自举分析（Monte Carlo bootstrap analysis）脚本。这些脚本依托LISST-200X激光衍射仪测得的粒径分布（Particle Size Distributions, PSDs）数据，用于确定最优测量时长。最优测量时长是测量精度与时间投入之间的权衡结果。脚本借助配套函数从LISST原始输出数据读取粒径分布数据。随后随机抽取一组测量样本并计算其D50值，样本量范围为1次测量至全量测量。针对每一种样本量，均执行1000次蒙特卡洛自举分析，以计算该样本与数据集平均D50值的偏差。每次运行均记录偏差的最小值与最大值。测量频率（单秒可多次测量）可用于将蒙特卡洛自举分析得到的测量次数转换为实际测量时长。通过考察随测量读数增加（即测量时长延长）时，数据集平均D50值的最大偏差变化情况，即可估算获得具有代表性的样本D50值所需的测量次数（亦即实际测量时长）