dynamicsoar's log

主に研究関係のメモ

読んだ：A general scaling law reveals why the largest animals are not the fastest

journal club biomechanics swimming flight

論文の情報

前書き

僕は動物の飛行の力学（主に流体力学）が専門で、最近は遊泳にも手を出し始めた程度で、歩行・走行は専門ではない。

統計については、最近ようやく勉強を始めたところだけれども、全くの素人であり、この論文の統計手法に関する評価をすることはできない*1。

それから、筋肉の話がでてくるけど、生理学についても素人なので、以下では正確性を欠く記述をしている可能性がある。

内容まとめ

概要

いろいろな動物について、水平移動の最大ダッシュ速度（≠巡航速度）を縦軸、体重（質量）を横軸にとった散布図を描く。そうすると、だいたい、小さいほど遅くて、大きいほど速い傾向がある*2。この傾向を、うまく表すような、シンプルな曲線（を表す数式）は、従来の理論では、power law（べき乗則？）といって、「大きければ大きいほど速い」ような数式が提案されていたようだ (第1式; Bejan & Marden, 2006, J. Exp. Biol., 未読)。本研究では、よりよいフィッティングとして、「最大の動物よりもちょっと小さいあたりの動物が一番速い」ような、形としては ∩ のようになる曲線（を表す数式）を提案している。実際に、複数のモデルを比較したところ、提案している新しいモデルが一番よいフィットになった。また、歩行・遊泳・飛行という3つの運動（移動）方法の全てについてこのモデルはフィットした。さらに、いまの世の中にいる生物だけでなくて、化石になった生物、たとえば恐竜の速度をも予測できるはずだとしている。実例として、いくつかの恐竜について、形態学をもとにした推測値と、モデルによる推測値がわりとよく合っていることを確認している。

もうちょっと踏み込んだ内容

当然だけど以下は僕の理解。生理学的な部分などで間違いがあるかもしれない。

速度に「小さいほど遅くて、大きいほど速い傾向」があるのはなぜだろうか。この論文では以下のように説明してる。運動には無酸素運動 (anaerobic) と有酸素運動 (aerobic) の2種類があり、それぞれ筋肉の速筋と遅筋に対応している。速筋は素早い運動を可能にするけれど継続時間が短い；遅筋は逆に遅いけれど継続時間が長い。いま注目しているのはダッシュ速度なので、無酸素運動・速筋が関係していそうだ。ポイントは、速筋は継続時間が短いということ。この影響を考慮したのが今回提案するモデル。

要点は Figure 1 にまとめられている*3。これは上段の左右 (a, b) と下段の左右 (c, d) がペアになっている。まず、上段は（従来の）理論的な予測で「無限に加速できるとしたら、速度は体重の順番になる」ということ。ここでもう一つのポイントは、最大速度に到達するのは体重が大きいほど遅いという点。これは一種の2乗3乗則ということでいいのだと思う。力は筋肉の断面積つまり長さの2乗に比例するが、質量（体重）は3乗に比例する。したがって、加速度 = 力/質量 は長さの1乗に反比例する…つまり、大きいほど加速が遅い。まぁ、直感的にも、小型飛行機と大型ジェット機の離陸とか見てればわかるんじゃないかと思う。

ところが実際には、速筋が蓄えられるエネルギ (ATP) には限界があるので、有限の、けっこう短い時間しか加速することはできない。この影響を考慮すると、体重がとても大きな（重い）動物は、加速度が小さいために、速筋のエネルギを使い切った時点ではまだ十分に加速しきれていない、ということになるようだ（下段）。逆に、これはハッキリとは書いていないのだけれど、あまりにも小さい（軽い）動物だと、速筋のエネルギを使い切るよりも前に加速しきってしまうようだ。こうしたことから、結果として「実現できる最大ダッシュ速度」は、体重の順番にはならず、「ちょうど速筋のエネルギを使い切る頃に十分に加速しきる動物」すなわち「一番重いやつよりもちょっと軽いやつ」が最速、ということになる。

このように、提案している新モデルでは筋肉のことを考慮しているけれど、地面から受ける反力とか水流や風の影響は全く考慮していない。個人的にはこのあたりが気になるのだけれども。

手法

  • 文献を調査して、622データ（474種）の最大（ダッシュ）速度と質量の情報を取得。
  • 提案モデルを含む、複数のモデル（数式）をこの速度-質量の関係にフィッティング（Rを使用）。
  • フィット具合は Bayesian information criterion (BIC) で評価。
  • Residual variation analysis → ここは僕の知識不足でよくわからず。内温性・外温性の比較もやっている (figure 3) ようだがいまいちわからなかった。

