To propagate the uncertainty of the parameters through the perfume diffusion model, 60,000 samples were drawn from each presented probability density function using both "crude" Monte Carlo and Sobol sequences as described in the methodology, except for those that are Dirac mixtures-specifically, the odor thresholds for limonene, geraniol, and ethanol. To estimate statistics of the odor values of these components, samples from the concentration are apportioned among each discrete possible value that can be assumed. For example, given that the limonene threshold can be either 0.21 or 1.7, the 60,000 trajectories are effectively doubled to 120,000 by dividing each trajectory by each threshold value. The quantity of samples was predetermined, taking into account the runtime of the perfume diffusion model. Nevertheless, to visually assess the appropriate cutoff value for the number of samples, a convergence verification of summary statistics (mean) was executed for both sampling methods. Upon reaching this number, it is reasonable to infer that the sampling statistics exhibit negligible variation.