ぎょーむ日誌 2012-02-(01-10)
2012 年 02 月 01 日 (水)
-
0600 起床．
昨晩やらせておいたアリ計算は
……
うーむ，
いまいちだ．
どう改善したらいいだろう．
-
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
寒波接近中．
0820 研究室着．
-
アリアリ．
-
アリにやられていようが，
生態学会さーばー雑用あれこれはとびこんでくる．
けっきょく，
私の「調査」はへなちょこで，
Alivenet のヒトたちがかなり有望そうな要因を発見してくれたわけだが
……
サーバーの状況を確認させていただきました。
drwebプロセスについては常時比較的消費メモリも高い状態ではありますが、
当該時間帯に特別なものは見受けられませんでした。
一方でバックアップの処理のようなものは設定されているようで
そのための負荷が高いように見受けられます。
AM10:30頃にtopコマンド等で状況を確認すると、CPU稼働率も高く
gzip, perl, tar が上位にくる状況でございましたが、
現在は落ち着いており、そのコマンド行は確認できません。
Pleskバージョンが異なる為、Plesk側で設定されている内容については
詳細を確認する事が出来ませんでしたが、
/var/lib/psa/dumps/
以下には10時30分から40分の間に作成されたた「.tgz」ファイルが
100以上確認出来ます。
たしかにめちゃくちゃたまってますね．
この毎日 10:20 に起動する莫大ばっくあっぷ作業，
設定したのは私ではないねえ
……
というか，
やはり細かいところまでいちいち確認しないとダメだ．
意外や意外．
-
昼飯．
-
調査地である公園数は 81．
あまりにもデータの少ない 4 公園は除去しよう
……
といろいろ調べてるうちに，
データ整形でちょっとしたまちがいを見つけた．
やれやれ．
-
ひたすら試行錯誤．
-
1820 研究室発．
買いもの．
1850 帰宅．
晩飯の準備のてつだい．
晩飯．
-
今夜も，
夜中までアリモデル試行錯誤．
計算まち時間は台所仕事とか．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0640):
米 0.5 合．
ダイコン・ニンジン・ナメコのスープ．
- 昼 (1200):
研究室．
米 0.7 合．
大金精肉ソーセイジ．
レンコン・ニンジン・豚肉の炒めもの．
加賀揚げ．
ミニトマト．
レタス．
リンゴ．
- 晩 (1900):
キャベツ・ホタテ・豚肉などのカレー鍋．
雑炊．
2012 年 02 月 02 日 (木)
-
0600 起床．
昨晩やらせておいた統計モデルの MCMC サンプリング
……
やはりダメか．
端的にいって，
データの中に情報が足りないわりに，
モデルが複雑すぎるということだろう．
-
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
まだ晴．
モデルを簡単化する方法をひたすら考えながら歩く
……
こういう詐術の発明みたいなことは，
職業的うそつきの本領が問われるところだ．
0820 研究室着．
-
今回のアリ動態もでる (空間構造・時間構造あり)
であつかう要素は以下のごとしだ．
以下を同時にあつかっていた．
-
各調査地におけるアリ採集データから，
その場所が低・中・高密度のいずれの状態に属するかを決める
(同時に各状態の「アリ採集確率」を探索しつつ)
-
状態の遷移を決める
(同時に遷移確率を決めながら)
-
各遷移確率に影響を与える要因の効果を調べる
これらをいっぺんにやってしまうのが理想なのだが
……
ここ数日の試行錯誤で「ぜんぶいっぺん」はムリらしい，
とわかってきたので
……
1. と 2. & 3. を分離しよう．
今まで作ったモデルとその推定計算プログラムを分割すればよい．
-
まず状態遷移だけを決める部分．
-
そしてここが詐術なのだが，
各場所・各時刻におけるアリ密度の状態を MAP でおきかえてしまう，
という小細工をする
……
今まで推定計算がうまくいかなかったのはどうもこのあたりにあるようで，
とくに欠測をウメているところなんかはもーしっちゃかめっちゃかで
……
-
かくのごとく「状態」が固定されてしまうと，
各場所・各時刻における遷移も決まってしまうので，
遷移確率のモデリングもできますということで．
-
予想どーりといいますか，
いろいろ局所的な環境要因の影響はよくわからん，
ということで
……
まあ，
これはしょうがないでしょう．
-
モデルを分割することで，
「とりあえず」の結果は得られたし，
「このデータはどういうものか」
という見とおしも良くなった．
計算時間は全部やっても 5 分弱
(WinBUGS による MCMC サンプリングを並列化した場合)
……
いままでひとつの推定で何時間も計算させていたのは，
いったい何だったのだろうか
……
-
寝不足ぎみでチカラつきてきた．
昼飯．
[窓の外の雪はかなり豪雪っぽい]
明日は朝から移動，
香川県に行って，
このアリモデルの発表をしないといけないんだけど，
どうなることやら．
鉄道・航空の運行は雪に阻害されないだろうか，
発表準備は気力ぎれとともに終息してしまうのではなかろうか
……
-
精神的よゆーがないので，
今回は LibreOffice
の Impress で発表準備．
R でどんどん PNG file を生成し，
Impress の中にとりこんでゆく
……
-
ぜんぜんオワってないけど，
帰宅時刻だ．
1730 研究室発．
雪．
もー勝手に降ってください．
買いもの．
1800 帰宅．
晩飯の準備．
晩飯．
-
その後も発表準備はつづく
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.5 合．
キャベツ・油揚の味噌汁．
カブの葉・油揚の炒めもの．
- 昼 (1230):
研究室．
米 0.7 合．
カブの葉・油揚の炒めもの．
レンコン・ニンジン・豚肉の炒めもの．
加賀揚げ．
サツマイモ．
ミニトマト．
リンゴ．
- 晩 (1900):
米 0.5 合．
キャベツ・油揚の味噌汁．
コマツナ・ハタケシメジ・卵・豚肉の炒めもの．
ゴボウの煮物．
2012 年 02 月 03 日 (金)
-
0600 起床．
今日は香川県までアリ出張．
新千歳空港のほうはそんなにふらなかったし，
今日はだいじょうぶだろう
……
しかし，
なんだか不安になって新千歳空港のサイトで東京行き
AIRDO 便の状況を調べてみた
……
ほーら，
どの便も無事にでてるぢゃん．
ただし私が搭乗予定の 0900 発 14 便は欠航
-
あとからわかったんだけど，
雪のせいではなく機材の故障だそーで．
-
えーい，
時刻までに到達するためにはひとつ早い便にするしかない．
9 時ではなく 8 時発になるのか．
いまから出れば新千歳空港に到達可能．
-
AIRDO のサイトで航空券変更を試みるが
……
なぜかうまくいかず．
時間をムダにした．
AIRDO 予約センターの電話番号を調べる．
-
0625 緊急発進的に自宅発．
雪の上を桑園方向にかけながら AIRDO 予約センターに電話．
この時間なのになかなかつながらない．
なんとか変更できた．
0635 JR 桑園着．
-
0642 快速エアポートにのって同発．
けっこう混んでる．
しかも札幌で降りるヒトは少なくて，
みなさん空港までのってる．
車内でぐったり．
-
0725 新千歳空港着．
AIRDO カウンターでクレジットカードだして，
変更便の二次元バーコードうけとる．
こういう場合は，
もちろん追加料金とか払わなくてもよい．
-
0750 AIRDO 12 便 18F に搭乗．
B-767-300．
予定よりちょっとおくれて離陸．
よーやく朝飯．
-
昨晩は 25 時すぎまでかかって発表準備してたんだけど，
もうちょっと修正したほうがよいんだけど
……
どーもヤル気がでなくて，
機内では作業しなかった．
-
0945 羽田空港着．
ひとつ早い便できたので，
ここでの待ち時間は 1.5 時間ぐらい．
いわゆる「パソコンデスク」で作業．
-
おや?
羽田に HANEDA-FREE-WIFI
なる新しい無料 WiFi ができたよーだ．
ThinkPad とかでも問題なくつながる．
暗号化されてない通信だけど．
-
発表内容の修正作業，
ちょっとだけ進捗．
-
「高松は雪なので」着陸できない場合は，
岡山か伊丹に降ります……ってなんだよー，
と思いつつ 1120 ANA 535 に搭乗．
また B-767-300．
機中で準備のつづき．
-
まあとくに問題もなく，
1245 高松空港着．
なるほど，
空港周辺には雪が残っている．
[恵方まきデザート]
瀬戸内の魚のさまざまな刺身 (とくにサワラ)，
牡蠣，
親鳥
という歯ごたえのある鶏肉のニンニクいため，
豚肉の角煮，
ローストビーフとか，
食べきれないほどおいしいものがいろいろと．
写真の白いふきんのように見えるものは，
マスターの遊びココロのあるデザート「恵方まき」，
イチゴを綿菓子でくるんだもの．
-
小林さんとひさしぶりにいろいと話したり，
JaSPa
というアイデアでこの地域を見ているといったことを教えていただいたり．
-
宿は JR 高松のちかく．
-
[今日の運動]
-
うんどうしてない
……
朝 JR 桑園まで走ったとか?
-
[今日の食卓]
- 朝 (0640):
羽田にむかう機内でサンドイッチ．
- 昼 (1200):
握り飯．
- 晩 (1950):
ごちそう．
2012 年 02 月 04 日 (土)
-
高松出張二日目．
0500 起床してしまった．
寝不足で眠いのに眠れん
……
ぬるい風呂に入ってみたけどだめ．
7 時ごろ，
無料朝食券をもって宿に併設されてる焼き鳥屋で
朝飯．
土曜日は「鶏汁」だそうで．
このへんの人たちは鶏肉をよく食べるのかな．
-
8 時ごろ宿発．
コスト削減というか，
私のとまってる「別館」みたいなところは完全に無人化されてて
(法律のぬけあな?)，
カギは 1F の箱にぽいっと返却するだけ．
駅前・朝食つきで一泊 5130 円だけど．
[JR 高松駅前]
宿の出口からみた光景．
なんだか不思議な建築物・オブジェクトがいろいろとあるんだけど
……
他のひへーした地方都市とかに比べると，
高松市はこういうところにお金かけるよゆーがある，
というかんぢもする．
-
今日はなぜか気まぐれで，
高松から神戸まで高速バスで移動し，
神戸空港から新千歳空港に飛ぶ．
-
で上の写真ののっぽビルの足もとにバスターミナルがある．
あらかじめネットで予約・支払いをすませていたので，
カウンターで何かをやることなし．
予約メイルの印刷出力を乗車・降車時に示せばよいとのことだが
……
印刷してないね．
iPhone とかの画面を示せば OK とのこと．
0815 高松発．
[高速バスの座席]
高速バスにのるのはすごくひさしぶりだ．
各座席に電源コンセントあり．
車内はけっこう満員だった．
女性客おおい．
神戸で買いもの?
北海道でいうと道北・道東から土日に札幌に買いもの日帰り旅行するようなかんぢ?
[大鳴門橋]
香川県から淡路島にわたる橋．
私にとっては淡路島初上陸．
[あわじしま]
休憩のため停車した緑 PA にかかげられていた地図．
高速道路は島の中央部を走っているので，
島というより，
なんだか山の中にいるようなかんぢで．
[西日本 JR バス]
高松-神戸 3600 円．
3 時間弱．
[淡路島北端]
明石大橋が見えてきた．
[明石海峡大橋]
航路高が 65 m もあるので，
橋の北側にある神戸市垂水区を見おろしつつ近づいていく．
-
1110 ごろ三宮バスターミナル着．
1120 ポートライナーに乗って三宮発．
1138 神戸空港着．
[たこ焼き]
480 円．
明石焼きではなくたこ焼きうってた．
手早くぱっぱっと食べようと思ってたんだけど，
関西風で外はぱりっとしてるけど中はクリーム状，
熱くて食べるのにけっこう時間かかった．
香川でうどん食べれなかったので，
空港内にあった讃岐うどん店にもかなりひかれたのだが
……
ここでうどんを食べるのはなんだか敗北感があったのでヤメ．
[B-737-700]
本日の搭乗機．
ANA ブルーではなく黄色っぽい塗装だったので，
不思議に思って
調べてみると
ゴールド塗装の記念機か何かみたいで (ダッシュ 700 の最初のやつ??)．
B-737 ってのも驚くべき長命機種だね
……
ダッシュ 700 ともなると，
機体・エンジン・アヴィオニクスともどもまったく別ものみたいだけど．
1245 定刻どおり出発．
-
後記:
「外はかりっ」たこやきは関西風ではなく「築地銀だこ」ふうではないかという
ご指摘をいただいたのですが
……
調べてみると
銀だこのたこ焼きでは，
最後に表皮に油を塗布してフライするので「かりっ」となるとのこと．
神戸空港ではそれはやらなかったんではないかな
……
ということで，
おそらく銀だこほど「かりっ」としてません．
-
1435 新千歳空港着．
また寒い雪国にもどってきました．
別世界だな．
1449 JR 快速エアポートにのって同発．
1525 札幌着．
のりかえて，
1535 桑園着．
AEON で買いものして，
1610 帰宅．
-
コーヒーのんでひとやすみ．
-
晩飯の準備のてつだい．
晩飯．
-
ぐったり．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
宿の朝飯．
- 昼 (1150):
神戸空港のたこ焼き．
- 晩 (1950):
バゲット．
カブ・ゴボウ・ニンジン・セロリ葉・大金ソーセイジのスープ．
鶏もも肉のタンドリー．
サツマイモ．
2012 年 02 月 05 日 (日)
-
0730 起床．
よくねた．
香川アリ発表ばてからかなり回復してきた．
ラジオ体操．
朝飯．
紅茶．
-
怠業．
1/7-2/3 までの電気使用量 122 kWh．
2869 円．
-
1140 自宅発
ミニ大通公園
氷点下雪上走．
1220 帰宅．
体重 73.6 kg．
香川アリうんどう不足だ
……
-
かくのごとく「体重が落ちない状態」であるにもかかわらず，
なぜかしら昼飯は外で．
-
ゆっくりとした昼飯からもどる．
また怠業．
-
先日の自分のアリ発表をみなおす．
他のヒトの発表メモもよみなおしてみたり
……
-
昨年，
廿日市市のいろいろな公園で (アルゼンチンアリのみつかっていない公園も含めて)
ベイトトラップをつかった一斉駆除みたいなことが実施されたと伊藤さんが言ってたが
……
廿日市市のサイトをみると，
たしかにそういう記録
があるな．
「2011年9月 廿日市地区町内会連合会が実施するアルゼンチンアリの一斉防除に併せ公共施設の防除に取り組む」
か．
-
などと高松アリ研究会について検討しているうちに，
金曜日の夜に「総合討論」やった瓦町ちかくの店が気になってきた
……
そう，
ばててて店名確認とかしなかったところ．
-
地図上で町名・建物名を確認しつつ，
あれこれ検索してみると
……
うーむ，
どうやら
「おきる」
だったのかな．
位置的にはここらへんのはず
……
あ，
画像検索
をみたら正解だと確信できた．
お，
YouTube に TV 出演時の動画ファイルもある
……
新感覚！和食とサイエンスの融合．
店主・小島さん，
やっぱり研究者的な料理人だったんだ
……
(金曜日にいただいた料理の工夫あれこれをふりかえりつつ)．
-
まあ，
これからアルゼンチンアリ研究をがんばって，
来年も高松の研究会に行こうかな
……
-
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0800):
バゲット．
チーズ．
カブ・ゴボウ・ニンジン・セロリ葉・大金ソーセイジのスープ．
カブ・セロリのマリネ．
- 昼 (1200):
知事公館南のラ・トルテュでエゾシカとかいろいろおいしいものを．
- 晩 (2000):
米 0.5 合，
明石のタコめし．
ネギ・ハタケシメジ・油揚の味噌汁．
カブ・セロリのマリネ．
セロリ・ニンジン・豚肉のソース炒め．
2012 年 02 月 06 日 (月)
-
0620 起床．
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
今日はすでに気温がプラス 2 度．
ひさしぶりの「脱・真冬日」．
0820 研究室着．
-
気象庁公表の過去データを調べてみると，
先月
は 23 日が最後の「脱・真冬日」だった
……
ということで，
今日も真冬日だったら「二週間連続真冬日」となっていたわけだが，
脱・真冬日になりましたということ．
予報では，
明後日からまた真冬日がつづく．
-
どうも，
先週までのどたばたもえつき，
ぼけぼけな状態．
-
アリ情報を香川と共有するための wiki を何か新しく
……
などと考えていろいろ調べてみた．
dokuwiki
とか twiki
なんかがおもしろそうだなぁ，
などと．
-
しかし，
結局のところ，
以前に vikiwiki で以前につくったアリ wiki があることを
「再発見」してしまったので，
そちらを惰性で流用することに．
いやはや．
-
昼飯．
-
午後は project 管理ソフトウェアの試行錯誤にはまっていた
……
これは，
院生たちの野外調査スケジュール管理に使えないかと考えて．
-
Redmine
の設定に苦闘してみたり
……
-
1800 研究室発．
買いもの．
1830 帰宅．
晩飯の準備．
晩飯．
-
工学部の人から「この台車はそちらのものでは?」
という写真が送られてきた
……
そうなんですよね，
この台車は，2009 年 7 月から 2010 年 3 月までの
環境科学院耐震工事のせいで工学部 Q 棟・R 棟に疎開してたときに行方不明になったやつで
……
翌日いただいたメイルによると，
工学部 R 棟にある段ボール箱集積室に長期間放置されてたそーで．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.5 合．
ネギ・ハタケシメジ・油揚の味噌汁．
カブ・セロリのマリネ．
- 昼 (1230):
研究室．
食パン．
セロリ・ニンジン・豚肉のソース炒め．
リンゴ．
- 晩 (1930):
米 0.5 合．
ニンジン・タマネギ・豚肉のスープ．
鶏肉のシュプリーム，
インドなスパイス版．
カブ・セロリのマリネ．
2012 年 02 月 07 日 (火)
-
0620 起床．
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
曇．
今日も気温高い．
0820 研究室着．
-
今朝もプロジェクト管理ソフトウェアをちょっと調べてみたり．
dotproject
は PHP なのでラクだ．
-
いまどきのプロジェクト管理ソフトウェアなるものがどんなものか，
少しばかりわかってきたような
……
おもしろいソフトウェア (web application software)
だけど，
院生たちの野外調査「でます」「もどりました」管理に使えるようにするには，
改造がいろいろ必要だし，それは複雑すぎる．
この程度の web application software
ぐらいなら自分で作ったほうがマシなように思えてきた．
cookie の勉強とかすればよいのかな
……
まあ，
これは 3 月後半ごろの作業になりそう．
-
窓の外はかなりの降雪．
二日間「脱」真冬日だったけれど，
これからまた寒い日々がつづくもよう．
-
昼飯．
[そして台車回収]
昨日かいてたように，
工学部で行方不明の台車が発見されたので，
雪の降るなか回収の旅にでた．
-
廿日市アリデータのもうひとつの側面
--- アリ群集モデルについて考えてみた．
時系列データなので，
案外簡単なものになりそうな気もする．
もちろんほとんどの種で「他種との相互作用があるのかどうか，
よくわからない」という結果がでるのは必然だが
……
伊藤さんとか「アリ名人」になってしまうと，
5 分間というかぎられたサンプリング時間内に，
その異常なる瞬間識別能力で
「これはさっき採ったから，こっちだ」
と標本種数最大化をやってしまうので，
珍しいアリがどんどんとれてしまうんだよね．
-
まあ，
種ごとの効果とかはでないだろうけど，
「自分 vs 自分以外連合」
の対立ぐらいが見つかればよいかな
……
-
14 時より植物コースの修論直前な発表練習．
まあ，
けっこうやばいところが多々あるわけだが
(データを「圧縮」して 4 点にして 3 パラメーターな曲線をフィットするとか，
説明変数に対数をつけ忘れているとか)，
まあ限られた時間での修正をがんばってもらいたいものです
……
-
またアリ群集モデルについて検討．
先日のアルゼンチンアリ拡大もちょっと再検討する必要あるのだが
……
-
1720 研究室発．
買いもの．
1750 帰宅．
晩飯の準備．
晩飯．
-
弁当の準備．
また野菜をきざんだり，
イモをふかしたり．
マメ科の樹木．
西アフリカでは料理の材料 (スパイス) をくだく臼として，
この樹種を使うそーで．
写真にうつってるやつはレプリカみたいなもので，
実物はおそらくもっと巨大なのだろう
……
-
[今日の運動]
-
[今日の食卓]
- 朝 (0650):
米 0.5 合．
ニンジン・タマネギ・豚肉のスープ．
モズク酢・カニカマボコ．
- 昼 (1200):
研究室．
米 0.7 合．
鶏肉のシュプリーム，
インドなスパイス版．
セロリ・ニンジン・豚肉のソース炒め．
ゴボ天．
プチトマト．
ミズナ．
リンゴ．
- 晩 (1950):
米 0.5 合．
ネギ・シジミの味噌汁．
キャベツ・ニラ・ブナシメジ・豚肉の炒めもの．
シシトウの炒めもの．
2012 年 02 月 08 日 (水)
-
0620 起床．
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0800 自宅発．
また真冬日の日々がつづく．
0830 研究室着．
-
他にもやらねばならないことあるのだが
……
先日の廿日市アルゼンチンアリの解析をみなおす必要があるので，
まずはそのあたりをうろうろ．
-
10 時に神戸大の佐倉さんがきてくださったので，
いろいろ相談．
来年 2 月 (一年後ですな) に理学部で集中講義することになったので．
生物学科の各学年の人数は 25-30 人ぐらい．
学科の演習室があり PC が使える．
講義は 15 時間で 2-2.5 日ぐらいの日程．
その他あれこれ
……
-
昼飯．
-
いやはや，
先日の香川出張の航空券の領収書のうちひとつをとりわすれていた
……
1730 までに AIRDO 営業とやらに出頭しなければならぬ仕儀に．
-
アルゼンチンアリ解析みなおし作業，
やっぱり公園番号の齟齬があった．
システム全体で統一する．
いまさらながら．
-
統計モデルあてはめもやりなおし．
まず第一段階，
各公園のアリ密度状態の時間変化．
左がもとデータの割算値で，
右は推定．
一致してない点もあるけれど，
もとデータで欠測していたり，
割算値にしているといったことが原因．
-
状態遷移の統計モデルあてはめ．
上の修正によって，
まあ，
前よりは妥当な結果になった．
侵入定着確率は，
近傍アリ密度と幹線道路密度によって高められる，
と．
なぜアルゼンチンアリが太い道路が好きなのか，
まだわからないけれど．
-
1700 研究室発．
AIRDO 航空券の領収書を受けとるべく，
ひたすら南下．
昨年，
東日本大地震の翌日に開通した (開通式典はとりやめ)
全長 500 m 超の地下歩道．
JR 札幌駅と大通り公園を南北に直結する．
信号まちをくらわないし，
積雪の影響を受けないので，
いつでもすばやく移動できる．
外はマイナス 6 度ぐらいだけど，
気温は 10 度以上．
けっこう皆さんあるいてるなあ．
こういうところにありました．
とりあえず 無人のカウンターがあった．
声をかけると奥からヒトがでてきて，
手書きの領収書を発行してくれた．
経費節約的といいますか．
-
N 1 W 2 を出発して，
こんどは西北方向にススみ，
いつものごとくコープ札幌植物園店 (N 5 W 11) で買いもの．
1755 帰宅 (N 4 W 13)．
晩飯の準備．
晩飯．
-
明日の飯の準備．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.5 合．
ネギ・シジミの味噌汁．
モズク酢・カニカマボコ．
- 昼 (1230):
研究室．
米 0.7 合．
シシトウの炒めもの．
セロリ・ニンジン・豚肉のソース炒め．
ゴボ天．
プチトマト．
ミズナ．
リンゴ．
- 晩 (1900):
ハクサイ・ネギ・ブナシメジ・カキ・鶏だんごの鍋．
卵雑炊．
2012 年 02 月 09 日 (木)
-
0620 起床．
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
0820 研究室着．
-
いま検討すべき問題なのかどうかはともかくとして
……
アルゼンチンアリ (アリ) モデルの改良が気になっている．
調査は廿日市市の公園で実施されたわけだが，
公園だけがアリにとっては必ずしも特殊な生息地ではない
はず (いや，もしかしたら特殊なのかも)．
ただ単に誰の許可も得なくてよいので，
そこで調査した
……
ぐらいに考えてみよう．
-
つまり公園でつかまえたアリの個体数は
「そのへん一帯のアリ密度の一部を代表するもの」
ぐらいの意味あいと考えればよい．
ということで，
このような局所的なアリ密度なるものが，
まあ動態モデルの基盤となるわけで．
このあたり，
やや空想的ではあるけれど．
-
しかし，
シミュレイションによる予測などといったものを提示したいのであれば，
そういうふうに考えるしかない．
今回のデータのあつかいかたでいうと，
500 m 四方の「グリッド」のアリ密度を考えなさいということで．
-
問題は，
このグリッドのアリ密度と，
観測地点であるグリッド内の公園のアリ密度の対応関係だ．
いまやってるモデルとしては，
グリッドアリ密度が公園のアリ密度に影響するといったハナシになっている．
-
しかし，
人間 (造園業者?) によって，
アリ無し地域の公園にアリがもちこまれた場合
……
いや，
人間によるアリ輸送はとくに公園だけに限定されるわけではないな．
-
パラメーター推定だけやって終わるなら，
このあたり何とでもごまかせるんだろうけれど
……
予測までやるとなると，
きちんとグリッドとグリッド内公園の関係，
グリッド間の関係などを考えないといけない．
-
グリッドというのは何だかばかげたデータ構造ではあるが
……
ともかくアリ数データが限定的であり，
しかもどのような空間スケイルでみても局在しているとしかいいようがないので，
グリッドのような粗視化が必要とされるわけで．
これはこれで便利なのだけど，
このような空間における autonomous な動態モデルの設計がよくわからない．
そして，
グリッドのアリ密度はどうやっても観測できないので，
公園アリ密度とうまく対応づけないといけない．
-
などとうだうだ考えているうちに午前おわってしまった．
昼飯．
-
午後もちょっとアリ問題を考える．
1990 年代初めのアルゼンチンアリ侵入時期にさかのぼってモデリングする，
という方法もありえるかなあ
……
現在のモデルでは，
各公園の状態遷移率が環境あれこれに影響されるわけだが，
そうではなくグリッドの状態遷移をそのようにして，
各公園の状態変化はグリッド密度だけに影響されるといった方式も可能だな．
少なくともこう変えれば，
多少は首尾一貫するんで，
説明の面倒は軽減される．
-
ちょっと気分を変えるために
……
3 年間弱にわたってココロの重荷になってた
石狩浜アリ論文の再投稿問題について検討するか．
これは Behavioral Ecology and Sociobiology
で reject されて，
さっさと修正して別 journal に投稿するつもりだったんだけど，
岩波本すとーむにまきこまれて
……
いや，
これは 2009 年後半からで，
前半はその他あれこれ (耐震工事疎開とか)
にヤラれていたな．
ともかく，
不本意ながら放置されていたもの．
-
この reject 時のぎょーむ日誌をみると，
何やらモデル改良を考えているみたいだな
……
まあ，
よけいなことはヤメといたほうがよさそうだ，
というのが 2012 年の心境．
-
ともかく 2009 年・2010 年・2011 年はいろいろたいへんな
3 年間だったということなのでしょう
……
-
さて，
次の投稿先としては
Animal Behaviour
はどう?
とかいうハナシがあるんだけど，
うーむだいじょうぶかなぁ
……
-
ちょっと検索してみると，
この 3 年間に
nestmate な識別の論文がいろいろ出てるぢゃん
……
とゆーうつな気分になったのだが，
意外と interspecific なハナシが多い．
石狩浜のそれは，
いちおう intraspecific だからね．
アリからみれば，
どっちもケンカのハナシではあるが，
人間はそこに勝手な文脈をつけて理解したふりをするので，
あたかもかなり異なる現象であるかのように作文される
……
-
まあ，
まじめに最新文献の探求とかやめとこう．
かなりしんどいし．
しかし，
文献よみが足りなくてまた rejct をくらう可能性，
かなり高いような気がする．
まあ，
計算傭兵の弱点ということで．
-
ともかく，
現状の原稿をよくよみ，
よけーな discussion を削るのがよかろう
……
と前回 reject 時の私も言ってくれてるわけだし．
-
1750 研究室発．
買いもの．
1820 帰宅．
晩飯の準備．
晩飯．
-
またいろいろ飯の準備．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.5 合．
ネギ・豆腐の味噌汁．
モズク酢・カニカマボコ．
- 昼 (1210):
研究室．
米 0.7 合．
レンコン・ニンジン・豚肉の炒めもの．
ゴボ天．
ジャガイモ．
ミズナ．
リンゴ．
- 晩 (1910):
米 0.5 合．
ネギ・豆腐の味噌汁．
キャベツ・タマネギ・ピーマン・豚肉の蒸しもの，
ネギソース．
ナガイモすりおろし．
2012 年 02 月 10 日 (金)
-
0620 起床．
かなり寒い．
ラジオ体操．
弁当の準備．
朝飯．
コーヒー．
0750 自宅発．
曇．
0820 研究室着．
-
廿日市アリモデル改良で，
ベータ分布
についてちょっと考える．
α, β というパラメーターを与えるのがふつーのやりかたなのだが
……
ちょっと別の与えかたもありえる．
たとえば，
平均を m とすると分散が m (1 - m) に比例すると考えることもできる．
このときの比例定数は 0 より大きく 1 より小さい．
つまり，
ベータ分布の分散の上限は m (1 - m)
……
これはちょっと以外というか，
α と β がゼロに近づくときには，
もっと大きくなると思っていたんだけど．
m が 0 とか 1 にちかいときでも分散の上限は 1 か?
α > 1 かつ β > 1 の場合
(密度関数が「ひと山」型になりうる場合)
の分散の上限は m (1 - m) / 3 となるわけか．
まあ support がせまいせいなのかな．
-
石狩浜アリ原稿みなおし
……
3 年間ちかくも放置してたので，
やっぱりわけわからん．
-
窓の外は雪．
今日はぜんぜん気温が上がらないとゆー予報．
雪まつりにきてるヒトたち，
かなりたいへんかも．
-
論文原稿みなおしがススまず，
ちょっと現実逃避にはまってしまった
……
自分用の top page に javascript 時計を追加．
-
しかし，
考えてみれば，
上のごとく見苦しく限定的な時計を表示せずとも，
Chrome Extensions で
overlay clock
と検索したら，
いろいろでてくるぢゃん
……
一見あやしげなれど
じつはあれこれカスタマイズできる
``Overlay Hanzi | Kanji Standard Clock''
をインストールしてみる．
これでいいぢゃん
……
Gmail 上では表示が抑止されるけど．
-
昼飯．
雪はやんでる．
ちょっと散歩．
-
寒い中をふらふらとさまよっていると，
廿日市アリモデル改良案で検討していた状態遷移の対応関係が整理できた．
|
公園アリ密度状態
|
グリッドアリ状態
|
| 高密度 |
侵入・定着ずみ |
| 中密度 |
侵入・定着ずみ |
| 低密度 |
まだ侵入してない |
公園アリ密度状態はグリッドには影響しない
--- 公園とはグリッド内の「一点」にすぎないので．
「まだ侵入してない」グリッドでもアリは公園内で「ごくまれに」見つかることにする．
いっぽうで，
グリッドアリ状態が「侵入・定着ずみ」であることは，
各公園における「低密度→中密度」遷移の必要条件とすればよい．
-
以前のモデルでは公園の状態変化の速度が周辺の環境に依存していたわけだが，
改良版ではグリッドの状態変化だけが依存すると仮定すればいいんじゃないかな．
これで，
モデルの「わけのわからなさ」が減ってくれたらいいんだけど
……
-
ひとつ重要なのは，
このようにモデルを変更することによって，
autonomous & forward & grid-based なシミュレイションができるようになる
…… はず．
-
さてさて，
また石狩浜アリけんか実験の「脱」アリ作業にもどりますか
……
-
自分で書いたアリ論文の敵国語作文，けっこう意味不明だ…… reject もむべなかるかな．
(2012-02-10 15:03:43)
-
なんだか読んでてキモチわるいんだよね．
(2012-02-10 15:04:11)
-
discussion 後半の 4-5 パラグラフは削除してもよさそう．
2008 年ごろに苦闘しながらアリ文献を解読しつつ，でっちあげた考察ではあるのだけれど．
(2012-02-10 16:02:57)
-
ともかく現在 9000 words ちかい本文を 6000 まで圧縮しなければならないわけで．
(2012-02-10 16:03:44)
-
introduction から統計モデリングのご利益を説いた 3 パラグラフも削除するか…… Animal Behaviour ではあいかわらず Bayes な統行動解析の統計モデルは使われてないみたいだけど．
(2012-02-10 16:08:04)
-
とりあえず，ためしにばさばさと切ったら，ちょうど 6000 ぐらいになった．
コンパイルされないメモ書きとかがあちこちにあったので．
(2012-02-10 16:20:26)
-
パラグラフをけずるということは……修正すべき文章の分量が減るわけで，その点は気楽になったかな．
(2012-02-10 17:08:07)
-
なぜかしら晩飯は北大のちかくで食べることに．
[ごはんやデンスケ]
N 13 W 3 にある定食屋．
毎週金曜日はばら寿司の日．
-
2100 研究室発．
2140 帰宅．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.5 合．
ネギ・油揚げの味噌汁．
ゴボウ・ベイコンのバルサミコ炒め．
- 昼 (1230):
研究室．
米 0.7 合．
ゴボウ・ベイコンのバルサミコ炒め．
カニかまぼこ．
トマト．
ジャガイモ．
レタス．
リンゴ．
- 晩 (1900):
デンスケのばら寿司．
