前回は御病気の統計でしたが、今回はducksです。duckといえば北京ダック食いて～などと思う不埒な年寄です。しかしマガモとアヒル(duck)の関係も知らなかった年寄です。合鴨もいたな～（旨いな。）しかし今回のサンプルデータ・セットは羽毛を採る方かもしれません。オナガガモとマガモの「ハイブリッド」種についてです。
※「データのお砂場」投稿順Indexはこちら
このところ、R言語のパッケージbootに含まれているサンプルデータセットを abc 順にみています。今回のデータセットは ducks です。
マガモ、アヒル、オナガガモ
duck（アヒル）とマガモさらに言えば合鴨との「関係」、まったくもって理解しておりませなんだ。
アヒルはマガモを家禽化したもの。生物学的には同じ種。
なんだそうです。見た目の色は違うが同じ生物種の「品種」レベルの差であると。さらに合鴨は「野生の」マガモとアヒルの交雑品種だということなので込み入ってます。調べるに以下の『大阪市立自然史博物館』様のページから１か所引用させていただいたお言葉がその辺の関係をあらわしておりますで。
ご質問の点に関しては、そこで言われている“カモ”が何を指しているのかがよくわからないので、きちんと答えるのは困難です。カモの養殖場で飼育しているのが、野生のマガモを捕獲したものだとしても、飼育化で何代も繁殖させれば“アヒル”になるのは当然だと思います。
そういう関係だったのね。どうも動物を家畜化すると色が変わったりするものらしいし。。。図鑑的な絵がみたいところね。アヒル(duck)の絵が『畜産Z00鑑』様の以下のページにありました。
一方、今回のサンプルデータセットで扱っているのは、マガモとオナガガモの「ハイブリッド」アヒルです。英名と列挙すれば以下の通り。
-
- mallard マガモ
- pintail ducks オナガガモ
『サントリー』様 日本の鳥百科に絵と解説記事がありました。
ハイブリッド・アヒルの行動と羽毛の特徴
さて、サンプルデータセットの解説ページは以下です。
Behavioral and Plumage Characteristics of Hybrid Ducks
これによるとマガモとオナガガモのハイブリッド・アヒル共についてあるインデックスで評価しているのです。そのインデックスは以下のとおり
-
- plumage、外観的（羽毛）が似ている ０から２０
- behaviour、行動的に似ている ０から１５
いずれも数が小さい方がマガモに近く、数が大きいとオナガガモに近いという値になっているらしいです。詳しいことはまったく書いてないです。
先ずは生データ
library(boot)を忘れずに、サンプルデータセットをロードしたところが以下に。形式は普通のデータ・フレームであります。ごくごく小さいデータセットなので、全データ表示しても大した量じゃないです。
とりあえず、２つのインデックスがどんな感じなのかboxplotしてみました。
boxplot(ducks)
プロットが以下に。
解説によるとplumargeの方は０から２０と書いてあったけれども、どっちも１５までしかないように見えるんだが。ハイブリッドだと典型的なオナガガモの２０というインデックスは現れないのか？
上記のboxplotでは、plumageとbehaviourの間の関係が分からないので２次元プロットしてみました。
plot(ducks, col=2, pch=16)
プロット結果が以下に。
上記をみるとなんとなく「線形な回帰」できそうな雰囲気。そこでlibrary(ggplot2)して、以下のように回帰直線を引いてみましたぞ。
p0 <- ggplot(ducks, aes(x=plumage, y=behaviour)) + geom_point() + geom_smooth(method="lm", formula=y~x)
p1 <- p0 + xlab("plumarge") + ylab("behaviour")
p2 <- p1 + ggtitle("Behavioral and Plumage Characteristics of Hybrid Ducks")
p2
グラフが以下に。
だいたい１直線に載っている感じがしないでもない。。。
ブートストラップ法を適用してみる
このplumargeとbehaviourインデックスについて以下のように無理やりブートストラップ法を適用してみました。いいのかそんなことで。
ducks.work <- ducks
ducks.work$ratio <- ducks$plumage / ducks$behaviour
colnames(ducks.work) <- c("plumage", "behaviour", "ratio")
ratioMean <- function(d, i) mean(ducks.work$ratio[i])
rm1000 <- boot(ducks.work, ratioMean, R=1000, stype = "i")
plot(rm1000)
結果のグラフが以下に。
plumarge対behaviourはだいたい1.2くらいの比率になる感じね。
羽毛を採るのにこういう「統計を駆使し」ているのかしら？知らんけど。でも行動指標でなんだ？