Matplotlib-এ হিস্টোগ্রামের বিনগুলো এমনভাবে বানান যাতে প্রতিটি বিনের এলিয়া (ফ্রিকোয়েন্সি × বিন-ওয়IDTH) সমান হয়। কাস্টম equal_area ফাংশন ব্যবহার করে বিন এজ ক্যালকুলেট করুন।
অ্যালগরিদম
- ডেটা সোর্ট করুন এবং ডিফারেন্স (dx) বের করুন
- dx^0.5-এর ক্যুমুলেটিভ সাম (tmp) বের করুন
- tmp-কে ইন্টারপোলেট করে সমান এলিয়ার বিন এজ পান
উদাহরণ কোড
import numpy as np
from matplotlib import pyplot as plt
plt.rcParams["figure.figsize"] = [7.50, 3.50]
plt.rcParams["figure.autolayout"] = True
def equal_area(x, nbin):
pow = 0.5
dx = np.diff(np.sort(x))
tmp = np.cumsum(dx ** pow)
tmp = np.pad(tmp, (1, 0), 'constant')
return np.interp(np.linspace(0, tmp.max(), nbin + 1), tmp, np.sort(x))
x = np.random.randn(1000)
n, bins, patches = plt.hist(x, equal_area(x, 20), edgecolor='black')
plt.show()আউটপুট
২০টি বিনের হিস্টোগ্রাম দেখাবে যেখানে প্রতিটি বিনের এলিয়া (আয়তন) প্রায় সমান। সাধারণ হিস্টোগ্রামে বিন-ওয়IDTH সমান হয়, এখানে এলিয়া সমান।
