То, что ниже скрипт делает: среднесуточные с 2011 года от blockchain.info, сделать некоторые преобразования с пандами, а затем построить его. Это почти подходит, но немного слишком много для 1 и не достаточно для 2 и 3.
Исходя из этого, я предположить, что мы либо видим длительный период $ 200-299 или $ 2000-2999.
Идеи приветствуются, улучшение, конечно, возможно ...
Код: (питон)
#! / USR / бен / окр питон
импорт NumPy в нп
импорт панды как э.р.
pricedata = "http://blockchain.info/charts/market-price?showDataPoints=false&TimeSpan = все&show_header = верно&daysAverageString = 1&масштаб = 0&Формат = CSV"
Данные = pd.read_csv (pricedata, имена = ["день", "цена"], Parse_dates = [0], index_col = 0)
#start на 2011-01-01
Данные = data.ix [ '2011-01-01':]
печать data.head ()
price100 = данные ["цена"] * 100
pcolstr = price100.map (ул)
высокий = pcolstr.map (лямбда _: _ [0]). Отображение (целое)
печать highest.describe ()
вал, _ = np.histogram (highest.as_matrix (), мусорные ведра = 9)
цифр = np.arange (1,10)
Benford = np.log10 (1 + 1 / цифры)
Vals = val.astype (np.float) / highest.size
печати (np.c_ [Benford, Vals])
импорт matplotlib.pyplot в PLT
plt.plot (цифры, Benford)
plt.plot (цифры, Vals, цвет ="зеленый")
plt.show ()
импорт NumPy в нп
импорт панды как э.р.
pricedata = "http://blockchain.info/charts/market-price?showDataPoints=false&TimeSpan = все&show_header = верно&daysAverageString = 1&масштаб = 0&Формат = CSV"
Данные = pd.read_csv (pricedata, имена = ["день", "цена"], Parse_dates = [0], index_col = 0)
#start на 2011-01-01
Данные = data.ix [ '2011-01-01':]
печать data.head ()
price100 = данные ["цена"] * 100
pcolstr = price100.map (ул)
высокий = pcolstr.map (лямбда _: _ [0]). Отображение (целое)
печать highest.describe ()
вал, _ = np.histogram (highest.as_matrix (), мусорные ведра = 9)
цифр = np.arange (1,10)
Benford = np.log10 (1 + 1 / цифры)
Vals = val.astype (np.float) / highest.size
печати (np.c_ [Benford, Vals])
импорт matplotlib.pyplot в PLT
plt.plot (цифры, Benford)
plt.plot (цифры, Vals, цвет ="зеленый")
plt.show ()