研究を行う上で非常に重要なのは，やはりどのような問題に決着をつけるべきなのかその問題を見極めることだと思う。問題を設定するために，はじめは指導教員や先輩の研究者などに教えを請うのが望ましい。それによって，何が研究としてあるべき形なのかを身…

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枚のグラフかのように見せかけるのがポイント。 けっこう可愛い感じのグラフを描画できる。 # # 縦向きの，信頼区間の図。 # ポイントは，レイアウト…

最近思うこととしては、とにかく毎日進み続けることが何よりも大切。 ルーチンワークでやらなくてはならない雑務も多く、体力や気力も削られるかもしれないが、その中でも研究をし続け、論文をなんとしても書き続ける必要がある。 書き続けていれば、書くこ…

今回は前回紹介しなかったコマンドなどについて，紹介する。しかしあまり気合が入っていないのでグダグダな感じになっている。 SCALE OF DEPENDENT VARIABLES FOR ANALYSIS OPTIONS FOR DATA ANALYSIS CLASSES AUXILIARY SURVIVAL VARIABLES WITH SPECIAL FU…

「MONTECARLO:」コマンドを使った場合に一緒に利用するコマンド・オプションの一部を一覧にしておく。 GENERAL SPECIFICATIONS DATA GENERATION CUTPOINTS GENCLASSES NCSIZES CSIZES HAZARDC MISSING DATA GENERATION 次の記事で紹介する予定？ GENERAL SPE…

Mplusを使ったモンテカルロ・シミューレーションを実行する際に使えるコマンドを，マニュアルに沿って紹介する。 Missing 単純なランダムな欠測 MAR カテゴリカル潜在変数 マルチレベル 探索的因子分析 実際のデータの利用 マニュアルの例では生存時間解析の…

潜在変数モデルを統合的に扱うことができる統計分析ソフトであるMplus(https://www.statmodel.com/)はかなり広範な分析を実行可能である。 基本的な使い方については，インターネット上に多くの例が示されている。 例えば，以下の清水先生のサイトなどに日本…

適切な思考を展開するために、いくつか必要となる要素が存在する。それは、言葉の定義、前提、根拠だと思われる。 まず、言葉の定義を明確にすること。これにより、用語の境界を明確にできる。言葉を定義しないと、主題のズレや言葉の用法のズレに気がつくこ…

以下のページから引用qiita.com ポイントは次の部分： ”Google日本語入力のプロパティを開き、「一般」タブの「キー設定」に「キー設定の洗濯」という項目があると思います。ここで「MS-IME」を洗濯した上で横の「編集」ボタンを押します。この中の「モード:…

1ページあたりの設定について。 1行の文字数，1ページあたりの行数の競って地には，形式→文書のレイアウトから設定すること。後ほど追記すること。スタイルは使わないとシヌ。

潜在クラス分析(有限混合モデル）やら，k-meansとかいろいろと初期値に依存する分析するがある。また，Latent variableモデリングではEMアルゴリズムを使った最尤推定を行うことが多いが，複雑なモデルになるほど解が収束しないこともしばしばある。 とくに…

「研究とは何か」ということについて。研究とはどんなゲームか。 分野にとって解く事ができる，価値のある重要な問題に効率的に解を与えるゲーム。 人を納得させればよい ゴールの見極めが大事。 報告はフォーマットに則って行う。 ルールをキチンと理解する…

Rのリスト操作について改めてまとめる。 特に，purrrなどを用いた操作を理解すること。purrrの使い方についてのページ purrr: ループ処理やapply系関数の決定版 - Heavy Watal Hadley神のListについてのページ Lists · R for Data Scienceforループ的なlappl…

一番有名な，関数名をタイプするものでは見れないものについて。 あまり日本語の記事が無いようなきがしたので，メモを残してておく。 library("Matrix") Matrix # Matrix関数の中にある謎の関数。このままだとエラーが返ってくる。 spV2M # コロン３つなの…