Past exam of the mathematics course of the University of Cambridge 2016 iii Paper 205 4 Solution Created 2026-10-03 Updated 2026-10-06
Write for a positive-definite matrix. Up to a constant independent of the parameters, the log-likelihood of the multivariate normal distribution isThe identityshows that, for each , the unique optimizing mean is . Using the matrix trace to rewrite the first quadratic sum, the profile likelihood isTherefore maximum likelihood estimation of the precision matrix, whenever the maximum exists, is exactlyThis is the profile likelihood for a Gaussian precision matrix. The sample covariance matrix here uses divisor , as appropriate for maximum likelihood estimation.
There is an important rank qualification in the printed assumptions. If , differentiation gives , so the unique minimizer is . Uniqueness follows because the negative log-determinant is a strictly convex function on the positive-definite matrices. But full column rank of the uncentered alone does not imply : the rank of a centered sample covariance matrix depends on the centered design. For example, with and raw data matrix ,Although has full column rank, setting gives objective . Thus an unpenalized maximizer need not exist under the literal raw-rank assumption. The likelihood reduction remains valid; existence requires full column rank of the centered data, equivalently . For the nonsingular Gaussian model this holds almost surely when .
Use the entrywise-penalty convention for the Graphical Lasso:The sum is the entrywise matrix L1 norm. Its Graphical-Lasso Karush-Kuhn-Tucker conditions areHere is symmetric; its nonzero-entry values are , and its zero-entry values are any numbers in . These are the subgradient optimality condition for this convex optimization problem. Since a positive-definite matrix has positive diagonal, and . Counting both off-diagonal entries in the penalty and in the trace makes the same appear in their stationarity equations. The common variant penalizing only off-diagonal entries instead sets ; the off-diagonal equations and the column argument below are unchanged.
For the given column problem, abbreviateIts smooth part expands to plus a constant, because is symmetric and . The matrix is a positive-definite matrix, so this quadratic is a strictly convex function. Adding the convex L1 norm penalty preserves strict convexity, and the objective tends to infinity as . The minimizer exists and is unique. Its Karush-Kuhn-Tucker conditions are
We identify that minimizer directly from the precision matrix. Write and defineThe denominator is positive. The off-diagonal part of column of givesOn the other hand, the off-diagonal Graphical-Lasso Karush-Kuhn-Tucker conditions readBecause , the coordinates of have the opposite signs from those of and vanish at exactly the same positions. The subgradient of the absolute value therefore satisfiesIt follows thatThese are exactly the Lasso column problem's Karush-Kuhn-Tucker conditions. Its unique minimizer must be , yieldingThis is the Graphical Lasso column regression identity. The minus sign is essential: it converts the precision-entry subgradient into the regression-coefficient subgradient. The argument uses the inverse identity directly, so it does not require proving the supplied block inverse formula.