# Two sample problem(2): kernel function, feature space and reproducing kernel map

Find Two sample problem (1) here.

We will take a look at RHKS (Reproducing Hilbert Kernel Space ) in this post. You might think of it a very statistical term but it is amazing because of various applications. You will need to refresh your mind for some linear algebra computations. We start with some basic terms and definitions.

#### Statistical Distance

It measures the distance between two statistical objects, i.e two random variables, two probability distributions or samples.
It is simple for you to calculate the distance of two vectors, for example, cosine similarity is a common way. But when it comes to the divergence of two probability distributions, things are more complicated, the Kullback-Leibler Divergence(also called Relative Entropy) is a method. We use $D(P||Q)$ to denote the amount of information lost where P is the true distribution and Q is used to approximate P, $D(P||Q) \neq D(Q||P)$, note that it is not symmetric. Another impressive term is Maximum Mean Discrepancy (MMD), which is defined according to kernel embedding of distributions.

#### Kernel Function and Feature Space

A feature mapping is a function $\phi$ that maps data from an original space X to another space, where usually we call it feature space H. Keep in mind that we have $x \rightarrow \phi(x)$. And x, x’ refer to different samples or dots.
From the above picture, we are mapping data from a 2-D to a 3-D feature space. Given the certain mapping function, if we calculate the inner product of two dots in the feature space, we surprisingly find that it is just the square of the inner product of the two dots (x and x’) from the original space. We could simply define a kernel function $\kappa(x,x') = ()^2$ and eliminate the mapping function $\phi$.

##### Proposition

Also called Cauchy-Schwarz Inequality for kernels:

$|\kappa (x,z)|^{ 2 }\leq \kappa (x,x)\kappa (z,z)$

Proof: The Gram Matrix/ Kernel Matrix for x and z is positive semi-definite matrix.

$\kappa =\begin{bmatrix} \kappa (x,x) & \kappa (z,x) \\ \kappa (z,x) & \kappa (z,z) \end{bmatrix}=\begin{bmatrix} \kappa (x,x) & \kappa (x,z) \\ \kappa (x,z) & \kappa (z,z) \end{bmatrix}\geq 0$

So determinant is non-negative:

$det(\kappa )=\begin{vmatrix} \kappa (x,x) & \kappa (z,x) \\ \kappa (z,x) & \kappa (z,z) \end{vmatrix} = \kappa(x,x)\kappa(z,z)-\kappa^2(x,z)\geq0$

##### Distance and angle in feature space

In feature space, we could define distance between two sample:

$||\phi (x)-\phi (x')||^{ 2 }=(\phi (x)-\phi (x'))^{ T }(\phi (x)-\phi (x'))\\ =<\phi (x),\phi (x)>-2<\phi (x),\phi (x')>+<\phi (x'),\phi (x')>\\ =\kappa (x,x)-2\kappa (x,x')+\kappa (x',x')$

Inner product is defined as:

$<\phi (x),\phi (x')>=||\phi(x)||\cdot||\phi(x')||\cos\theta$

So simply calculate:

$\cos \theta =\frac { <\phi (x),\phi (x')> }{ ||\phi (x)||\cdot ||\phi (x')|| }\\ =\frac { <\phi (x),\phi (x')> }{ \sqrt { <\phi (x),\phi (x)> } \sqrt { <\phi (x'),\phi (x')> } }\\ =\frac{\kappa(x,x')}{\sqrt{\kappa(x,x)}\sqrt{\kappa(x',x')}}$

So we will notice that both distance and angle are in the representation of kernel functions only, without mapping function $\phi$

#### Reproducting Kernel Map

A set of functions from feature space $R^d$ to $R$.

$R^{ R^{ d } }=\{ f:R^d\rightarrow R \}$

The Reproducing kernel map is a map: $\Phi: R^d \rightarrow R^{R^d}$, while $\Phi(z) = \phi_z = \kappa(\cdot,z)$. More concretely, $\Phi(z)(x)=\phi_z(x)=\kappa(x,z)$. So “reproducing” simply “replace” a vector to get new functions,which is a process of re-producing new relations. We can treat it as, given a function who has a lot of parameters, we replace some of them but keep the rest. Every time we speficy a group of parameters, we will have a new function.
Here is an example of how it works:

Given the kernel function $\kappa(x,z)$ to be a Gaussian function who has two parameters x and z. In practice, x and z could be vectors, but in this example, we think them to be only a scalar. If we let $z=0$, we will have a function $exp(-\frac{x^2}{2})$. If we let $z=1$, we will have a function $exp(-\frac{(x-1)^2}{2})$. If you set a group of z values, you will get a set of corresponding functions.