Periodicity of random kitchen sink feature mappings

Question

In various papers, e.g. Random Features for Large-Scale Kernel Machines, Rahimi and Recht introduce the now popular methodology wherein a "low rank" approximation to a stationary, PSD, kernel $K(x,y) = <\phi(X),\phi(y)>$ is constructed by randomly sampling a fourier basis from the spectral density $p$ of the kernel K such that: $K(x,y) \simeq <Z(x),Z(y)>$, where $Z(x) = \frac{2}{\sqrt{d}}[\cos(w_{1}^Tx + b_{1}),\cos(w_{2}^Tx + b_{2}),...,\cos(w_{d}^Tx + b_{d})], w_i \sim p, b \sim U(0,2\pi)$.

I am confused about the potential fact that the approximation to the potentially non periodic feature mapping $\phi$ is now approximated by a vector of periodic features. Thinking of this as a fourier series, I would think someone might specify an interval over which this approximation is valid, though no such discussion appears in the literature so I think I'm missing something. Should I be thinking of this method as approximating the kernel (i.e. the inner product between approximate feature mappings) rather than as well approximating the actual feature mapping itself?

Answer 1

Indeed, this approximation is only valid over some region near the origin. I don't have time to do so right now, but it might be informative to simply plot the value of $K(0, x)$ in $\mathbb R$ for a moderate number of random features and a wide region.

I'd dispute, though, that "no such discussion appears in the literature." Even in the original paper, the approximation theorem is for a compact region, and contains an error term based on the diameter of that region.

This plot from my 2015 paper demonstrates the increase in empirical error; here $\breve z$ gives the embedding you used above and $\tilde z$ the (better) embedding which uses sin and cos of half as many frequencies, rather than cos with an offset.

Sriperumbudur and Szábo (2015) gave a bound with the (tighter) optimal-rate dependence on the radius of the approximation set.