استاد متوجه نشدم. این کد رو نوشتم. ممنون میشم روی کد راهنمایی بفرمایید:

import numpy as np

import pandas as pd

from scipy.stats import norm

from scipy.special import kl_div

def kl(p,q):

    result= np.sum(np.where(np.logical_and(p!=0,q!=0),p*np.log(p/q),0))

    return result

data= pd.read_csv('E:\\\\Python\\\\DataSets\\\\weight-height.csv')

males= data[data['Gender']=='Male']['Height']

m_mean= males.describe()['mean']

m_std= males.describe()['std']

count, division= np.histogram(males, bins=100)

count= count/len(males)

idata= norm.rvs(size=len(males),loc=m_mean,scale=m_std)

icount, idivision= np.histogram(idata, bins=100)

icount= icount/len(idata)

print(kl(count, icount))

print(kl_div(count, icount))