【第615回】 年明け初ミーティング（2021年1月7日）

• 8:00起床。80.95 kg，97％，35.8℃。外は雨。
• 高山先生が諮問メンバーになっているなら，寄付先はここがいいかな。
• ふと思い立ってRStudioのリリースノートのページを見て気づいたが，バージョン1.2の2019年11月1日版から花の名前のコードネームがついているんだな。最新版はApricot Nasturtiumなので，杏色のキンレンカと思われる（WikipediaによるとNasturtiumは本来オランダガラシを指す名前だったそうだが，杏色の花ならばキンレンカで良いのだろう）。R本体のコードネームがPeanuts!由来なのも洒落ているが，RStudioの花の名前も粋だな。
• Zoomで打ち合わせとか来学期の講義テキストの準備をしていたら17:00になったので帰途に就く。年明け初ミーティングは自宅からの予定。
• 東京の新規感染報告数が2400人を超えたが，今日までほぼ何も対策しなかったのだから，それはこうなるよな。
• 適当に晩飯を済ませてからミーティング。たぶん早く文章を書き始めた方が良いので，そうアドバイスしておいた。
• 欠損値処理の質問に答えるため，mvnmleパッケージについて調べてみた結果わかったこと（なお，無料で読める日本語による欠損値処理の説明としては，覚え書きとのことだが，村山航 (2011) 欠損データ分析(missing data analysis) ―完全情報最尤推定法と多重代入法―がわかりやすいと思う。英語だとcranのMissing DataのTask Viewやfinalfitというパッケージのvignetteがコンパクトでわかりやすい）。mvnmleは2012年頃にlme4の開発者であるDouglas Batesの協力を得てKevin Grossという方が開発しメンテナンスしていたが，既にメンテナンスを止めていて，cranではORPHANEDという扱いで，去年5月に削除済み。しかしインストールはGithubから可能で，devtoolsパッケージがインストール済みなら，

  library(devtools)
  install_github("cran/mvnmle")

  とすればインストールできる。ちなみにKevin Grossは2015年からテニュア付きの教授職を得ているが，webサイトにはmvnmleのことはカケラも載っていないのでもうメンテする気はないのだろう。LittleMCAR()についてはこの説明によると，BaylorEdPsychパッケージにも入っているようだが，このパッケージも既にcranから削除されている。LittleのMCAR検定については，メンテされていないパッケージを使うより，本当は，Little RJA (1988) A test of missing completely at random for multivariate data with missing values. Journal of the American Statistical Association, 83: 1198-1202.を読んで自分でコーディングした方が良いと思う。Stataには実装されているので，Rに現行のものがないとは思えないのだが……と思ってTask Viewを見直してみたら，RBtestパッケージでMCAR性の検定ができるそうだが，典拠となる文献が書かれていないし，回帰に基づく方法ということでLittleの方法とは違うようだ。なお，質問者が陥ったというエラーは，たぶんここに書かれている例みたいに，すべて欠損な変数が入っているのではないかと推察している。

Tweet(list)

▼前【614】（遠隔会議2つ，いや3つ（2021年1月6日） ）　▲次【616】（講義準備（2021年1月8日） ）　●Top

Notice to cite or link here | [TOP PAGE]