2023-06-01から1ヶ月間の記事一覧

ガンマ関数(再訪)

6月25日(日)にShibuya.lispのもくもく会がありました。14:30にオンラインで集合して、自己紹介と質問タイム(勉強になった!)があり、15時頃~17時頃にもくもくと作業をしていきました。(もくもく会の後の雑談は予定を1時間近くオーバーして盛り上がりまし…

本日の給油(CRF1100Lアフリカツイン/2BL-SD10)

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 2860 12.74 160 18.21 8.78 2022-07-31 3128 13.63 157 19.66 7.99 2022-08-27 3401 16.63 157 16.42 9.56 2022-11-03 372…

applyとreduceの違い

リスト内の数値の合計を計算しようとして (defun sum (lst) (apply #'+ lst)) (sum '(1 2 3 4 5)) ;=> 15 というような関数を書きました。普通に使うぶんには全く問題なかったのですが、長いリスト(具体的には480,000個のdouble-float)を引数にしたら「ス…

二項検定と二種類の過誤

先日、二項分布の累積分布関数を計算できるようになりました。 marui.hatenablog.com これを使って、簡単な二項検定をやってみたいと思います(あまりLISPとは関係ない内容ですが)。題材としては、ジャンケンにやたら強い人がいて「俺には天から与えられた…

二項分布の累積分布関数(と正則化不完全ベータ関数)

統計に使える関数の実装を続けています。2月のLISPもくもく会では二項分布の確率質量関数(dbinom)を実装しました。今度は累積分布関数(pbinom)が欲しいところです。 marui.hatenablog.com 原理どおりに動くのであれば確率質量関数を累積加算すると累積分…