def gaussian_2d_plain(): lim = 3.5 indexes = utils.normalize_to_indexes(low=[-lim, -lim], high=[lim, lim], n=100) p = support_gaussian_2d(indexes=indexes) plotting.plot_combined(p, indexes[0], indexes[1], k=[3, 5, 7, 10])
def gaussian_2d_central_splike(): lim = 3.5 indexes = utils.normalize_to_indexes(low=[-lim, -lim], high=[lim, lim], n=100) p1 = support_gaussian_2d(indexes=indexes, sigma=[[1.0, 0], [0, 1]]) p2 = support_gaussian_2d(indexes=indexes, sigma=[[0.005, 0], [0, 0.005]]) plotting.plot_combined(p1 + p2, indexes[0], indexes[1], k=[3, 5, 10])
def gaussian_2d_shifted_spike(): lim = 3 indexes = utils.normalize_to_indexes(low=[-lim, -lim], high=[lim, lim], n=100) p1 = support_gaussian_2d(indexes=indexes, sigma=[[1.0, 0], [0, 1]]) p2 = support_gaussian_2d(indexes=indexes, mu=[0.25, 0.25], sigma=[[0.02, 0], [0, 0.02]]) plotting.plot_combined(p1 + p2, indexes[0], indexes[1], k=[3, 5, 10])
def titanic_kde(kernel_bandwidth=None): # get data titanic = pd.read_csv('./data/titanic_mixed.csv', index_col=None, usecols=['Age', 'Fare']) data = titanic.values.transpose() # derive pdf mykernel = pdf_kernel(data, kernel_bandwidth=kernel_bandwidth) # get support indexes = utils.normalize_to_indexes(data=data) p = support_kde(mykernel, indexes) # plot plotting.plot_combined(p, indexes[0], indexes[1], k=[3, 5, 7, 10])
def gausssian_2d_three_gaussians(): lim = 2.5 indexes = utils.normalize_to_indexes(low=[-lim, -lim], high=[lim, lim], n=100) p1 = support_gaussian_2d(indexes=indexes, sigma=[[1, 0], [0, 1]]) p2 = support_gaussian_2d(indexes=indexes, mu=[0.25, 0.25], sigma=[[0.02, 0], [0, 0.02]]) p3 = support_gaussian_2d(indexes=indexes, mu=[-0.35, -0.35], sigma=[[0.07, 0], [0, 0.1]]) plotting.plot_combined(2 * p1 + p2 + p3, indexes[0], indexes[1], k=[3, 5, 10])
def titanic(): df = pd.read_csv('./data/titanic_age-vs-fare.csv', index_col=False) titanic_p = df['p'].values plotting.plot_combined(titanic_p, k=[3, 5, 7, 10])
def allbus(): df = pd.read_csv('./data/allbus_age-vs-income.csv', index_col=False) allbus_p = df['p'].values plotting.plot_combined(allbus_p, k=[3, 5, 7, 10])