Как написать решение (X'X) в R более лаконично?

Question

Я думаю, что помню, что я видел сокращение от solve(t(X) %*% X) в R, но я не помню, что это было. Есть что-то в этом роде? Только способ сделать это за меньшее количество нажатий клавиш?

Answer 1

Возможно, вы думаете, crossprod()? Это не так много нажатий клавиш, но немного более изящно и, according to its help file, это может быть немного быстрее, чем наивная конструкция.

x <- matrix(rnorm(9), ncol=3) 
solve(crossprod(x)) 
#    [,1]  [,2]  [,3] 
# [1,] 1.34638151 -0.02957435 0.8010735 
# [2,] -0.02957435 0.32780020 -0.1786295 
# [3,] 0.80107345 -0.17862950 1.4533671 
solve(t(x) %*% x) 
#    [,1]  [,2]  [,3] 
# [1,] 1.34638151 -0.02957435 0.8010735 
# [2,] -0.02957435 0.32780020 -0.1786295 
# [3,] 0.80107345 -0.17862950 1.4533671