感想

たまには時系列で感想を綴ってみるのもいいかもしれないと思うのでやってみる。

twitter で存在を知り、AFPBB のニュース記事を読み、論文のアブストと図をチラッと見た時点では「おっすごそうじゃん」という印象だった。「歩行・遊泳・飛行という複数の移動様式における最大速度」がひとつのモデル（数式群）で表現できるという売り込み。どうやら、筋肉の収縮かなにか、生理学的な発想がもとになっているようだ。「あぁ、そっち系なのかな。流体力学はどれくらい考慮されているのだろう？」僕は力学屋なので、このモデルというのが力学をどれほど考慮しているのか、あるいはデータに合うようにフィッティングしただけで力学的メカニズムは考慮していないのか、という点がまず気になる。歩行では地面との接触による地面反力 (ground reaction force, GRF) だけを考慮して流体力学的な考慮はしなくていいのかもしれないが、遊泳では水の抵抗（抗力）を受けるし、飛行では抵抗に加えて重力に対抗するだけの揚力を生み出す必要がある。そしてこれら流体力学的な力というものは、物体のサイズによって影響を受ける。そこでまず、論文本文をダウンロードして、いくつかのキーワード（やその断片）で全文を検索した。具体的には、drag（抗力）, lift（揚力）, force（力）, fluid（流体）, aero*4, hydro*5 あたりで検索した。すると force は多少ヒットするが3/5が muscle force であり、aero はヒットした9つ全てが aerobic か anaerobic, hydro は参照文献のタイトルで1ヒットのみ。よって「あぁ、やはり筋肉の生理学は考慮しているようだが、流体力学は考慮していないようだな」という事がわかった。この時点で熱狂は少し落ち着き、「うーん、深く読むのは後でいいかな」となった。僕が興味あるのは、流体力学的な考慮がなされたモデルだからだ。

ところがアブストをよく読むと、このモデルは「474種 (species)」ものデータで支持されているという。これがひっかかった。少なくとも飛行に関して、正確な最大速度というのはそんなに多くは計測されていないという印象が僕にはあった。「このうち飛行ってどれくらいあるんだろう？」という疑問がわいた。そこで、イントロダクションも結果もすっ飛ばして手法 (Methods) を見に行った。ここで最大速度は具体的には文献から “maximum speed”（最大速度）、“escape speed”（逃避速度）、“sprint speed”（スプリント速度…おそらく短距離ダッシュのことだろう）のキーワードで検索したとある。さらに以下の2つの速度を除外している：

  1. 垂直方向の速度。主に鳥の文献に多い。要は急降下で最速になるということだろう。重力加速度により速度を得てしまうが、これは今回のモデルに入っていないため除外するとのこと。
  2. 通常速度のうちの最大値。これは「including also dispersal and migration（分散や渡りを含む）」とあるので、要は、長時間に渡って維持が可能な最大巡航速度のことだろう。

本文をすっ飛ばしていたので、ここで初めて最大速度というものが実際には「水平方向の」「ダッシュ速度」だとわかった。前者はなんとなくそうかなと思っていたが、後者が意外で、僕はここを読むまでは maximum speed とは「最大の巡航速度」のことだと思っていた*6。

続けて、"622 data points for 474 species (see Supplementary Table 1 for an overview)" とある。どうやら本文中にはデータ数の内訳がなく、supplementary material を見ろということのようだ。

…なんか時系列で書くのに疲れてきたので、ここからはまた要点のみの記述にします…

出典について

本文中で使用されたデータセットの出典は supplementary information*7 に明記されている。しかし、

  • データ数の少なさ
  • 出典調査のずさんさ

の2つの問題があると思う。

追記：出典の表をPDFでのみ提供していることも大きな問題だわ…。CSVにすべきだ。

遊泳と飛行のデータ数

走行（歩行）に比べて遊泳と、とくに飛行のデータ数が少ない。具体的には、

  • 走行：458データ（節足動物50・鳥類3・哺乳類261・爬虫類144）
  • 遊泳：109データ（節足動物1・鳥類5・魚類81・哺乳類16・軟体動物5・爬虫類1）
  • 飛行：55データ（節足動物19・鳥類29・哺乳類9）

となっている (Supplementary Table 1)。統計学的にどうなのかはわからないけど、まぁ、当然、多いほうがよりよいフィッティングになるんだろうってことくらいは想像できる。

