Also BTW (since it's taking a little while creating an account at peakoil.com) I wanted to make a comment on your post over there adding random Gaussian noise and then doing kind of a bootstrap type error estimation. I think that approach will underestimate the size of the error bars because you're not taking account of the structure of the errors. I don't know if the residuals of the various models are just auto-correlated or actually systematic, but they are very clearly not iid random. You might want to plot the autocorrelation r^2 as a function of lag in the residual time series from your model fit. If the residual autocorrelations fall off exponentially with lag, IIRC, very roughly you can consider that you have one independent observation per lifetime of the falloff.
Thanks for the comment. You'right, the noise is clearly not white and should be characterized.

Since, I've performed a more rigourous bootstrapping analysis using the R software (I also put the code):

Bootstrapping Technique Applied to the Hubbert Linearization

I looked at the world production (BP data) but I will post results on the US production probably tonight.