In this notebook I describe evaluating the rank of a model matrix X by forming a pivoted Cholesky factorization of X'X. I use this form because the X'X matrix will be evaluated whether or not the rank is needed. I would appreciate if those who know the numerical linear algebra involved could tell …

Your link isn’t rendered correctly (use nbviewer), but realistically if you want a response you should post the code directly

You have: As explained in the documentation for dpstrf, this value [tol] should be negated to be used as a relative tolerance. Otherwise it is used as an absolute tolerance which usually is not what you want. This is not quite right. If tol is less than zero, dpstrf uses the recommended default v…

Usually one cares more about identifying linear dependent variables in the model matrix rather than just the rank. This way, if there are linearly dependent variables, these can be dropped such that one obtains a model matrix that is full rank. If the approach to solving the problem is through inver…

[image] Nosferican: This way, if there are linearly dependent variables, these can be dropped such that one obtains a model matrix that is full rank. This approach can be somewhat brittle (the estimates change discontinuously as a function of the data). In practice, using a vague but proper p…

I use it mostly for instrumental variable estimators such as 2SLS or after transformations such as fixed effects. In those cases, it makes sense for the user to specify a model which for some estimations the model matrix has to be trimmed to a subset which is full rank.

@Ralph_Smith Thank you for the correction on my reading of the documentation. The reason that I chose the square root of eps() is because the positive semidefinite matrix is constructed as X'X and it is rank deficiency of X that I want to detect. This also related to the next comment by @Nosferica…

@Tamas_Papp I agree that the use of hierarchical/multilevel models or, more generally, mixed-effects models is more robust. This calculation is being performed within the MixedModels package whose purpose is to fit exactly those models. I have been contacted by users who have some factors they ar…

Tolerance specification for pivoted Cholesky

Specific Domains Numerics

Topic		Replies	Views
Tolerance for numerical-rank computation Numerics linearalgebra	9	326	December 3, 2024
Cholesky Decomposition of a Sparse Symmetric Positive Semidefinite (SPSD) Singular Matrix Numerics linearalgebra , numerics , sparse , factorization	19	592	September 23, 2024
Cholesky decomposition of low-rank positive-semidefinite matrix General Usage linearalgebra , numerics	3	2822	October 26, 2021
Pivoted Cholesky: Why is rank a private field/property? Numerics linearalgebra	6	234	November 23, 2024
Singular exception when it shouldn't happen, tol for inv? Optimization (Mathematical) linearalgebra	1	1698	October 8, 2018

Tolerance specification for pivoted Cholesky

Related topics