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])