Matrix ATA has a special name in the context of least squares problems. This matrix is called the Fisher information matrix, named after the famous statistician.
Imagine measuring how sharp the peak of a mountain is. The sharper the peak, the easier it is to determine the exact location of the peak. Similarly, the Fisher information matrix provides a measure of how well we can determine the optimal parameters.
Matrix C=(ATA)−1 is the covariance matrix of the parameter estimator x^=(ATA)−1ATb. This matrix applies when we assume that components bi for i=1,…,n are independent values that are standard normally distributed.
With this assumption, the estimator x^ follows a multivariate normal distribution
x^∼N(xtrue,C)
where xtrue∈Rn is the unknown true parameter as the expected value and C∈Rn×n as the covariance matrix.
Diagonal elements cii describe the variance of parameters, like measuring how far parameter estimates can deviate from their true values. From these values, confidence intervals for the parameters can be calculated. Off-diagonal elements cij with i=j are covariances that show how the uncertainties of two parameters are related. From these covariances, correlations cij/cii⋅cjj between parameters can be obtained.
What matters in parameter estimation is not only the estimator x^ itself, but also its statistical significance as described by the covariance matrix C. Like a doctor who not only provides test results, but also explains the level of confidence in those results. In statistics courses, these concepts are discussed in more detail.
To meet requirements regarding measurement errors and provide appropriate weights to measurement data, weighted least squares problems are commonly used
xmini=1∑mσi2(h(ti)⋅x−yi)2
=∥Ax−b∥22
This problem can be transformed by defining
A=Σ−1h(t1)⋮h(tm)
b=Σ−1y1⋮ym
with
Σ−1=1/σ10⋮0⋱⋯1/σm
Here σi2 is the variance of measurement errors yi that are independent and normally distributed. Additionally, it is assumed that measurement errors have expected value 0, so there are no systematic errors. Thus bi is standard normally distributed.
In weighted least squares functions, measurement values with large measurement errors are given weaker weights compared to measurement values with small measurement errors. Think of listening to opinions from several sources. We give greater weight to more reliable sources and smaller weight to less accurate sources.