this is for holding javascript data
adam greenberg edited method.tex
about 10 years ago
Commit id: 77cb886467da6a72fde0b36994bcc154219f27d2
deletions | additions
diff --git a/method.tex b/method.tex
index edf314a..79f130d 100644
--- a/method.tex
+++ b/method.tex
...
\subsection{Steepest Descent Routine}
A classical steepest descent routine (SDR) minimizes the weighted residuals between a model and data with Gaussian noise by determining the direction in parameter space in which the $\chi^2$ is decreasing fastest. Specifically, suppose one has a set of $m$ observables, $\vec{z}$ with weights $W$, and a model function $\vec{m}(\vec{x})$, where \vec{x} is an $n$-dimensional parameter vector. Assuming independent data points with Gaussian-distributed errors, the probability of the model matching the data is given by \[p(\vec{m}(\vec{x}) | \vec{z}) \propto p(\vec{z} | \vec{m}(\vec{x})) \propto \exp(
\frac{-1}{2}\vec{R}^T \frac{-1}{2}\vec{R}^\intercal W \vec{R})\] where $\vec{R} = \vec{z} - \vec{m}(\vec{x})$