個別鵯記

【第1399回】 東京出張前にゼロ割対応版拡張マクネマー検定の関数を書いてみた（2016年9月13日）

• 前回書いた拡張マクネマーの問題だが，全部のセルに0.01を足すのは，検定したい行列オブジェクトをxとして，mcnemar.test(x+0.01)で済むので，分母がゼロになる組合せをスキップする最低限の関数を書いてみた。使い方はmct(x)とするだけ。

  mct <- function(x) {
  L <- NROW(x)
  X <- 0
  N <-0
  for (i in 1:(L-1)) {
  for (j in (i+1):L) {
  s <- (x[i,j]+x[j,i])
  if (s>0) {
  N <- N+1
  X <- X + (x[i,j]-x[j,i])^2/s }}}
  return(list(X2=X, df=N, p=(1-pchisq(X, N)))) }

• 今日は公衆衛生学会の編集委員会のため10:50神戸空港発の飛行機で羽田に飛ばねばならないが，その前にいったん出勤して事務仕事をしていくため，2番のバスに乗らねばならない。6:00起床し，ご飯に刻みショウガと牛そぼろを振りかけた主食に，おかずは鶏ムネ肉とオリーブの実と冷凍野菜をオリーブオイルで炒めて，味醂とスウィートチリとアンデスの岩塩と妻から貰った東欧のミックススパイスで味付けしたものを作った。想像以上に美味。
• いったん名谷キャンパスに出勤して事務書類を提出してから神戸空港へ。飛行機の機材到着の遅れということで25分ほど離陸が遅くなったが，とりあえず問題はない。12:30過ぎに羽田に着き，ローソンでお握りを買って昼飯を済ませてから新宿へ。15:00から公衆衛生学会編集委員会。いろいろ議論はあったが無事に終わった。来月も大阪で行われる年次総会に併せて拡大編集委員会が開催されるので，9月から11月は毎月編集委員会があるという状況。
• 編集委員会の後は，新宿御苑前から赤坂見附で銀座線に乗り換え，新橋から都営浅草線・京急直通に乗って羽田空港へ。18:00頃着いたので，まずSkyMarkの自動発券機でチェックインした。最終の神戸空港行きは満席のようだ。晩飯でも食べてから手に持つ検査場を通ろうと思っているところ。
• 羽田カレーを食べてから手荷物検査場を通り，ほぼ予定通りの時刻にボーディングし，予定通りに離陸。最新ネット犯罪のリアルを感じさせ，テクニカルなアイディアは凄いような気がするが，それ以外は実にラノベっぽいラノベな怪作『天才ハッカー安部響子と五分間の相棒』に続く一田和樹の新作『女子高生ハッカー鈴木沙穂梨と0.02ミリの冒険』を機内とポートライナーとバスの中で読了。この作家さんは（少なくともこのシリーズは）こういう作風なんだな。これら２作品は同じ世界線の続きの話なので，前作を好きだった人ははまると思うが，決して万人向けではなく，例えば川端君の『The S.O.U.P.』とか，幸森軍也『あなたの待つ場所』とか，『なれる！SE』シリーズみたいな意味でのテクニカルな記述はされてないし，たぶんそんな説明はしたくないのだと思われた。

Tweet

△Read/Write COMMENTS

▼前【1398】（家で洗濯しながらレポートを見る（2016年9月11日-12日） ）　▲次【1400】（事務仕事が終わらない（2016年9月14-15日） ）　●Top

Notice to cite or link here | [TOP PAGE]