幾霜

　あー俺こういう話だめなんだわ。分かっていても耐えられない。

　最終回が何故か録画されていない・・・orz。よりによって2時間スペシャルを失敗なんて欝だ。早く再放送してくれませんかねぇ。

　温度変わりすぎなんだよバカヤロー。今日はあったかいのに・・・。

　MacバイナリをリコンパイルしてTigerで動作するようにした。蓑さんでcodemlによる尤度計算をサポートした。デフォルト設定でInvariantを無効にした。デフォルト設定で塩基・アミノ酸頻度パラメータの最適化を無効にした(データの平均値に固定)。蓑さんでcodemlで尤度計算する場合にSite-Specific RateやN-GAMやAutocorrelated Discrete Gammaをサポートした。とりあえずダウンロードできるようにはしましたが、解説PDFやWebページ中の説明はまだ更新していません。

　-isysrootオプションを付けるだけでLeopardでは通っていたのですが、Snow Leopardでは色々ファイルが無いと言われて困っていました。また、コンパイルできたバイナリも実際にはLeopard上で走らなかったりしました。が、-mmacosx-version-minオプションを併用することでうまくいくようです。具体的には以下のようにすることでTigerで実行可能なバイナリが作成できました。

-isysroot /Developer/SDKs/MacOSX10.4u.sdk -mmacosx-version-min=10.4

　Leopard以降で動作するバイナリを吐かせるには以下。

-isysroot /Developer/SDKs/MacOSX10.5.sdk -mmacosx-version-min=10.5

　なお、gcc-4.2ではTiger用バイナリが作成できないことが多いようですがgcc-4.0を使えばうまくいきました。

　予報では朝には止んでいるはずだったので外に出てみましたが、自転車で帰るのはまだきつそうでした。というわけで居室へ逆戻り。

　名古屋市立大と福岡女子大(いずれもリンク先はPDFです)で教授・准教授・講師いずれかの公募が出ていますね。名古屋市立大はパーマネントっぽい。分野は進化学・遺伝学・生態学などみたいです。いずれも31日まで。

　福岡女子大の「英語で授業できること」っていうのは厳しいなぁ。準備にものすごい時間と労力をかけられるのならできなくはないでしょうけどねぇ。

　さすがにポスターそろそろやばいなと思い始めた。ので今から作り始めます。

　4年生が卒研発表お疲れ様会で余った食物をしこたまタッパに詰め込んで持ってきてくれた。えらい! NMZさんの日頃の教育の賜である。

　6時半、みぞれが降る強風の中計算科学まで行ってミーティング。しかし、うちのボス(棟続きの部屋の住人)しか来ていなかった・・・。何のために我々はあそこへ行ったんだ。

　昨日練習していた3人は午前の発表。ちょーっと直ってない部分とかもありましたが、滞りなく終了。3人とも進学なので研究は続けるわけで、頑張っていただきたいものです。

　先日購入したUSB 3.0接続のHDDが、1週間ぶりのバックアップをしようとして電源入れてもPCに認識されませんでした。ケーブルかケースかHDDのどれかまたは複数が壊れたのでしょう。カードは別の外付けHDDを接続して問題なくアクセスできたので大丈夫でしょう。

　まずはケーブルを交換することでケーブルをテスト(USB 2.0ケーブルしかありませんが)。次にHDDを入れ替えることでケースをテスト。最後にHDDを別の方法で接続することでHDDをテスト。の予定。

　PAUP*やTNTのTBRは、内分枝を一つ選んで抹消することで系統樹を2つの部分木に分割し、その2つの部分木から枝を1つずつ選んで接続するというもの。TBR探索では、これらが全ての内分枝・部分木の全ての枝の組み合わせで候補樹形が作成され検討される。

　MrBayesのExtending-TBRでは、内分枝を無作為に1つ選び長さを変え、一方の部分木との接点をその接点のある枝に隣接する枝に無作為に移し、反対の接点でも同じことを行う。

　以上のように、MrBayesのExtending-TBRは明らかに通常のTBRよりも樹形探索範囲が狭い。というのも通常のTBRでは接点は隣接する枝に限らず部分木上の全ての枝に移すことを検討するからだ。また、Extending-TBRは1回のproposal中に、実際には「内分枝の伸縮」、「一方の接点を隣に移す」、「もう一方の接点を隣に移す」の3つのproposalを含んでいるように思われる。ということは、Extending-TBRでは通常のMetropolis-Hastings法よりもacceptのハードルを下げる必要があるのではないだろうか。MrBayesがそんなことをやっているのかはわからない。だがやっていないのなら過度にrejectされやすくなり結果として事後確率は最大の樹形でより大きな値となると考えられる。

　やっぱりMac版バイナリのテストにはMacが要るよなーということでMac miniを検討。スペックからしたら高いのですが、まぁ6万円なら許容範囲か。でもほぼこのテスト専用機になってしまうのに、そんなことに6万円を費やすというのもばかげている気がする。うーどうしよう。いっそのこと中古品でも買うか。

