概要 新しいマシンにwin11を新規インストールした時に遭遇したエラーとその対処を以下に書きます 簡単に解法だけ先に書くと，win11を新規に入れようとして install driver to show hard wear という画面になって先に進まなくなったら， そのマザボに対応する…
授業で使ったベイズファクタ計算用のコードのメモです*1． 必要なパッケージのロード library(rstan) library(bridgesampling) パッケージがなければ， install.packages("bridgesampling") でインストールします． stanでのコンパイルの効率化のために，以…
作業メモ 必要なこと Heroku アカウント登録（アプリをアップするためのサーバ） Git CMDのインストール（Git を動かすためのシェル） Heroku CLIのインストール（Heroku関連コマンドをGit CMDで使うためのツール） otree hubのアカウント登録（必須ではない…
やること Claudeを使って，全自動でotreeトリートメントつくってもらう． 例によって，忘れるので自分用のメモです． 準備 MacでJSをインストールしてからターミナルでclaudeをインストールする こちらの記事を参照しました なおサファリでQiitaが読めないと…
よく分からない注 パール『入門統計的因果推論』でDAG上のノードの従属性について説明している箇所に，「$Y$が$X$の影響を受けるのに，$X$と$Y$が独立になる」という意味の注があって，しばらく意味が分からないまま放置していた． 別の文献で忠実性について…
やりたいこと pdfファイルから複数ページにわたって特定位置にあるテキストを抽出したい pdfファイルからのデータ抽出 方法 PythonのPyMuPDFパッケージを使う． 手順 PyMuPDFをpipでインストールする． 指定pdf内の指定矩形領域からテキストを抽出する関数を…
はじめに たとえば次のような分析結果の図を作りたいと考えたとしよう． アウトプットの例 このとき「実際の計算」と，「図の作成」は別々にコード化すると効率が良い． ひな形の作成 上記の図のイメージができたら，「その結果を生成するようなダミーのdata…
この記事の内容 運と才能のシミュレーション(Pluchino et al 2018)を数理モデルとして一般化した論文 Hamada, Hiroshi, 2023, Luck of Outcome in the Talent Versus Luck Model, Advances in Complex Systems Vol. 26, No. 04n05. をもとにjuliaでパラメー…
三日後に忘れるのでメモを残しておく 前準備 VSCodeインストール LaTeX-workshop extensionインストール TeXliveインストール Github copilotインストール 設定 こちらのサイトを参考にして.latexmkファイルを作成し，ルートフォルダに置く． Winの場合，C:\…
ラウンド毎にパラメータをランダム化する方法 最初はよく分からないので def set_parameters(group: Group):#グループオブジェクトを引数にした関数，引数groupのタイプ（クラス）がGroupであることを指定 rvec = C.Rvec.copy() random.shuffle(rvec) group.…
冬場になると車の補機バッテリーの警告灯が毎朝点滅する． 自粛要請時代に一度バッテリーが上がってしまったので，残量を確かめるために市販の車用バッテリーチェッカーを購入した． 以下その方法のメモ 前準備 うちの車は補機バッテリーがトランク側に設置…
特定のアプリケーションのエラーでwindwosが強制終了し場合に，再起動すると「起動ドライブ内にOSが見つからない」というエラーがしばしば発生する． 以下，その対処方法． MSIのロゴが出たらf11キーを連打してUFEIを起動する． 起動ドライブの優先順序が表…
久しぶりにTeXLiIveをアップデートした． 3日後に忘れるのでメモを残しておく． TeXStudioの設定 設定→コマンド→latexを platex.exe -synctex=1 -interaction=nonstopmode %.tex に変更する．この変更で日本語文書クラスのコンパイルができる． beamerでminj…
前回の記事でワンショットの剥奪ゲームのトリートメントを作成しました． 同じ被験者にゲームの設定を変えつつ複数回繰り返す方法を考えてみましたが3日後には忘れそうなので，例によって メモを残しておきます． グループに共通の値を設定する 結論から言う…
はじめに とある原稿を書いているうちに，昔作ったztreeトリートメントをアップデートしたい気持ちがふつふつと高まってきたので，otreeで新たに実験用コードを作成しました．以下忘れた時用のメモです． otree準備 アナコンダの最新版をインストール pip in…
Win10からUbuntu経由でPythonを使おうとしたら思いのほか手間取ったので作業メモを残しておく． 問題：WindowsPowerShellからWSLをインストールできない 原因：AMDのCPUが仮想環境を使える状態になっていない 解決法：BIOS→アドバンス設定→OC→CPUの機能→SVM …
OSがインストールされたM.2SSDドライブを，異なるマザーボードにCPUと一緒に移植するとOSはそのまま起動できるのでしょうか？ またチップセットとCPUが異なる場合はどうでしょうか？ 実験した結果を記録を残しておきます． 実験1 （CPUは同じだがチップセッ…
すぐ忘れるのでメモ． library("ggplot2") y1 <- rnorm(10,5,1) y2 <- rnorm(10,2,1) y3 <- rnorm(10,3,1) years <- 2001:2010 data <-data.frame(years,y1,y2,y3) ggplot(data, aes(years)) + geom_line(aes(y = y1, colour = "y1")) + geom_line(aes(y = y…
確率変数の組の独立性 とあるエコノメのテキストに確率変数の組 $$(X_i,Y_i), i=1,2,...,n$$ が独立ならば $$X_j, j=1,2,...,n$$も独立である（大意），と書かれていた． まあようするに確率変数のペア$(X_i,Y_i)$と$(X_j,Y_j)$が独立ならば，ペアの片方だけ…
はじめに 介入の割り当て確率が共変量に依存して異なるとき，統制群と介入群の応答変数の平均の差によって介入効果を推定するとバイアスが生じます． もし，共変量の値毎に割り当て確率を決定する関数が分かっていると，推定量のバイアスを取り除くことがで…
はじめに 2018年12月に『その問題，数理モデルが解決します』をベレ出版さんから刊行しました． 当blogではこの本の内容紹介をしていなかったので，以下にざっと解説します*1． その問題、数理モデルが解決します作者:浜田宏ベレ出版Amazon 本書は数学，統計…
実験2 前回の続きです． データを生成する分布として次のような複合分布を考えます． $$ X \sim \mathrm{Bernoulli}(q) $$ $$ q \sim \mathrm{Beta}(a, b) $$ これはベルヌーイ分布のパラメータ$q$の事前分布がベータ分布$\mathrm{Beta}(a, b)$にしたがう複…
イントロ 以前，Mathlogに自由エネルギーについて書きました． mathlog.info 自由エネルギーは，データへのモデルのフィットを確認するために便利な指標の一つです． 特に，階層モデルでWAICの計算が煩雑な場合に，相対的に近似計算が簡単な自由エネルギーは…
久しぶりにSourcetree（GitをGUI操作するためのローカルクライアント）を経由してGithub上のファイルを更新しようとしたら，「Githubの認証方式が変わったのでアクセスできません」というエラーが発生しました． qiita.com こちらを参考にトークンでの認証に…
イントロ この記事はベイズ予測分布を応答変数データのヒストグラムに重ねてplotする方法のメモです． 目的 応答変数が離散的なスケールで測定されている場合，bayesplotのppc_dens_overlay関数が使えなくて面倒くさいってことがあります． 無理矢理やると，…
イントロ 微分方程式をベースにした統計モデルに関して興味深い論文を読んだので紹介します． Naik, Prasad A., Murali K. Mantrala and Alan G. Sawyer, 1998, Planning Media Schedules in the Presence of Dynamic Advertising Quality, Marketing Scienc…
PCがうるさい なんだかPCがうるさい．HWMonitorで確認したらケースファンの回転数が常時1000RPMを超えていました． マザーボードの設定を変えたら静かになったのですが，多分3日後にはやり方を忘れてしまうので記録を残しておきます． 4ピンではなく3ピンの…
ファイル名を直書きするのがめんどうくさい 50個のデータセットを使い複数のモデルで推定して自由エネルギーを計算する必要があったのでコードを書いてみました． 多分忘れるだろうから，ここに貼り付けておきます． fevalues.m0 <- c() for (i in 1:50) { s…
はじめに 『その問題，やっぱり数理モデルが解決します』をベレ出版さんから刊行しました．以下，その内容をざっと紹介します． その問題、やっぱり数理モデルが解決します: データ時代を生き抜くための数理モデル入門作者:宏, 浜田発売日: 2020/09/14メディ…
内容 ジニ係数をPythonを使って計算します．よく知られた定義ではなく，簡略化された定義式を使うと計算が速くなる，という例を紹介します． ジニ係数の直感的意味 ジニ係数は，所得や資産などの非負の資源分配の不平等度を測定する指標です． たとえば，社…
