Unsupervised Learning: Randomized Optimization
Tags: machinelearning, unsupervisedlearning
Published: 2019-02-23 12:32:43

---
title: "Unsupervised Learning: Randomized Optimization"
published: true
description: Hill Climbing, Simulated Annealing, Genetic Algorithms, oh my!
tags: Machine Learning, Unsupervised Learning
---

*This is the 12th in a series of class notes as I go through the [Georgia Tech/Udacity Machine Learning course](https://www.udacity.com/course/machine-learning--ud262). The class textbook is [Machine Learning by Tom Mitchell](https://www.cs.ubbcluj.ro/~gabis/ml/ml-books/McGrawHill%20-%20Machine%20Learning%20-Tom%20Mitchell.pdf).*

> This marks the start of a new miniseries on Unsupervised Learning, [the 2nd of 3 sub disciplines within Machine Learning](https://dev.to/swyx/machine-learning-an-overview-216n).

## Our Goal: Optimization

Given an input space `X`, and an objective (aka fitness) function `f(x)`, Find `x` from `X` such that `f(x)` is the maximum possible it can be. 

This is similar to the [argmax](https://en.wikipedia.org/wiki/Arg_max) function we have used in prior formula, but now we are actually exploring methods to do it.

## When Traditional Methods fail

Given a small `X`, we can simply run a `for` loop [like this one](https://repl.it/@swyx/ML-optimization-practice) and keep track of the highest values.

Given an infinite `X` but with a solvable `df(x)`, we can use calculus to help us find optima.

If `df(x)` can't be solved, we can use the [Newton Raphson method](https://en.wikipedia.org/wiki/Newton%27s_method) to iteratively get closer and closer to the optimum. However, it can get stuck in local optimum.

But what happens when we don't have a derivative to work with? (for example, you don't even know what f(x) is, you just know the output without knowing its theoretical maximum) and possibly many local optima?

## Method 1: Random Restart Hill Climbing

The simple solution is, to guess a few starting points and then guess a direction and try to do some [random restart hill climbing](https://en.wikipedia.org/wiki/Hill_climbing). However, to do it right, you'd have to do a number of restarts and it turns out to be not a lot better than just looping through the entire space especially for spaces with narrow global optima basins.

Can we do better than that?

## Method 2: Simulated Annealing

Random Restart Hill Climbing combines **exploring** (random restart) with **exploiting** (climbing). Maybe we can split that up. 

[Simulated Annealing](https://en.wikipedia.org/wiki/Simulated_annealing) takes its name from metallurgy where repeatedly heating up and cooling down a sword makes it stronger than before.

![http://2.bp.blogspot.com/--kOlrodykkg/UbfVZ0_l5HI/AAAAAAAAAJ4/0rQ98g6tDDA/s1600/annealingAtoms.png](http://2.bp.blogspot.com/--kOlrodykkg/UbfVZ0_l5HI/AAAAAAAAAJ4/0rQ98g6tDDA/s1600/annealingAtoms.png)

The pseudocode for this looks like:

```
Start with a temperature T
For a finite set of iterations:

- sample new point x+
- jump to a new point with probability given by an acceptance probability function P(x, x+, T)
- decrease temperature T

where 

P(x, x+, T) = 1 if f(x+) >= f(x) else e^((f(x+) - f(x))/T)
```

So a high T makes the algorithm behave like a random sampler, and as it cools, it starts to behave more like a hill climber.

![https://slideplayer.com/slide/8038378/25/images/6/Convergence+of+simulated+annealing.jpg](https://slideplayer.com/slide/8038378/25/images/6/Convergence+of+simulated+annealing.jpg)

A remarkable analytical result about SA is that the probability of it ending at X is equal to `e ^ (f(x) / T) / Z(T)` ([a Boltzmann distribution](https://en.wikipedia.org/wiki/Boltzmann_distribution)) which makes the likelihood of its ending point directly related to its ending point's fitness.

## Method 3: Genetic Algorithms

Evolution is a pretty good optimizer, and we can simulate that by "breeding" individual points to produce hopefully better offspring and "mutating" them by doing local search. Over multiple "generations", the hope is that we evolve towards the globally optimum result. One feature that helps a lot is if the individual dimensions/attributes can be additively combined, so that breeding two parents can regularly produce better offspring.

![https://image.slidesharecdn.com/gasbytheseatofyourgenes-160629032611/95/genetic-algorithms-programming-by-the-seat-of-your-genes-6-638.jpg?cb=1467171105](https://image.slidesharecdn.com/gasbytheseatofyourgenes-160629032611/95/genetic-algorithms-programming-by-the-seat-of-your-genes-6-638.jpg?cb=1467171105)

To contrast this with the other methods, GAs are like random restarts done in parallel (since each individual in a population is like a thread, and each generation is like a restart), EXCEPT that the "crossover" of information when "breeding" happens helps to direct the general direction of evolution towards better places.

The pseudo-code looks like:

```
Start with an initial population of size K
Repeat until converge:
  - compute fitness of all population
  - select "most fit" individuals
  - pair up individuals, replacing "least fit" via crossover/mutation
```

As you might imagine, the secret sauce is in the breeding. Two methods of doing preserving individuals are:

- truncation selection - for example, only keeping the top half of the population
- roulette wheel selection - taking a probabilistic approach to keeping each individual according to their fitness

Crossover/Mutation also has nuances. Deciding on what attributes can be treated as a group involves domain knowledge and implicit assumptions that may or may not reflect reality.

## More Randomized Optimization

As you can see, there are many ways to tackle the problem of optimization without calculus, but all of them involve some sort of random sampling and search. 

Those we have explored don't have much in the way of memory or of actually learning the structure or distribution of the function space, but there are yet more algorithms that explore those abilities:

- [Tabu Search](https://en.wikipedia.org/wiki/Tabu_search) adds memory to the search (stay away from "taboo" regions)
https://www.cc.gatech.edu/~isbell/tutorials/mimic-tutorial.pdf
- [MIMIC: Finding Optima by Estimating Probability Densities](https://www.cc.gatech.edu/~isbell/papers/isbell-mimic-nips-1997.pdf) estimates a probability distribution in lieu of a genetic algorithm metaphor.

## MIMIC

MIMIC has a simple idea at its core - define a uniform probability distribution of points with a fitness over a threshold `theta`:

![https://i.ytimg.com/vi/49Y2C_mTWD0/maxresdefault.jpg](https://i.ytimg.com/vi/49Y2C_mTWD0/maxresdefault.jpg)

When theta is at a minimum, that means we uniformly sample over the entire space.

When theta is at a maximum, that means we converge on one or more global optima.


```
at every iteration:

- Generate samples from population consistent with the distribution so far
- set `theta` to the nth percentile
- retain only those samples where f(x) > `theta`
- estimate the new distribution based on these samples,
- repeat
```

Therefore, under MIMIC, our task is to continually raise `theta` until we just have small regions around our optima left. In the process of doing so, we infer a (uniform) probability distribution around our optima. Because we use the `nth` best percentile as our cutoff factor, this is similar to genetic algorithms, however there is no longer a concept of breeding or genes. Instead, we stay laser focused on performance against the fitness function as our one metric.

The comparison with genetic algorithms isn't accidental - the best version of MIMIC estimates distributions by assuming that the conditional distributions inherent in the space fit to "dependency trees", which are [Bayesian Networks](https://dev.to/swyx/supervised-learning-bayesian-inference-4l72) where every node (but one) has exactly one parent. Finding these trees is out of scope but [check](https://www.youtube.com/watch?v=R-Mf9-tKC5o) [the](https://www.youtube.com/watch?v=u2V7esTMhKc) [videos](https://www.youtube.com/watch?v=rCxGiKQoe4w) for the math, but the theoretical grounding of this approach is rooted in maximizing [Mutual Information](https://en.wikipedia.org/wiki/Mutual_information) (in information theory terms) and building a [Maximum Spanning Tree](http://shodhganga.inflibnet.ac.in/bitstream/10603/33822/4/chapter4.pdf). While not the only way to estimate the next probability distribution, Dependency Trees allow capturing of parent-child structure without having exponential cost (Max Spanning Trees are only polynomial in cost).


Benefits:

- The fact that MIMIC's "rising theta" approach can be used with any underlying probability distribution is pretty cool - here's a [simple applied example](https://www.youtube.com/watch?v=4DtVLCC73bA).
- Being able to learn structure via the distribution assumption helps with generalizability
- Empirical tests show MIMIC using 2-3 orders of magnitude less iterations to converge compared to Simulated Annealing and other alternatives, however, this is offset by each iteration taking much longer because of the need to estimate and generate from probability distributions. Thus, in terms of overall time taken, MIMIC works better when cost of evaluating `f(x)` is high (e.g. rocket simulation, antennae design, human test subjects)



## Next in our series

Further notes on this topic:

- [Wikipedia on Randomized Optimization](https://en.wikipedia.org/wiki/Random_optimization)
- [The original MIMIC paper](https://www.cc.gatech.edu/~isbell/papers/isbell-mimic-nips-1997.pdf)

Hopefully that was a good introduction to Randomized Optimization. I am planning more primers and would love your feedback and questions on:


- [Overview](https://dev.to/swyx/machine-learning-an-overview-216n)
- Supervised Learning
    - [Decision Trees](https://dev.to/swyx/machine-learning-classification-learning--decision-trees-1mbh)
    - [Regression](https://dev.to/swyx/supervised-learning-regression-4d17)
    - [Neural Networks](https://dev.to/swyx/supervised-learning-neural-networks-mpo)
    - [Instance Based Learning (K Nearest Neighbors)](https://dev.to/swyx/supervised-learning-instance-based-learning-and-k-nearest-neighbors-kge)
    - [Ensemble Learning (AdaBoost)](https://dev.to/swyx/supervised-learning-ensemble-learning-lim)
    - [Kernel Methods & SVMs](https://dev.to/swyx/supervised-learning-support-vector-machines-3mgk)
    - [Computational Learning Theory](https://dev.to/swyx/supervised-learning-computational-learning-theory-160h)
    - [VC Dimensions](https://dev.to/swyx/supervised-learning-vc-dimensions-10b)
    - [Bayesian Learning](https://dev.to/swyx/supervised-learning-bayesian-learning-403l)
    - [Bayesian Inference](https://dev.to/swyx/supervised-learning-bayesian-inference-4l72)
- Unsupervised Learning
    - [Randomized Optimization](https://dev.to/swyx/unsupervised-learning-randomized-optimization-4c1i)
    - [Information Theory](https://dev.to/swyx/unsupervised-learning-information-theory-recap-4iem)
    - Clustering - week of Feb 25
    - Feature Selection - week of Mar 4
    - Feature Transformation - week of Mar 11
- Reinforcement Learning
    - Markov Decision Processes - week of Mar 25
    - "True" RL - week of Apr 1
    - Game Theory - week of Apr 15