I would like to fit data with integral function(truncated gamma distribution).
I tried following code, but errors occur. I am appreciate if you would kind help me. Thank you very much in advance.
%matplotlib inline
import numpy as np
from scipy import integrate
import scipy.optimize
import matplotlib.pyplot as plt
xlist=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14]
ylist=[1.0, 0.7028985507246377, 0.4782608695652174, 0.36231884057971014,
0.2536231884057971, 0.1811594202898551, 0.12318840579710147,
0.08695652173913046, 0.057971014492753645, 0.04347826086956524,
0.02173913043478263, 0.007246376811594223]
xdata=np.array(xlist)
ydata=np.array(ylist)
parameter_initial=np.array([0.0,0.0,0.0])#a,b,c
def func(x,a,b,c):
return integrate.quad(lambda t:t^(a-1)*np