個別に見ると、たとえば遊泳する節足動物と爬虫類のデータが1ずつしかなくて、これらはモデル生成に使っていない、といったことも disputable かもしれない（ないものはないのでしょうがないが）。

出典調査のずさんさ

孫引きやひ孫引き（？）が多く、出典調査がずさんである（手間を惜しみすぎている）。したがって出典の信頼性に対する配慮が充分なされているかが疑わしい。

オンライン・ウェブページを出典とすること

出典 96 から 103 はオンラインのウェブサイトとウェブページである。基本事項として、ウェブページの情報は永続性に乏しく、ほとんどの場合は自らの研究ではなく他人の情報の転記や引用であり、査読を経てもいない。当然、扱いには注意が必要だ。

たとえば出典101は www.speedofanimals.com というサイトで、おおもとの出典はなんだろうと思うと、一番下の方に "Animal descriptions are from Wikipedia" と書いてある。英語版ウィキペディアは比較的出典がきちんと整備されている印象があるので、値としては間違っていないのかもしれない、と思いたくはなる。ところが、speedofanimals.com のリンク先URLは www.wikipedia.org とだけなっており、個別の動物についてのリンクではない。英語版であるという保証もない。いつ時点という情報もない。Wikipedia は誰でも書き換えられる wiki なので、「いつの時点の版（バージョン）」という「版指定」がないと、出典としての正確性が劣ることになる*8。

しかしそもそも、研究者として、ウィキペディアから転記したサイトをさらに孫引きするというのは正直にいって驚愕した。これはつまり、ウィキペディアにある情報を自分たちでチェックすることを怠り、このサイトに丸投げしているということだ。正しい態度は、ウィキペディアへは当然あたり、さらにそのウィキペディア記事が論拠としている原著論文へあたることだろう。こんなのは基本中の基本で、研究室なら学部生がゼミで指摘されるくらいのことではないだろうか。それを、データをフィッティングしている、データに強く依存した論文で、そのデータをきちんと遡っていない、なんていうのは、本当に、ちょっと、どういうことだ…。しかも、Nature 関連ジャーナルに載るような論文で*9。

「そうはいっても、そういう疑わしいデータは、450とかあるデータ点のうちの、ごくごくわずかなんでしょ？」と思うだろうから、速度だけについて、いくつあるかを数えてみた*10：

  • 走行：27/458 (5.9%). 内訳は、哺乳類23/261*11, 鳥類2/3, 爬虫類2/144.
  • 遊泳：17/109 (15.6%). 内訳は、鳥類3/5, 魚類9/81, 哺乳類4/16, 軟体動物1/5.
  • 飛行：10/55 (18.2%). 内訳は、全て鳥類10/29.

となっていた。特に飛行では20%近い。だから直ちにダメだというわけではないが、もともと少ないデータのうちの2割をきちんと原典まで遡って確認していない、というのは、信頼性を損ねる行為だと思う。

書籍を出典とすること

これ自体は批判されるような行為ではないのだけれど、やはり査読を経ていない（だろう）し、孫引きなので、もっと sharper な引用の仕方をして欲しいところだ。つまり、多くの本というのは自分で独自に計測したデータを載せることはなくて、さらなる出典があるだろうから、そちらの方を引用すべきではないか、ということ。たまにそうやって調べていると伝言ゲームの過程でデータの転記ミスに気づくこともあるし。

