我正在尝试使用PyMC3,试图在PyMC3文档中适应的修改版本。假设您有两个煤矿(mine1和mine2),每个煤矿在相同的年份范围内都有类似的灾难计数。
然而,mine1在实施减少灾难数量的安全程序方面晚了5年:
import numpy as np
import matplotlib.pyplot as plt
mine1=np.array([0,4,5,4,0,1,4,3,4,0,6,3,3,4,0,2,6,3,3,5,4,5,3,1,4,4,1,5,5,3,4,2,5,2,2,3,
4,2,1,3,0,2,1,1,1,1,3,0,0,1,0,1,1,0,0,3,1,0,