Replace `np.sum(a * b)` with `a @ b` in lecture-python.myst

[oyamad]

やり方

1. まだ作っていない人は，https://github.com/QuantEcon/lecture-python.myst の "フォーク" を作る ("Fork" ボタンを押す)．
2. 自分のフォークで，編集したいファイル (たとえば lectures/mccall_model.md) を開く．
3. 鉛筆マーク ("Edit this file") を押して，編集モードにして，ファイルを編集する．
4. 終わったら "Commit changes..." ボタンを押す．
5. 適切に "Commit message" を書く．(QuantEcon/lecture-python.myst の過去の Pull requests を参考にするとよい．)
6. 重要："Create a new branch for this commit and start a pull request" の方を選んで，新しい branch 名を書く (デフォルトのものでもよい)．
7. Pull request を送る．(ここは，最初は教員といっしょにやる．)

[UdohNakamura]

担当範囲

中村

• lectures/career.md: self.F_mean = np.sum(self.θ * self.F_probs)

• lectures/career.md: self.G_mean = np.sum(self.ϵ * self.G_probs)

• lectures/lake_model.md: β * γ * np.sum(np.maximum(U, V) * p_vec)

• lectures/lake_model.md: w = np.sum(V * p_vec * (w_vec - τ > w_bar)) / np.sum(p_vec * (w_vec -

• lectures/mccall_model_with_separation.md: d_new = np.sum(np.maximum(v, u(c) + β * d) * q)

• lectures/mccall_model.md: reject = c + β * np.sum(v * q)

• lectures/mccall_model.md: return (1 - β) * (c + β * np.sum(v * q))

• lectures/mccall_model.md: h = np.sum(w * q) / (1 - β)

• lectures/mccall_model.md: h_next = c + β * np.sum(s * q)

須田さん

• lectures/mccall_q.md: reject = c + β * np.sum(v * q)

• lectures/mle.md: return np.sum(y * np.log(μ) - μ - np.log(factorial(y)))

• lectures/mle.md: return np.sum(y * np.log(μ) + (1 - y) * np.log(1 - μ))

• lectures/mle.md: return np.sum((X.T * y * ϕ / μ - X.T * (1 - y) * ϕ / (1 - μ)),

• lectures/multivariate_normal.md:Σθ_hat_arr_C = np.array([np.sum(C[n, i+1:n+1] ** 2) for i in range(n)])

• lectures/opt_transport.md:total_cost = np.sum(X * C)

• lectures/prob_matrix.md:mean = np.sum(Benford_pmf * k)

• lectures/prob_matrix.md:var = np.sum([(k-mean)**2 * Benford_pmf])

[oyamad]

今日やった lectures/mccall_model.md は機械的に置きかえるだけどよかったけれど，必ずしもそうとは限らないので，内容をチェックしてください．

たとえば，lectures/opt_transport.md:total_cost = np.sum(X * C) の X と C は2次元配列なので，X @ C とすると別のものになります (サイズが合わないのでエラーになるはず)．

[suda-yuga]

以前送ってくださったやり方通りにやったところ、自分に対してpull request しただけになってしまい、少し違う方法でやることになってしまったのですが、試しに一個やって見たので問題がないか確認してくださると助かります。

[oyamad]

コード修正問題ないです．よくできました．

• こっち Replace np.sum(a * b) with a @ b in lecture-python.myst.mccall_q.md QuantEcon/lecture-python.myst#472 はそのまま生かして
• こっち Replace np.sum(a * b) with a @ b in lecture-python.myst.mccall_q.md suda-yuga/lecture-python.myst#1 はクローズしましょう ("Close pull request" を押す)．

[suda-yuga]

ありがとうございます。
度々申し訳ないのですが、lectures/mle.mdに関して
return np.sum(y * np.log(μ) - μ - np.log(factorial(y)))　は
return (y @ np.log(μ)) - np.sum(μ) - np.sum(np.log(factorial(y)))とすれば
書き換えることができなくないですが、そうした方がいいのでしょうか？chatGPTは計算速度は変わらないし、見にくいというのですが。

[oyamad]

@suda-yuga lectures/mle.md の np.sum(y * np.log(μ) - μ - np.log(factorial(y))) はそのままにしておいた方がよいですね．

• それより np.log(factorial(y)) は問題ですね．Overflow のリスクがある．ChatGPT に聞いたら scipy.special.gammaln(y + 1) にすべしとのこと．
• TeX の式にタイポがある．\log y! は \log y_i! としないといけない (2か所)．こっちを pull request してくれてもよいです．

[suda-yuga]

@oyamad
ありがとうございました。該当箇所を訂正しました。

[oyamad]

@suda-yuga mle の PR にコメントしました．

• Add gammaln to from scipy.special import factorial
• scipy.special.gammaln -> gammaln

だけやってもらって (ファイルを編集して，現在のブランチ suda-yuga-patch-2 に新たなコミットとして加える)，あとは中の人のコメントを待ちましょうか．

[suda-yuga]

@oyamad
承知しました。ありがとうございます。上記の修正含め、私の担当箇所は一通りpull requestしました。

[oyamad]

lecture-python-advanced.myst にもいくつかあるので (検索結果)，こっちもやっていただいてけっこうです．