Tutorials › Multivariate Statistics (Graduate) › Parameter Counts Under Covariance Constraints

Model-Based Clustering · Tutorial 659 of 975

Parameter Counts Under Covariance Constraints

Learn to count the free covariance parameters in Gaussian clustering models and to distinguish covariance counts from full mixture-model counts.

Advanced 10 min read

What You'll Learn

  • Count covariance parameters for component-specific and shared spherical models.
  • Derive the parameter counts for diagonal covariance structures.
  • Count parameters when cluster shapes are shared but orientations vary.
  • Explain the dimension of an unrestricted covariance matrix and an orthogonal orientation.
  • Separate covariance counts from means and mixture proportions in a full model.
  • Recognize when repeated eigenvalues make orientation parameters non-identifiable.

Why Covariance Parameter Counts Depend on Constraints

A Gaussian clustering model assigns each cluster its own mean and covariance structure, or imposes selected constraints to reduce complexity. Counting the covariance parameters makes those constraints precise: it tells us how many independent quantities must be specified, not how many entries happen to appear in a matrix. This distinction matters when describing a model and when assessing its overall dimension.

Let \(p\) be the dimension of each observation and \(G\) the number of clusters. Write \(\Sigma_g\) for the positive-definite covariance matrix of cluster \(g\). Unless a constraint is stated, a covariance may differ across clusters. We will count covariance parameters first, then show how to add the means and mixing proportions in a finite Gaussian mixture.

Definition: The covariance-parameter count is the number of independent continuous parameters needed to specify all the covariance matrices under the stated model constraints. It is a dimension count on the model's regular, identifiable part; discrete choices such as eigenvector signs do not add continuous parameters.

Three Basic Dimensions

A symmetric \(p\times p\) matrix has \(p\) diagonal entries and \(p(p-1)/2\) distinct off-diagonal entries. Thus an unrestricted symmetric matrix has \(p(p+1)/2\) free entries. Positive definiteness restricts which values those entries may take, but does not change this dimension: positive-definite matrices form an open subset of the vector space of symmetric matrices.

A diagonal covariance matrix has only its \(p\) diagonal entries free. A spherical covariance matrix is more constrained still: it is a positive scalar times the identity, so it has just one scale parameter.

Result: For one cluster, an unrestricted covariance has \(p(p+1)/2\) parameters, a diagonal covariance has \(p\), and a spherical covariance has \(1\). If these structures vary independently across \(G\) clusters, multiply the corresponding one-cluster count by \(G\).

The word shared changes the count. If all clusters share one diagonal covariance, there are \(p\) covariance parameters total, not \(Gp\). If they share one spherical covariance, there is one parameter total, not \(G\).

Worked Example: Spherical and Diagonal Counts

Suppose there are \(G=3\) clusters in \(p=2\) dimensions. Compare component-specific spherical covariances, one shared spherical covariance, component-specific diagonal covariances, and one shared diagonal covariance.

For component-specific spheres, write \(\Sigma_g=\lambda_g I_2\), with \(\lambda_g>0\). Each cluster contributes one scale, so the count is \(3\).

If the sphere is shared, \(\Sigma_g=\lambda I_2\) for every \(g\). There is one common scale, so the count is \(1\).

A component-specific diagonal covariance has the form \(\Sigma_g=\operatorname{diag}(a_{g1},a_{g2})\), with both entries positive. That gives \(2\) parameters per cluster and \(3\times2=6\) total.

For a shared diagonal covariance, \(\Sigma_g=\operatorname{diag}(a_1,a_2)\) for all \(g\), giving \(2\) total. Sharing reduces the count because the same entries are used for every cluster; it does not create a separate copy for each cluster.

Counting Orientation and Shape

For covariance models described through volume, shape, and orientation, a useful parameterization is \(\Sigma_g=\lambda_g D_g A_g D_g^{\mathsf T}\). Here \(\lambda_g>0\) is a volume or scale factor, \(D_g\) is an orthogonal orientation matrix, and \(A_g\) is diagonal with positive entries and determinant \(1\). The determinant constraint separates scale from shape: \(\det(\Sigma_g)=\lambda_g^p\). The shape matrix therefore has \(p-1\) free entries.

For a generic covariance with distinct eigenvalues, the orientation has \(p(p-1)/2\) continuous parameters. The reason is that the orthogonal matrices form a space of that dimension. Although \(D_g\) has \(p^2\) entries, its columns must be orthonormal, imposing \(p(p+1)/2\) independent local constraints.

Derivation: At an orthogonal matrix \(D\), a small change preserving orthogonality has the form \(D K\), where \(K\) is skew-symmetric: \(K^{\mathsf T}=-K\). Such a matrix has zero diagonal and \(p(p-1)/2\) freely chosen entries above the diagonal; the entries below are then fixed. Thus the local dimension of orientation is \(p(p-1)/2\). The \(p\) positive diagonal entries of \(A\) obey one independent constraint, \(\prod_{j=1}^{p} A_{jj}=1\), leaving \(p-1\) shape parameters.

Signs and permutations of eigenvectors can give different matrix representations of the same covariance. These are discrete ambiguities, not additional continuous degrees of freedom. When eigenvalues are repeated, rotations within the repeated-eigenvalue subspace also leave the covariance unchanged; the orientation is then not fully identifiable. The generic orientation count applies on the distinct-eigenvalue part of the model.

Shared Shape and Unrestricted Component Covariances

