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

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

千葉でのキャンプのときに大荷物を乗せてのバイク走行が大変だったので、積載の方法を考え直さないといけないなーと思っているうちに、福島でのキャンプ予定が決まってしまいました。せっかくの大型連休ですしね。 走行の当日は雨になる予報だったので、東京…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-04-10 2283 14.23 163 19.40 8.40 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 340…

日本のAmazonで「Julia」と検索すると（「18歳以上ですか？」に「いいえ」と答えても）セクシーな画像が多く並んでしまうので、これから購入したい本のメモも兼ねて一覧を作りました。とは言っても出版された件数はまだまだ少ないので作成は楽ちん。リンクは…

youtu.be 給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-04-10 2283 14.23 163 19.40 8.40 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-…

以下のエントリに、二項分布の確率質量関数をCommon Lispで計算する関数を紹介したのですが、どうも意図した通りに動いてくれないときがありました。少し時間を作って調べてみたら、コードに間違いがあるのが明らかになりました。 marui.hatenablog.com Load…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39 3.43 2022-07-27 13405.0 2.21 157.92 57.29…

ごくたまにオーディオシステムの準備をすることがありますが、パワーアンプからスピーカへの結線は赤黒ケーブルが逆になってしまったり、左右スピーカへの接続を逆にしてしまったりということがあります。そこで、ステレオスピーカへの結線がうまくいってい…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-23 12910.1 2.55 158.04 37.57 4.21 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39…

【2023-04-17】このページの最下部に載っているコードでは、特定の引数の時に間違った答えが出力されたり、無限ループに入ってしまうことが分かりました。たとえば(dbinom 4 10 0.5) ;=>0.557となってしまいます（正しくは0.205）。元のアルゴリズムには問題…

車輪の再発明を続けています。今日は画像の書き出しです。 スペクトログラムを計算できるようになったものの、それを画面表示するためにはなんらかの画像にしないといけないんですね。画面にウィンドウを作って表示するよりもファイルにするほうが簡単なので…

Common Lispで振幅スペクトルの計算ができるようになり、Hann窓もかけられるようになったので、次はスペクトログラムの実装にいきたいと思います。 スペクトログラムは、振幅あるいはパワーの強さを時間と周波数ごとに計算したものです。以下の画像ではAudac…

macOSにはafで始まる名前の、オーディオファイルを操作するコマンドラインアプリが入っています（afclip、afconvert、afhash、afida、afinfo、afplayの6つ）。僕がよく使うのはafinfoとafconvertです。 afinfoは以下のようにオーディオファイルの情報を表示…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-23 12910.1 2.55 158.04 37.57 4.21 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-05 1740 18.22 158 18.77 8.42 2022-03-07 2007 16.46 162 16.22 9.99 2022-04-10 2283 14.23 163 19.40 8.40 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 286…

非負整数nの階乗n!はn×(n-1)×(n-2)×(n-3)×...×2×1と定義されています。また、0!=1とします。例えば6!であれば、6×5×4×3×2×1=720です。これを、整数以外にも拡張したのがガンマ関数で、nが整数の時にはΓ(n+1) = n!となります。例えばΓ(6)=120、Γ(7)=720となり…

目次 正規分布の確率密度関数 累積分布関数の近似式（Abramowitz and Stegun） 累積分布関数の逆関数の近似式（Abramowitz and Stegun） 累積分布関数と逆関数を誤差関数から計算する（Numerical Recipes） Common Lispを音だけでなく、統計にも使いたいです…

youtu.be キャンプ道具を揃えたシリーズの番外編です。本編は以下の3つ。 キャンプ道具を揃えた2022（バッグと小物類） キャンプ道具を揃えた2022（調理器具） キャンプ道具を揃えた2022（テーブルとチェア） 最初のキャンプ道具を揃えた2022（バッグと小物…