今回は前回紹介しなかったコマンドなどについて，紹介する。しかしあまり気合が入っていないのでグダグダな感じになっている。 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 # コロン３つなの…
ちょっと気になったので，近年一部界隈で話題になっているベイズ情報仮説の評価をやってみる。 普通に計算すると大変なベイズファクターも，情報仮説の元では簡単に計算できる。元ネタは以下の資料 http://www.tandfonline.com/doi/abs/10.1080/17405629.201…
尊敬するHoxo_mさんらが翻訳された”みんなのR”にあるGelman's secret weaponを真似して自分の武器を作成してみようかと思う。 d.hatena.ne.jp https://www.amazon.co.jp/みんなのR-データ分析と統計解析の新しい教科書-Jared-P-Lander-ebook/dp/B012Z86Q6Qや…
Rから今流行のスクリプト言語Juliaを使いたい。 普通にJuliaにパスを通してRのsystem関数を使えばいいのだけれど，パスを通さずにRで下のような関数を書いて作業をする。jlfileにはJuliコードを書いた拡張子.jlであるファイル名を入れる。jlfilepathは作業フ…
このノートの荒っぽいまとめ。 https://www.statmodel.com/examples/webnotes/webnote14.pdf３ステップで潜在クラスの数を決める。 最適な対数尤度探し optseedとtech 11の挿入 optseedとtech 14の挿入 ブートストラップ尤度比検定を行うには結構時間がかか…
Mplusはミクスチャーモデルの分析もらくらくできるし，計算も結構速い。しかし，人間業が深いもので，もっと早く分析して欲しい，なんていう欲望が出てくる。 そんな時はamalysisコマンドの中にprocessors = 4;などと記述すると，複数コアで計算を回してくれ…
MCMCサンプラーのうちの1つであるJAGSをRから利用する新しいパッケージができたようだ。基本はrjagsのラッパーのようだが，並列計算の指定などがやりやすそうだ。今後，これを中心に使っていくかもしれないのでヘルプにあった例をコピペしておく。初期値の指…
学会に行く意味とはなんだろうか。今回はこのことに関して自分にとって何が良いことなのか考えてみる。 まず，第1に締め切りがあることによる効果が挙げられる。締め切りがあることによって，適当なところで研究をまとめなければならない。コレによって，1つ…
いきなり新しい趣味を始めたくなった。昔から見るのが好きだったジャグリングをやろうと思う。一番基本的なボールというやつから手をつける。3つのボールを小気味良くポコポコと投げ続けられるようになりたいものだ。 ネットでやり方を調べてみたり，youtube…
自分で文章を書くときにやりがちなことをまとめておく。文章を書くのは苦手なので，注意点をキチンと書いてみる。この文章の目的は，自分にとっての現在把握している気をつけるべき点を明確にし，書き残すことで自分の作文能力を向上させるための資料を作成…
Mplusは潜在変数モデリングに適した解析ソフトで，心理学，社会学，疫学などの分野で近年人気が出ているソフトウェアである。ホームページは次のものがある。 Muthén & Muthén, Mplus Home Pageこのソフトは非常に便利で，文法も簡潔なのであるが，文法が特…
3レベルの階層線形を行う時のlmerのメモ学校のなかに，個人がネストされていて，個人の中に測定時点がネストされている場合。学校IDをschool_id, 個人IDがind_id，あとは独立変数をx, 従属変数をyとする。簡単のため，中心化は行われているとして，レベル2以…
今日やったこと 某ポスター作成 発表資料の作成 MCMCの観察 なす購入 1.2.については直前に再度確認を行うこと。3. について，バーンインは結構早い。サンプリングがえらく時間がかかる。いらないパラメータはサンプリングしないでいたほうがいいか？ 欲しい…
最近読んだ良記事 最近のRの界隈で流行っている（気がする）purrrパッケージの使いかたについて。 よく参照させてもらうので，とりあえずまとめておく。sinhrks.hatenablog.com sinhrks.hatenablog.com uribo.hatenablog.com {purrr} による非テーブルデータ…
前回から随分間が空いているが，気にしない。 思ったらやってみるのが大事な気がする。最近の困り 研究のオリジナリティが良くわからない。何を研究したら面白いのか全くわからない。 査読のコメントへの対応が困った。オリジナリティの主張の仕方，伝え方。…
ノートパソコンで作業をしていると，そのノートパソコンに付属しているキーボードをそのまま使いがちです。 それはそれでいいのですが，外付けのキーボードを使ってみるとコレが意外といいもので，癖になります。 キーボードは文章を書くにしても，何かプロ…
何事も続けるのが大事ということですが，淡々と継続し続けるのも難しいものですよね。 本当にいまやっていることに飽きないで続けられればいいのですが，私は飽きっぽいので何かと続けることができません。 大事なことを続けるために大事かもしれないと最近…
日記を書き始めて二日目をすでにすっ飛ばして三日目に突入しています。 しかも，三日目もすでに終わりが近づいてきているという。。。 さて，今日は日記に何を書けばいいのかということについて，少し考えてみました。 私は昔からなにかと文章を書くのが苦手…
新年あけましておめでとうございます。 2015年という新しい年に合わせて，新しいことを始めてみようかと思いブログなんてものを始めてみました。 いままで興味はあったものの，どうやっていいのか余りわかっていなかったのですが一念発起して行動しようかと…
