Massey 行列の調整
WNBA の Massey 行列モデルでは、レーティング問題に解が存在し一意になるよう、いくつか調整が必要です。現在の行列は（計算上）可逆ではないためです。
これを変更する一つの方法は、行列 \(M\) の一番下に 1 の行を追加し、\(M\) の一番右に -1 の列を追加し、得点差ベクトル \(\vec{f}\) の一番下に 0 を加えることです。この演習では、これらすべてを実際に行列 \(M\) に追加します。
この演習はコースの一部です
Rで学ぶデータサイエンスのための線形代数
演習の手順
rbind()を使って、行列Mの一番下に1の行を追加し、結果をM_2に代入します。cbind()を使って、M_2の右側に-1の列を追加し、結果をM_3に代入します。M_3の 13x13 要素を1に設定します。- 得られた行列を出力します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Add a row of 1's
M_2 <- rbind(M, ___)
# Add a column of -1's
M_3 <- ___(M_2, rep(-1, 13))
# Change the element in the lower-right corner of the matrix
M_3[___, ___] <- 1
# Print M_3
print(___)