The expected rate for 100m miles for Uber is 33.333...
But how confident can we be? To answer that let's compute a poisson confidence interval around that rate, as in https://stats.stackexchange.com/questions/10926/how-to-calcu....
Let's see what a 95% confidence interval for 1 death in 3m miles looks like:
> poisson.test(1,conf.level = 0.95)$conf.int
[1] 0.02531781 5.57164339
attr(,"conf.level")
[1] 0.95
Multiply that by 33.333 to convert to deaths per 100m miles: > 33.333333*0.02531781
[1] 0.843927
> 33.333333*5.57164339
[1] 185.7214
So 95% confidence that the rate per 100m miles is from 0.84 to 185.72. That's pretty wide! And since the lower bound crosses 1.18, the difference is not significant at the .05 level (if we must make that particular comparison).
However, let's look at 90% CI: > poisson.test(1,conf.level = 0.9)$conf.int
[1] 0.05129329 4.74386452
attr(,"conf.level")
[1] 0.9
Which gives a CI of 1.71 to 158.13. So with 90% confidence we can say Uber is less safe than sober drivers. Ok.Now let's look at 93% CI:
> poisson.test(1,conf.level = 0.93)$conf.int
[1] 0.03562718 5.17251332
attr(,"conf.level")
[1] 0.93
That gives a CI of 1.188 to 172.417. The lower bound being just a bit worse than sober drivers.So we can conclude with 93% certainty from this data that Uber is less safe than sober drivers. Probably a LOT less safe. Although the CI is really wide, this is shocking data for Uber, in my opinion.
People with more knowledge about statistics than me might be able to explain why.
But the poster did something different. He took it one level further and attempted to calculate this confidence number for different spans in which the actual "deaths per 100 million miles" number of Uber's current cars would fall into, given an ideal world (from a data perspective) in which they would have driven an infinite amount of miles. But he actually did it the other way round - he modified the confidence variable and calculated the spans, and then he adjusted the confidence until he arrived at a span that would put Uber's cars just on par with human driving in the best case.
The fact that a fatal incident happens that early (at 3 million, and not closer or past the 86 million that a statistical human drives on average until a fatal incident occurs) does not allow us to extrapolate a sound number per 100 million miles, but it tells us something about the probability by which the actual number of fatalities by 100 million miles that we'd get if Uber continued testing just like it did and racked up enough miles (and killed people) for a statistically sound calculation will fall into different margins. Sure, Uber could have been just very, very unlucky - but that's pretty unlikely, and the unlikeliness of Uber's bad luck (and conversely the likeliness of the fact that Uber's tech is just systematically deadly) is precisely what can be calculated with this single incident.
The assumption of using Poisson seems pretty solid to me, given we are talking about x events in some continuum (miles traveled in this case), but always happy to hear any cogent objections.
The bigger problem though is what you are doing with your confidence interval. A CI is a statement about replication. A 95% confidence level means that in 100 replications of the experiment using similar data, 5 of the generated CIs -- which will all have different endpoints -- will _not_ contain the population parameter, although IIRC this math is more complicated in practice, meaning that the error rate is actually higher. As such, if you generate a CI and multiply the endpoints by some constant, that's a complete violation of what is being expressed: there is vastly more data with 100m driving miles than 3m miles, which will cause the CI to shrink and the estimate of the parameter to become more accurate. There is absolutely no basis for multiplying the endpoints of a CI!
Ultimately, given that the size of the sample has an effect on CI width, you need to conduct an appropriate statistical test to compare the estimated parameters between the 1 in 3m deaths for Uber and whatever data generated the 1.18 in 100m deaths for sober drivers. There's a lot more that needs to be taken into account here than what a simple Poisson test can do.
For an analysis of how AVs with various safety levels perform in terms of lives saved over time, I recommend https://www.rand.org/blog/articles/2017/11/why-waiting-for-p...
Edit: Note the default values of the T and r parameters when you run poisson.test(1, conf.level = 0.95), and also that the p-value of the one-sample exact test you performed is 1. Also, since this is an exact test, the rate of rejecting true null hypotheses at 0.95 is 0.05, but given my reservations about the use of a Poisson distribution here, I don't think that using an exact Poisson test is appropriate.
Null hypothesis: The true rate of events is 1 (r) with a time base of 1 (T).
Alternative hypothesis: The true rate is not equal to 1.
The reason that you end up with a p-value of 1 is because you've said that you've observed 1 event in a time base of 1 with a hypothesized rate of 1. So given this data, of course the probability of observing a rate equal to or more extreme than 1 is 1! As such, you're not actually testing anything about the data that you claim you are testing.
I'm not trying to be harsh here, but please be careful when using statistics!
> poisson.test(c(1, 11800), c(3, 1000000), alternative = c("two.sided"),conf.level = .93)
Comparison of Poisson rates
data: c(1, 11800) time base: c(3, 1e+06)
count1 = 1, expected count1 = 0.035403, p-value = 0.03478
alternative hypothesis: true rate ratio is not equal to 1
93 percent confidence interval:
1.006334 146.142032
sample estimates:
rate ratio
28.24859
The lower bound of the CI approaches a rate ratio = 1 for a 93% confidence interval.Interestingly, if you multiply the CI I claimed before by the rate ratio instead of the expected rate, you get almost exactly the same CI as here.
> ci <- c(0.03562718, 5.17251332)
> 28.24859 * ci
[1] 1.006418 146.116208
* Note 11800 is about two years of pedestrian deaths and time units are in millions of miles.
https://crashstats.nhtsa.dot.gov/Api/Public/ViewPublication/...> Uber is less safe than sober drivers
But the 1.18 deaths in 100m miles is for all drivers, not just the subset of sober drivers. Not quite sure why you are claiming it is only sober drivers.
It’s more like 10 miles of observation 300,000 times. Or rather an attentive human can be 50x better than average.
That's like saying: "I rolled this red d20 twenty times before I rolled a 1, whereas I rolled a 1 the first time on this blue d20, so the red d20 is obviously better and I'm rolling all my saves on it".
Or, I don't know- "I rolled three 1s on this d20 in twenty rolls so it's obviously not a fair d20".
If you have a strong prior the dice are equivalent, then yes, the rolls shouldn't change your mind.
If you have a prior that the dice are weighted in an unknown way, then yes, the rolls really should change your mind.