以下の記事を読んで個人的に思ったことをまとめておく。 hidekatsu-izuno.hatenablog.com 個人的な思いは，タイトルの通りだけどそれ以外にも色々思ったことがある。この文章はなにかの戦いに決着を付けるためのものでもないし，現状の交通整理をするための…
Windowsへのgfortranのインストール www.compe.media.kyoto-u.ac.jp makeの仕方 www.matsumoto.nuem.nagoya-u.ac.jp Lapacの導入 qiita.comFortran www.str.ce.akita-u.ac.jp
テラモナギさんの記事と，PRMLの11.4節を参考にしてみた。 https://teramonagi.hatenablog.com/entry/20130203/1359888495#fn-b0191759 本当はもとの論文を見るべきだが，大体あっていそうなのでOKなような気がする。 （1変数の確率密度からの）スライスサン…
この雛形を使うと並列計算ができる。個人的には返り値などは基本的に全部書きだしてあとで分析する形にしておくと安心。 シミュレーションを関数にしてしまうことと，foreach関数の指定のときに必要なパッケージや環境の変数（自作の関数とか）を全部持って…
また読んでいても気になる論文があったので，備忘録的にツッコミをメモしておく。 http://team1mile.com/sjpr61-1/shimizu.pdf 第2章のベイズ統計学の流れの整理区分は後の議論を明確にする上でもわかりやすいように思う。さて，第2章2節から本題。 統計モデ…
仕事の仕方について自分なりやっていること工夫など自分のためにメモを残しておく。言語化することによって，より精緻化したり，不要な部分を削除したり工夫ができるようになると思われる。論文を書く手順 アイデアを思いつく 先行研究を調べる（このときに…
参考： https://mc-stan.org/docs/2_18/stan-users-guide/multivariate-hierarchical-priors-section.html 1.13 Multivariate Priors for Hierarchical Modelsより 高速なコードにはなっていないが，これで基本的にはOK。 添字が慣れたものではないが，読め…
下のリンクの論文を読んで，ベイズ統計について，少し気になる記述が多かったので，気になった点や考えたことをメモしておく。 https://www.jstage.jst.go.jp/article/psychono/37/2/37_37.25/_pdfベイズ統計学については東京工業大学の渡辺澄夫先生のWebペ…
F列は全体のfactor score，C_Fはmost likelyなクラスに対するfactor score。 あと，fscoreはEAP。http://www.statmodel.com/discussion/messages/13/741.html Rick Sawatzky posted on Tuesday, June 03, 2008 - 10:10 am I am fitting an polytomous unidim…
お金関係で便利なサイトhttp://disclosure.edinet-fsa.go.jp/http://tskl.toushin.or.jp/FdsWeb/view/FDST000000.seamhttp://csr-toshokan.net/index.phphttp://www.jsda.or.jp/http://www.jsda.or.jp/shiru/index.html冊子わかりやすい http://www.jsda.or.…
GitHubにある，crayonパッケージを使えば簡単 github.com 詳しくはヘルプを見たり，exampleを走らせればすぐにわかる。 簡単に使えるので，良さそう。 library(crayon) cat(blue("Hello", "world!\n")) # Crayon defines the %+% string concatenation opera…
フィールドについてなど，このページに書いてあることが非常に詳しい。 もう少しきちんとWordの使い方について勉強する必要がある。山葵のアプリのサプリ
ショートカットキーをいじるといろいろと都合が良い。 → よく使うやつは，登録しておくこと。慣れてきたら増やすこと。 （Windowsのやつとの整合性をつけること。覚えることが多くなるのは嫌だけど。。。）Word ショートカットメモ command + cont + 1 黒 (C…
短期的目標(今年中) 今年中にあと1，2本の論文を投稿する。可能であれば英語で。 新しい武器を身につける。 中期的目標(向こう5年) 英語論文で10本論文を出版する。 コンスタントに英語論文を投稿できる基盤とリズムをみにつける。 新しい方法論について、さ…
例によってダラダラと気になることをなんの脈絡もなく書き留めておく。昨今のベイズ統計学を推進しようとしている人たちがやっていることは，ベイズ統計学というより，実際にはパラメタ推定方法として，事後分布を使うということに過ぎない。 実際には予測分…
最近，ベイズ統計学が流行っていることに関して，個人的に思ったことをとくに引用などもせずにただつらつらと書いておく。学問的な正しさは一切保証しない。また，私はベイズ統計学や頻度論的統計学の区別なく，統計学を愛しているし，不当な争いや煽りあい…
みんなのRcpp teuder.github.io Rcpp 使い方 - Rcpp 入門アルマジロの使い方 aaaazzzz036.hatenablog.comaaaazzzz036.hatenablog.comaaaazzzz036.hatenablog.comaaaazzzz036.hatenablog.com
全般的な入門 rishida.hatenablog.comhttps://nbviewer.jupyter.org/github/bicycle1885/Julia-Tutorial/blob/master/Julia高速チュートリアル.ipynbThe Fast Track to Juliawww.geocities.jpLinear algebra · The Julia Language 行列の使い方 行列 | Julia…
以下のページがすばらしい。 今までの表を作る形式とは何だったんだろうか。https://std-experimental-optional.blogspot.jp/2017/10/word-word-2016.html
問題を解く前に考えるべきこと。第１に何を目的にしているのかを明確にしなければならない。大きな目標からしたら，現在目の前にある”問題のようなもの”は本当の問題ではないかもしれない。実際には取り組んだり解決する必要がないものの可能性がある。もっ…
研究を行う上で非常に重要なのは，やはりどのような問題に決着をつけるべきなのかその問題を見極めることだと思う。問題を設定するために，はじめは指導教員や先輩の研究者などに教えを請うのが望ましい。それによって，何が研究としてあるべき形なのかを身…
IRTモデルで予測分布や対数尤度の計算などをおこなってみる。 対数密度の記述は以下の頁の質問に少し書かれている。 https://sourceforge.net/p/mcmc-jags/discussion/610037/thread/c728e87b/ The general rule is that if y ~ dfoo(alpha, beta, gamma, ..…
やり方は簡単で，予測分布を未知の変数としてモデルに書けばいい。 Stanではgenerated quantitiyブロックでrng系の関数を使えばいいが，JAGSではそのまま"x_pred ~ dnorm(mu, sigma)"などとするとRでいうところのx_pred sigma)のように乱数を生成してくれる…
WAICなどを計算するには，対数尤度のMCMCサンプルが必要。 個々のMCMCサンプルからがんばって計算する方法もあるが，より簡単にはMCMCサンプラーのなかで対数尤度を得ると楽。 Stanではgenerated quantityが容易されているが，JAGSでは少し工夫が必要。参考 …
気が向いたらStanやJAGSなどとの対応関係をきちんと取りならが検算を行っていきたい。 MCMCサンプルは「.gh5」ファイルに入っているので，これを読み込む。以下の資料のp.65あたりの「mplus.get.bayesian.」を利用すると，ギブスサンプラーのMCMCサンプルが…
文字（数字）の桁数を揃えたい。 参考URL： http://d.hatena.ne.jp/teramonagi/20100902/1283429021 formatC("hoge", widh = 5, flag="0")flagは '0 + - #' だけが指定可能とのこと。 不要な空白を消したい。 # いらない空白を全て削除する # 全角半角の指定…
Analysisコマンドで，rotation = targetを指定し， modelコマンド内で目標とする値を~で指定する。 また，各因子の最後の(*t)などとつけて，各因子負荷が1つの因子負荷行列で表すことができることを示さないといけない。 Mplusのマニュアルは，versin 8のp68…
Wordで便利な小技を見つけた。 どちらも，「高度な検索と置換」を使う。 特に，ワイルドカードと，検索対象を指定するのがポイント。 「[０-９]」で検索を行うと，大文字の数字が選択される。その上で，「A->a」で「半角」を選択すればOK。テキストボックス…
Rで，信頼区間を書いていいプロットするコードのメモ。 layout関数で，上手くプット画面を分割して，さも1枚のグラフかのように見せかけるのがポイント。 けっこう可愛い感じのグラフを描画できる。 # # 縦向きの，信頼区間の図。 # ポイントは，レイアウト…
最近思うこととしては、とにかく毎日進み続けることが何よりも大切。 ルーチンワークでやらなくてはならない雑務も多く、体力や気力も削られるかもしれないが、その中でも研究をし続け、論文をなんとしても書き続ける必要がある。 書き続けていれば、書くこ…
