通常の逆行列に代わる方法
直前の動画では、一般化逆行列(疑似逆行列)について簡単に触れました。
この演習では、MASS パッケージにある Moore–Penrose の一般化逆行列を使います。扱っている行列が可逆であれば、結果が通常の逆行列と一致することを確認しましょう。
ginv() コマンドは、R で Moore–Penrose の一般化逆行列を計算します。
Moore–Penrose 逆行列の詳細は、次の リンク を参照してください。
この演習はコースの一部です
Rで学ぶデータサイエンスのための線形代数
演習の手順
- WNBA の Massey 行列
Mと得点差ベクトルfは読み込まれています。もう一度Mを表示してください。 Mの通常の逆行列をfに適用し、結果を表示してください。MASSパッケージはすでに読み込まれています。一般化逆行列ginv(M)をfに適用し、結果を表示してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Print M
print(M)
# Find the rating vector the conventional way
r <- ___%*%f
colnames(r) <- "Rating"
print(r)
# Find the rating vector using ginv
r <- ___ ___ ___
colnames(r) <- "Rating"
print(r)