R/spginv.R

In rezaalavian/LRQMM: Fitting Linear Quantile Regression Mixed Models with Relationship Matrix

spginv<-function (x)
{
  Xsvd <- sparsesvd::sparsesvd(x)
  Positive <- Xsvd$d > max(sqrt(.Machine$double.eps) * Xsvd$d[1L], 0)
  if (all(Positive))
    Xsvd$v %*% (1/Xsvd$d * t(Xsvd$u))
  else if (!any(Positive))
    array(0, dim(x)[2L:1L])
  else Xsvd$v[, Positive, drop = FALSE] %*% ((1/Xsvd$d[Positive]) * t(Xsvd$u[, Positive, drop = FALSE]))
}