A shared-shape model needs its orientation constraint stated explicitly. In the version considered here, every cluster has the same normalized shape \(A\), but each cluster may have its own scale and orientation: \(\Sigma_g=\lambda_g D_g A D_g^{\mathsf T}\). The scales contribute \(G\), the common shape contributes \(p-1\), and the \(G\) orientations contribute \(G p(p-1)/2\). Therefore the count is \(G+(p-1)+G p(p-1)/2\).

If orientation is shared as well as shape, all clusters use the same \(D\). Then the count is \(G+(p-1)+p(p-1)/2\). This is a different, more constrained model. If both scale and orientation are shared, only one common covariance shape remains, with \(p(p+1)/2\) parameters in total. Always translate a verbal constraint into an explicit equation before counting.

With unrestricted component covariances, each \(\Sigma_g\) is a free positive-definite symmetric matrix. The total covariance count is \(G p(p+1)/2\).

Result: For generic positive-definite covariance models, the principal counts are:
  • Component-specific spherical: \(G\); shared spherical: \(1\).
  • Component-specific diagonal: \(Gp\); shared diagonal: \(p\).
  • Shared shape with component-specific scales and orientations: \(G+(p-1)+G p(p-1)/2\).
  • Unrestricted component-specific: \(G p(p+1)/2\).

Worked Example: Shared Shape in Three Dimensions

Take \(p=3\) and \(G=2\), with a common shape matrix but cluster-specific scales and orientations. The scale count is \(G=2\). A determinant-one diagonal shape has \(p-1=2\) free parameters. Each orientation has \(p(p-1)/2=3\) parameters, so two orientations contribute \(2\times3=6\). The total is \(2+2+6=10\).

If orientation is common too, the count becomes \(2+2+3=7\): two scales, two shape parameters, and one three-parameter orientation. For comparison, two unrestricted covariance matrices require \(2\times3(4)/2=12\) parameters. These counts differ because the constraints differ, not because the number of clusters or dimensions has changed.

From Covariance Counts to Full Mixture Counts

In a \(G\)-component Gaussian mixture, each component also has a \(p\)-dimensional mean, giving \(Gp\) mean parameters. The mixing proportions are positive and sum to one, so they contribute \(G-1\) free parameters. Consequently, if the covariance count under a chosen structure is \(K\), the usual full parameter count is \(Gp+(G-1)+K\).

This addition assumes the ordinary interior parameterization: all mixture proportions are positive, and the component means and covariance parameters are counted as distinct continuous parameters. It is a bookkeeping identity, not a claim that every parameter is identifiable in every fitted mixture. For instance, coincident components can create non-identifiability even though the nominal dimension is unchanged.

Worked Example: Full Count for a Diagonal Mixture

Consider \(G=4\) Gaussian components in \(p=3\) dimensions, each with its own diagonal covariance. The covariance count is \(Gp=4\times3=12\). The means contribute \(Gp=4\times3=12\), and the mixing proportions contribute \(G-1=3\). Thus the nominal full count is \(12+12+3=27\).

If instead the diagonal covariance is shared across all four components, the covariance count is \(p=3\), while the means and proportions still contribute \(12\) and \(3\). The full count is \(3+12+3=18\). Only the covariance part changed; the shared diagonal restriction does not imply shared means or equal mixing proportions.

Common Mistakes and a Reliable Counting Procedure

A frequent mistake is to count every matrix entry as free. An orthogonal matrix is not an arbitrary \(p\times p\) array, and a symmetric covariance does not have \(p^2\) distinct entries. Another is to count a shared parameter once per cluster: if the model specifies one common shape or covariance, that quantity is counted once.

A subtler error is to count \(p\) shape parameters after imposing \(\det(A)=1\). The determinant constraint removes one degree of freedom, leaving \(p-1\). Conversely, if the scale has not been separated from shape by a normalization, counting both a free scale and an unconstrained shape can count the same degree of freedom twice.

1
Write the covariance equation.
Make clear which quantities vary by cluster and which are common.
2
Count each free block.
Use \(1\) for a scalar scale, \(p\) for a diagonal matrix, \(p(p+1)/2\) for a symmetric covariance, and \(p(p-1)/2\) for a generic orientation.
3
Apply constraints before multiplying.
For example, a determinant-one diagonal shape has \(p-1\) parameters; multiply by \(G\) only if that shape is component-specific.
4
Check the scope of the count.
State whether the result counts covariances alone or also includes means and mixture proportions.

Parameter counts make the meaning of a covariance constraint concrete. They also prevent misleading comparisons between models whose labels sound similar but whose equations impose different sharing patterns. The count is most reliable when derived from the covariance parameterization itself, with each equality, normalization, and cluster-specific quantity accounted for explicitly.

Check Your Understanding

For each question, identify the relevant constraint before calculating.

  1. How many parameters are in \(G\) component-specific spherical covariance matrices in dimension \(p\)? How many if the spherical covariance is shared?
  2. Why does a shared diagonal covariance contribute \(p\) parameters rather than \(Gp\)?
  3. For a common determinant-one shape and \(G\) component-specific orientations and scales, derive the covariance count in terms of \(p\) and \(G\).
  4. How many covariance parameters are there for \(G\) unrestricted covariance matrices in dimension \(p\), and what matrix property gives the one-cluster count?
  5. In a full Gaussian mixture, what terms must be added to a covariance count to include the means and mixing proportions?
Clinical Biostats