なお、飛行では10/55（節足動物2/19, 鳥類8/29）が出典11番の本に依存している (http://press.princeton.edu/titles/162.html)。出典11番のところには2015年発行とあるが、実際には初版1965年の本の復刻であり、データとしてはかなり古い。他の分野は知らないが、動物飛行のデータで1965年というのは、計測方法などの面から不安が残る。さらにこの本のタイトル "An Essay on" というのがまた不安を掻き立てるが、中身は確認していない。鳥類については、ウィキペディアとこの本とで18/29=62%のデータを占めることになるのだが、僕なら怖くて原典を徹底的に確認しないとフィッティングなんてする気になれない。

あとは書籍じゃなくて査読済み論文だけど孫引き、ってのもあって、それはさすがに（原典の年代や手法などを一度確認しておけば）いいかなぁ、と思います。走行の出典44番と80番がかなり多い模様。

twitter でのツッコミ

以下のようなツッコミが見られた（他にあったら追加するかも）。

ハヤブサ

うーむ。まぁニュース記事はしょせんニュースなので、「論文をよく読んでなかったんですねー」で終わるかもしれないけれど、ちょっとそれだけじゃないんだよなぁ…。このハヤブサの出典は AllAbout からの孫引き： www.allaboutbirds.org/guide/Peregrine_Falcon/lifehistory ←このページには出典が5つ挙げられていて、このどれかが原典だろうとは思うが…（そしてその本の中にまた本当に原典が…）。いやいや、待って待って、その前に、AllAbout には 110 km/h という数字はないぞ。引用すると、

reaching speeds up to 112 km/h (69 mph) in direct pursuit of prey

とある。これ自体がまずちょっと問題で、112 km/h = 約69.59 mph (miles per hour) または 69 mph = 約111.05 km/h だから、記述自体が不正確だ。ただ、いずれにしても 110 km/h にはならない。これが 110 km/h になった可能性は2つありそうで、ひとつは論文執筆者がアクセスした時点では 110 km/h と書かれていた可能性*12。もうひとつは、論文執筆者が勝手に1の位を丸めたか、69 * 1.6 = 110.4 -> 110 km/h というような雑な計算をしたか。なんにせよ、「雑だなぁ」という印象を受ける。

カジキ

これは気づいてなくて（遊泳ちゃんと追いかけてなかった）、「え！？桁違いじゃん」と思ったけど、これもまた孫引き問題でした。しかもけっこう根が深い。

  • 出典10番のその "Direct ... blue marlin" (Block et al., 1992, J. Exp. Biol.) は black marlin (Istiompax indica) = 130 km/h および blue marlin (Makaira nigricans) = 75 km/h の速度の出典として載っている。この論文はタイトルの通り blue marlin (Makaira nigricans) の計測をしている。そして指摘の通りこの論文における直接計測の最大は 8.1 km/h なのに、どうなっているのか。
  • 本文を「130」で検索すると Table 1 にある Walters (1962)。また、blue marlin の75 km/h は Walford (1937) とある。いずれも Estimate とある。実測値ではないようだ。しかし References を見ても Walters 1962 なんてものはない（たぶん載せ忘れ）。ググるしかない。こっからまた遡ります。
    • なお、この論文の本文で一番関連ありそうなところを抜き出すと、"The fact that they can strike a bait trolled at 800 cm s-1 leaves no doubt that they are capable of short bursts of high-speed swimming, but we saw no such rapid movements in 160 h of speed observations. High-speed swimming and agility must have a place in the repertoire of activities displayed by blue marlin but, during the summer near Hawaii, such behavior is rare. There were only infrequent 10 to 30 s bursts of speeds exceeding 200 cm s-1. Short periods of fast swimming are undoubtedly important for catching certain prey, and hence important in terms of survival, but the characteristic mode of swimming we observed for this species was quite slow." というわけで、800 cm s-1 (28.8 km/h) という値*13はあるものの、計測できたのは 200 cm s-1 程度だったということになっていて、どうも Walters らの estimation は怪しいのではという雰囲気…（ちゃんと読んでないので違うかも）。
  • おそらく Walters, 1962, Am. Zool. http://www.jstor.org/stable/3881203 がそれ。しかしこれもおおもとではない。 Nursall, 1962 が cite されている。しかも、"Nursall (this issue) questions the validity of the reported swimming speeds of 80 and 90 km/hr for tuna, 130 km/hr for marlin and swordfish and of 120 km/hr for sailfish" とある。Block et al., 1992 はこの記述を見ていたはずなのに、Nursall には言及すらしていない…。
  • Nursall, 1962 はこれ https://doi.org/10.1093/icb/2.2.127 *14 であり、ここには、barracuda（カマス）の計測値が異常ではないかという記述に続いて、"Also suspect then are unconfirmed figures for salmon (45 Km/hr = 28 mph), tuna (90 Km/hr = 56 mph) and martin and swordfish (130 Km/ hr = 81 mph) reported by Barsukov (1960)" とある。まだ遡るのか…。
  • ところが References を見ると、Barsukov, V. V. 1960. The speed of movement of fishes. [In Russian]; Priroda 3:103-104 とある。はい出ました、Russian. これはちょっと遡れないですね…。どうやらこれが原典っぽいけども。
  • というわけで、ようやく、この 130 km/h は "unconfirmed figures" である（っぽい）ことがわかりました。

えーっと、何の話だったっけ？あ、つまり 130 km/h が怪しい、と。たぶん 75 km/h も怪しい。そして今回の論文に戻って figure 2d を見ると、シルエットからしてもカジキっぽいやつがちょうどピーク付近にいる。フィッティングカーブの形に効いてるんじゃないの？という気がしてくる。遊泳も…大丈夫なのか…？

その他の感想

全体を通して

どうにも、全体を通しての僕の感想としては、「走行（歩行）だけにしとけばよかった（っぽい）のになぁ」ということになりそう。いや、走行もチェックしてないけども。そうとうに穿った見方をすると、「走行のデータから新モデルでーきたっと。でも『少なくとも哺乳類の少なくとも走行について成立する新モデルです』と言ってもなんかしょぼいし、ちょこちょこっと wikipedia とかオンラインサイトとかで飛行と遊泳のデータもサクッと集めてみるかー」→「おぉ、なんかよさげじゃん。すげぇ発見だよ！Nature（系列）だ！！」みたいな感じじゃないのかなぁ…って気もしてくる。うがちすぎかな。

ただ、流体力学的なものを入れたとしても、定量的に変わるだけで、定性的な曲線の形状自体はこんな感じなのかもしれない、とは思う。つまりフィッティングに使うパラメタが大きく変わるけど、一番の大枠は間違っていないのかもしれない。その「大枠」の提唱こそがすごく評価されることなのだ、そんな細かい出典のことやらなんやらにツッコんでいるからお前は万年ポスドクなのだ、とか言われたら死ぬ。

mechanistic ってなに？

これ連呼されてるんだけど、意味がよくわからなかった。Methods の Model fitting のところを読むと、どうやら「単なる数値へのフィッティングではなくて、生理学というメカニズムを考慮している」ということを言いたいようだ、たぶん（自信なし）。確かに、流体力学は考慮していないとしても、生理学は考慮しているのだから、一歩前進ということでいいのだろう。

恐竜

Table 1 と figure 4 がそれ。Figure 4 でグレーの円がモデル生成（フィッティング）に用いた現生動物のデータで、緑の塗りつぶし三角が別の人が計算した恐竜の速度と質量の関係（Table 5 でいう Morphological models の列の値）。出典は supplementary table 5 とその references にある3つの文献。これら恐竜の値は、今回新たに提案したモデル生成には使っていない、つまり独立に計算された値。その割にはずいぶんよく合ってるでしょ、ドヤ、ということ。一方、従来手法の power law を使うと、Table 1 の Power law の列にあるように、体重が小さいうちは Morphological models と合うけれど、体重の大きな種では全然合わない（過大に評価する）。なぜ figure 4 に power law の線を入れなかったのかはわからない。入れたほうがわかりやすいのに。両対数なので、要は右肩上がりの直線になるってことだろう。

*1:個々のデータ点の数値が公開されているので、統計の自習素材として使ってみるのはいいかもしれないなぁ…というのはちょっと思った。ただ、表がcsvとかではなくてPDFなのだよな…嫌がらせかよ… → あとで出典チェックする時にウンザリした…。

*2:速度を体長で正規化、とかは今は考えない。

*3:たぶんこの論文は open access (OA) じゃないと思うので、読めないって人も多いかもしれないけど、僕は非OA論文の図をネットにバシバシUPすることにはまだ抵抗があるので、図は載せません（おそらく多くの場合 fair use とみなされるのだろうという気はするが）。bioRxivにプレプリントがあるのでそちらをどうぞ→ https://doi.org/10.1101/095018 . パッと見で図はだいたい同じように見える。

*4:aerodynamics（空気力学）, aerodynamic force などにヒットする。

*5:hydrodynamics（水力学）, hydrodynamic force などにヒットする。

*6:僕だったら maximum dash speed とかなにか、定義を読まなくても maximum cruising speed とは明確に区別される用語を使ったかもしれない。

*7:紙の雑誌ではなくオンライン版のみに提供される付録という意味での Electronic Supplementary Material の頭文字から ESM とか呼ばれることもある。

*8:ある人が2071年に、2017年出版の論文の出典を確認しようとして、「ウィキペディアから」とだけ書かれていたらどう思うだろうか？「え、いつの情報だよこれ…」となるだろう。

*9:「査読者なにやってんの」というのはあるかもしれないが、基本的には supplementary material までは査読義務はないはずなので、責めるのは難しい（よくない）かもしれない。

*10:質量は速度に比べると計測が容易で、疑いの余地が相対的に小さいので、とりあえず置いておく。

*11:ただしこのサイトのデータのみに依存した種は5, 6程度で、それ以外は他の出典もある。

*12:AllAbout は wikiではないが、当然、記述は変わりうる。

*13:しかしこれも personal communication なのだが…。

*14:American Zoologist というジャーナルは2002年に名前が Integrative and Comparative Biology に変わっている。毎年1月に国際会議をアメリカのどこかで主催してるのがこの発行団体である Society for Integrative and Comparative Biology (SICB).