# Feedforward backpropagation

Feedforward backpropagation is an error-driven learning technique popularized in 1986 by David Rumelhart (1942-2011), an American psychologist, Geoffrey Hinton (1947-), a British informatician, and Ronald Williams, an American professor of computer science.[1] It is a supervised learning technique, meaning that the desired outputs are known beforehand, and the task of the network is to learn to generate the desired outputs from the inputs.

## Model

Model of a neuron. j is the index of the neuron when there is more than one neuron. The activation function for feedforward backpropagation is sigmoidal.

Given a set of k-dimensional inputs represented as a column vector:

$\vec{x} = [x_1, x_2, \cdots, x_k]^T$

and a nonlinear neuron with (initially random, uniformly distributed between -1 and 1) synaptic weights from the inputs:

$\vec{w} = [w_1, w_2, \cdots, w_k]^T$

then the output of the neuron is defined as follows:

$y = \varphi \left ( \vec{w}^T \vec{x} \right ) = \varphi \left ( \sum_{i=1}^k w_i x_i \right )$

where $\varphi \left ( \cdot \right )$ is a sigmoidal function. We will assume that the sigmoidal function is the simple logistic function:

$\varphi \left ( \nu \right ) = \frac{1}{1+e^{-\nu}}$

This function has the useful property that

$\frac{\mathrm{d} \varphi }{\mathrm{d} \nu} = \varphi \left ( \nu \right ) \left ( 1 - \varphi \left ( \nu \right ) \right )$

Feedforward backpropagation is typically applied to multiple layers of neurons, where the inputs are called the input layer, the layer of neurons taking the inputs is called the hidden layer, and the next layer of neurons taking their inputs from the outputs of the hidden layer is called the output layer. There is no direct connectivity between the output layer and the input layer.

If there are $N_I$ inputs, $N_H$ hidden neurons, and $N_O$ output neurons, and the weights from inputs to hidden neurons are $w_{Hij}$ ($i$ being the input index and $j$ being the hidden neuron index), and the weights from hidden neurons to output neurons are $w_{Oij}$ ($i$ being the hidden neuron index and $j$ being the output neuron index), then the equations for the network are as follows:

\begin{align} n_{Hj} &= \sum_{i=1}^{N_I} w_{Hij} x_i, j \in \left \{ 1, 2, \cdots, N_H \right \} \\ y_{Hj} &= \varphi \left ( n_{Hj} \right ) \\ n_{Oj} &= \sum_{i=1}^{N_H} w_{Oij} y_{Hi}, j \in \left \{ 1, 2, \cdots, N_O \right \} \\ y_{Oj} &= \varphi \left ( n_{Oj} \right )\\ \end{align}

If the desired outputs for a given input vector are $t_j, j \in \left \{ 1, 2, \cdots, N_O \right \}$, then the update rules for the weights are as follows:

\begin{align} \delta_{Oj} &= \left ( t_j-y_{Oj} \right )\\ \Delta w_{Oij} &= \eta \delta_{Oj} y_{Hi} \\ \delta_{Hj} &= \left ( \sum_{k=1}^{N_O} \delta_{Ok} w_{Ojk} \right ) y_{Hj} \left ( 1-y_{Hj} \right ) \\ \Delta w_{Hij} &= \eta \delta_{Hj} x_i \end{align}

where $\delta_{Oj}$ is an error term for output neuron $j$ and $\delta_{Hj}$ is a backpropagated error term for hidden neuron $j$.

## Derivation

We first define an error term which is the cross-entropy of the output and target. We use cross-entropy because, in a sense, each output neuron represents a hypothesis about what the input represents, and the activation of the neuron represents a probability that the hypothesis is correct.

$E = -\sum_{j=1}^{N_O} t_j \ln y_{Oj} + \left ( 1-t_j \right ) \ln \left (1 - y_{Oj} \right )$

Next, we determine how the error changes based on changes to an individual weight from hidden neuron to output neuron:

\begin{align} \frac{\partial E }{\partial w_{Oij}} &= \frac{\partial E }{\partial y_{Oj}} \frac{\mathrm{d} y_{Oj} }{\mathrm{d} n_{Oj}} \frac{\partial n_{Oj}}{\partial w_{Oij}} \\ &= - \left [ \frac{t_j}{y_{Oj}} - \frac{1-t_j}{1-y_{Oj}} \right ] \frac{\mathrm{d} \varphi }{\mathrm{d} n_{Oj}} y_{Hi} \\ &= - \left [ \frac{t_j}{y_{Oj}} - \frac{1-t_j}{1-y_{Oj}} \right ] \varphi \left ( n_{Oj} \right ) \left ( 1 - \varphi \left ( n_{Oj} \right ) \right ) y_{Hi} \\ &= - \left [ \frac{t_j}{y_{Oj}} - \frac{1-t_j}{1-y_{Oj}} \right ] y_{Oj} \left ( 1-y_{Oj} \right ) y_{Hi} \\ &= \left ( y_{Oj} - t_j \right ) y_{Hi} \end{align}

We then want to change $w_{Oij}$ slightly in the direction which reduces $E$, that is, $\Delta w_{Oij} \propto - \partial E / \partial w_{Oij}$. This is called gradient descent.

$\Delta w_{Oij} = - \eta \left ( y_{Oj} - t_j \right ) y_{Hi} = \eta \left ( t_j - y_{Oj} \right ) y_{Hi} = \eta \delta_{Oj} y_{Hi}$

We do the same thing to find the update rule for the weights between input and hidden neurons:

\begin{align} \frac{\partial E }{\partial w_{Hij}} &= \frac{\partial E }{\partial y_{Hj}} \frac{\mathrm{d} y_{Hj} }{\mathrm{d} n_{Hj}} \frac{\partial n_{Hj}}{\partial w_{Hij}} \\ &= \left ( \sum_{k=1}^{N_O} \frac{\partial E }{\partial y_{Oj}} \frac{\mathrm{d} y_{Oj} }{\mathrm{d} n_{Oj}} \frac{\partial n_{Oj}}{\partial y_{Hj}} \right ) \frac{\mathrm{d} y_{Hj} }{\mathrm{d} n_{Hj}} \frac{\partial n_{Hj}}{\partial w_{Hij}} \\ &= \left ( \sum_{k=1}^{N_O} \left ( y_{Oj} - t_k \right ) w_{Ojk} \right ) \frac{\mathrm{d} \varphi }{\mathrm{d} n_{Hj}} x_i \\ &= \left ( \sum_{k=1}^{N_O} \left ( y_{Oj} - t_k \right ) w_{Ojk} \right ) y_{Hj} \left ( 1 - y_{Hj} \right ) x_i \\ &= \left ( \sum_{k=1}^{N_O} - \delta_{Ok} w_{Ojk} \right ) y_{Hj} \left ( 1 - y_{Hj} \right ) x_i \end{align}

We then want to change $w_{Hij}$ slightly in the direction which reduces $E$, that is, $\Delta w_{Hij} \propto - \partial E / \partial w_{Hij}$:

\begin{align} \Delta w_{Oij} &= - \eta \left ( \sum_{k=1}^{N_O} - \delta_{Ok} w_{Ojk} \right ) y_{Hj} \left ( 1 - y_{Hj} \right ) x_i \\ &= \eta \left ( \sum_{k=1}^{N_O} \delta_{Ok} w_{Ojk} \right ) y_{Hj} \left ( 1 - y_{Hj} \right ) x_i \\ &= \eta \delta_{Hj} y_{Hj} x_i \end{align}

## References

1. Script error: No such module "Citation/CS1".