Project a matrix onto the hyperplane orthogonal to a unit-norm axis

Computes data - (data %*% axis) %*% t(axis). The caller is responsible for ensuring axis is unit-norm.

lq_deflate(data, axis)

Arguments

data

Numeric matrix (n_units x n_dims)

axis

Numeric vector of length n_dims, unit-norm

Value

Numeric matrix of the same shape as data