where q i and q j are the estimated emission rates of naive events E i and E j , q is a vector of estimated emission rates for all naive events, and P 5 () and P 95 () are functions returning the 5 th and 95 th percentiles. If E j has a different source location estimate than E i or is separated by a period of information, then we set P i,j = 0.