　フランス・スペインに卒業旅行に行っていたM2が帰筑。お土産のアメを口に入れたらみな顔をしかめることに。「香水の味」としか表現しようのない品でした。ものすごく不味かった・・・。

　午前の発表練習で以前からあまりに良くなりすぎていたのが一人いて驚いた。それなりに時間があったとは言えあれだけやれれば上出来。

　朝食用の卵を忘れていて慌てて買いに行ってきた。そしたらなんかめっちゃ外が寒かった。三寒四温とは言うもののこのところの変動の激しさはさすがに体に堪える。だからと言って暖かくなればその日は花粉?で激しい頭痛や肩こりに悩まされるという・・・。痛し痒し。どっちもつらい季節やな。桜咲く頃にはそれも治まるのですが、まだ20日くらいはかかるかな。気象協会の予報では25〜31日に開花とのこと。

　先月のアクセスは1日平均11,108ヒット、3,606ページビュー、1,303ユニークユーザーでした。まぁ検索エンジンのクローラーが多いので、仮に半分(さすがにこの見積もりは多すぎると思いますが)がクローラーだとすると約650人の人間が訪問している計算になりますね。その半分以上が研究者か学生でしょうからざっと300人はその手の人が見ているわけですか。もっとも、これは幾霜と第五次元の合計値なので、均等割すれば幾霜を見ているのは150人ということになりますか。

　結局、徒歩で買い物へ。まずは腹ごしらえということで大成軒で麻婆麺650円。それから歩いて駅前へ。Q'tの花水木でドイツの人へのお祝いの品を探すもあまりめぼしいものは無し(そう言えば他の二人にも何か用意した方が良いかなぁ)。ライスシャワーというサクランボと何かのフレーバーティーをいただきましたが、ものっすごい安物の茶葉の味にげっそりしました。香りだけを楽しむお茶ですね。

　それから西武5階の無印良品でヘアターバンをしこたま買い込んだ。かなり昔から愛用しているんですが、けっこう劣化が早いのが難点。今回は細い方を3つと太い方を1つ。あと、化粧水と乳液もよさげなものがあったので購入。小分け用のスプレーとボトルも一緒に買いました。乾燥肌なので毎年冬は悩んでいますが、今回はオーガニックなんとかいうやつにしました。レジでお金を払う際に値段に気付いて真っ青になったのは秘密。たっか。でももう冬終わりなんですよねー。

　次に西武1階のザ・ガーデンでペルーシュ・ブラウンとネクタフローの蜂蜜を購入。ペルーシュは良い明日に行けばでかい箱の割安なやつが売っているのですが、今の砂糖で入れたミルクティーがいまいちおいしくないのに我慢できず買ってしまいました。雨さえ降っていなければ・・・。蜂蜜は朝食用。もうすぐ前のやつが切れそうなので。今回は野草の蜂蜜にしました。

　最後にてらしまで牛乳と入浴剤とブルーレットドボン。

　また買い物に行けないじゃないの。

追記 - 16:10:51
　まだパラパラと降っている。ダメだこりゃ。うーん、歩いて駅前に行くしかないかなぁ。ちょっと荷物持って移動するには距離がある気がしますが・・・。

　今日はうまく入れられた。うまい。しかし袋に書いてある「冷たいミルクを2、3滴加えてどうぞ」というのが意味分からん。そんなちょっとだけ入れるのに意味あるのか? 今はちょうど牛乳を切らしてしまってできないが、今度やってみるか。

　先日ぎょーむ日誌に書かれていて鼻血噴いていたのですが、現物が出てきて大・興・奮!!

　BUGS言語をC言語に翻訳する機能があり、CPUネイティブなコードでMCMCを走らせることができるらしい。また、多系列対応で、系列数だけCPUを使えるとのこと。

　早速伊藤さんがお試しになった結果を掲載されています。今後が楽しみです。

　理由はそれかい!

　特命係がいなくてもトリオ・ザ・捜一が正解には到達していたでしょうね今回は。逃亡に間に合ったかどうかは別として。

　「捜査なんてのは、推理だとか何とか、そんなもんじゃありません。」

　この時点では、ここが笑いどころだってことが分かりませんでした。最後まで見てからこのシーンに戻るとかなり笑えます。語るねぇ伊丹さん。

　外に出かける気が失せてしまった・・・。買い物に行きたかったんですけどねぇ。やっぱり車は必要ですね。

　非常に新しい分岐の年代を推定する際の最大の問題点は、非常に古い較正点からあまりに離れていて精度があまり期待できないということでしょう。しかし、この論文では13〜11万年前の化石から得られたミトコンドリア全ゲノム配列データを較正点とすることでこの問題をある程度解決しています。なるほどその手があったかという感じです。実際に使うのはサンプルがないとどうしようもないのでかなり難しいでしょうが、昨年のメダカの分岐年代推定結果などをこの方法で検証できると面白いですね。

　そもそもハマダラカが臭いを感知していることも知りませんでした。生物の持っているセンサーのメカニズムって面白いなぁ。しかも役に立ちそうだ。ハマダラカの場合は攪乱できればマラリア被害を防げるし、人工的に作成したりコンピュータに繋げられれば人感センサーにも応用できそう。