You could use a monte-carlo simulation. That is, if you have a distribution over the possible values of parameters for your model, you can randomly select from that distribution for each parameter many, many times, then apply the model to those point values, and aggregate the results.
This is done automatically for you in the program guesstimate.
Can you estimate a 90% CI of candyfloss sold over a month based on that information?
Not without additional assumptions, some of which are "obviously" incorrect. In particular:
Unfortunately, in most cases the product of two distributions is a mess. (If you want a pointer, look here.)
One notable set of exceptions is log-transformations of various distributions. (This is because in logspace the multiplication turns into a convolution, which is often a whole lot easier to calculate.)
For instance: the product of two log-normal distributions is "easy". (Of course, then you need to with with the distributions not straight CIs). Beware correlations however.
I'm confused about what concrete question about candyfloss the example is trying to answer, But my usual heuristic for combining estimates is that in the absence of more information (or more realistically, more desire to investigate), I will assume a uniform distribution over some natural scale. For example 10k-40k is on magnitude, so pretend it's uniform over log(10k) to log(40k). 50-8000 is also on magnitude.
Unfortunately, the product of two log-uniform distributions is not a log-uniform distribution...
Suppose you want to take a guess at the number of candyfloss sold over a month in some area and you would like a 90% confidence interval(CI) intead of a point estimate. You fermi estimates for two central subcomponents are:
(1) a 90% CI of the number of candyfloss a single candyfloss seller sells per month, say 10k-40k.
(2) a 90% CI of the number of candyfloss that a professional candyfloss seller sells during a month, say 50-8000.
Can you estimate a 90% CI of candyfloss sold over a month based on that information? If not, could you if you made some assumptions about the distribution of (1) or (2)(e.g. could you do it if they were uniformly distributed)? Could you use the percentiles of the root or squared fo the extremes(e.g. combining either the 0.25th(5%^2) or the 23th(square root of 5%) percentiles)?
My intuition is that you can't just multiply the extremes of (1) and (2) but I'm not confident in what you need to make event approximately correct claims.
Edited to fix an error in (1)