個別鵯記

【第913回】 再び採点三昧（2015年2月5日）

• 6:20起床。昨夜は沸かしたのに風呂に入らなかったので，沸かし直して入浴。朝食は，昨夜セブンイレブンで買って帰った生野菜にバルサミコ酢と塩コショウを掛けて食べ，炊き立てのご飯に冷凍しておいた明太子載せという主食。
• 往路3番のバスで出勤。今日も採点三昧，時々院生の相談への対応，という予定だったが，メールなどで事務仕事がいくつか届き，まずはその対応に追われた。
• 昨日も書いたハドラーさん絡みの話だが，@kohsukeさんのtweetを見て，dplyrを使わないでやってみたらどうなるかと思ってやってしまったので，採録しておく。要は，自分が使っているRの環境で，インストールされている全パッケージのうち，作者にHadley Wickhamさんを含んでいるパッケージの割合を計算し，ハドリー指数(Hadley Index)と呼ぼうという提案であった。それくらいハドラーさんはR世界を席巻しているということだな。

  ■RTした@kohskeさん作ハドリー指数をdplyrの無い環境で
  mean(grepl("Hadley", sapply(rownames(installed.packages()), packageDescription, fields="Author")))
  0.065だった

  ---

  ■（承前）関数化してみる
  PAIDX <- function(X="Hadley") { mean(grepl(X, sapply(rownames(installed.packages()), packageDescription, fields="Author"))) }

  ---

  ■（承前）今作った関数を使ってみる sapply(c("Hadley", "Nakazawa", "Fox", "Kanda", "Hyndman", "Ripley", "Ihaka", "Dalgaard", "Hornik", "Murdoch"), PAIDX)

  ---

  ■（承前）xに付値してcat(paste(names(x),round(x,3)),"\n")
  Hadley 0.065 Nakazawa 0.009 Fox 0.028 Kanda 0.005 Hyndman 0.037 Ripley 0.074 Ihaka 0.009 後略

  ---

  ■（承前）確かに自分が作ったパッケージは2つだけで、全部で230余りのパッケージをインストールしてあるから、0.009という数字は当然なのだが、なんだか悔しい。まあ、Rの始祖Ihakaさんと同じ数字だからいいか。

  ちなみに，Hadley Wickhamさん自身は，0.235という，物凄い値を叩きだしている。
• 永田靖: 多重比較法の実際, 応用統計学 27: 93-108, 1998.は，統計的多重比較法の適用において実務家から寄せられた疑問点を材料に，「多重比較法に関する誤解・誤用」，「多重比較法を用いる際の注意点」などについて検討した論文。対照群が２つある場合の多重比較についての詳しい議論や，多重比較を考慮したサンプルサイズの設計などは大変参考になる。1998年なので，FDR法には言及がなかった。
• 昼過ぎにアリムさんが帰国する前にもう1本仕上げておきたいという論文のためにsemの相談を受けた。どうしてもfitを良くするには試行錯誤が必要だが，そのときの方針の参考として去年出たレビュー論文を読んで貰うことにした。役に立ちそうな内容だったが，当時はmod.indices()だった関数が既にdepricatedで，現在はmodIndices()を使うので，そこは修正が必要だった。
• その後，某学会の編集後記を書く仕事を済ませ，後は採点三昧だったが，やっぱり終わらないまま20:00を過ぎ，直通終バスで帰宅した。冷凍ハンバーグを解凍してから野菜と一緒に炒め，焼き肉のたれを掛けたものと，白い御飯で晩飯。簡単にできる割には美味だった。
• GM1がもはや4万円を切ることはなさそうなので，ダブルズームキットで34,400円だったGF5WAを注文してしまった。

Tweet

△Read/Write COMMENTS

▼前【912】（会議三昧（2015年2月4日） ）　▲次【914】（今日も採点と学生の相談に乗るのと事務仕事の日（2015年2月6日） ）　●Top

Notice to cite or link here | [TOP PAGE]