Pythonでの法則分布の解答

Question

私は、power lawの分布によく似たデータを持っています。パイソンを使用して、私は形で二つの式を解くことによって、データを近似する：

yは、y軸データです。 Pythonではdata[i]となります。 xはi + 1となります。ただ

解決まで

問題が来る：私たちが最初のデータのインデックスで、データでどこかに「ランダム」第二1で2つの未知の変数を持つ2つの方程式を得ることになります数学的簡略化のために

です。私はnumpy.linalg.solveのようなライブラリを使ってこのような方程式を解く方法を知らない。 Pythonを使用してaの値を見つけるにはどうすればよいですか？

Answer 1

さて、私はそれを得ました。

import math 

def get_power_law_variables(data): 
    c = data[0] 
    middle_index = len(data)/2 
    division = float(data[middle_index])/c 
    logarithm_base = middle_index + 1 
    a = math.log(division, logarithm_base) 
    return c, a 

# Example usage 
data = range(50, 150) 
c, a = get_power_law_variables(data) 
print c, a