Given a matrix of genotypes, this function estimates the genetic relatedness matrix (GRM,
also known as the RRM, see Hayes et al. 2009, tools:::Rd_expr_doi("10.1017/S0016672308009981")) among
the subjects: \(\frac{1}{p}(XX^T)\), where X is standardized.