拜估狗大神的幫忙 只要鍵入關鍵字 "大麥蟲繁殖"
就會有一大堆網友的分享文章 這些文章會告訴您一件事
「大麥蟲必須要住單人房才會有化蛹動作!!」
所謂的單人房就是單獨一隻大麥蟲關緊閉啦!!
因為大麥蟲化蛹時會靜止不動 這時會被同伴當成食物吃掉
所以快要化蛹的大麥蟲會像遊魂一樣 一直在飼料表面閒晃 一副想要爬出飼育盒
想躲起來安全化蛹 免得被同伴吃掉
所以群居的大麥蟲 是看不到他們化蛹的狀態的
可是不化蛹怎麼有成蟲來繁衍後代呢?
所以讓大麥蟲化蛹是個關鍵步驟
而本篇主題是要探討有沒有甚麼簡單條件可以讓大麥蟲快速化蛹呢?
我們很簡單的做比較: 大麥蟲分別在有/無麥麩的隔離盒中, 化蛹時間是否不同?
(上圖左是有加麥麩, 上圖右是沒加麥麩)
我們驗證方式就是隨機抓60隻大麥蟲, 分成兩組
30隻的隔離盒沒有放麥麩
30隻的隔離盒有放麥麩(麥麩重量約2-5克, 原則上就是隔離盒底部有一層薄薄麥麩即可)
然後從單獨隔離起算, 紀錄他們的化蛹時間
藉由統計驗證方法驗證兩件事:
*利用雙比例檢定驗證有/無麥麩在隔離盒的化蛹比例是不是不一樣?
*利用存活分析中的無母數方法計算,隨著時間過去,有/無麥麩隔離盒的大麥蟲化蛹比例變化
總結論是: 想要讓大麥蟲快速進入化蛹狀態, 請務必要在隔離盒內加入麥麩!!!
*有加麥麩的隔離盒化蛹率較高
*有加麥麩的隔離盒, 大麥蟲死亡率較低
*有加麥麩的隔離盒, 大麥蟲可較早進入蛹期
有無加入麥麩到底差異多大? 有興趣的讀者請繼續看下面的數值分析
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~正式數據分析起始線~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*利用雙比例檢定驗證有無麥麩在隔離盒的化蛹比例是不是不一樣?
我們最後收集到57比數據, 30筆是無麥麩的大麥蟲數據, 27筆是關於有麥麩的大麥蟲數據
無麥麩的化蛹比例是 18/30=60.00%
有麥麩的化蛹比例是 26/27=96.30%
用Minitab做雙比例檢定的結果是
双比率检验和置信区间
X N 样本 p
無麥麩 18 30 0.600000
有麥麩 26 27 0.962963
差值 = p (無麥麩) - p (有麥麩)
差值估计值: -0.362963
差值的 95% 置信上限: -0.204161
差值 = 0(与 < 0) 的检验: Z = -3.26 P 值 = 0.001
Fisher 精确检验: P 值 = 0.001
由上結果顯示, P值達0.001, 這表示有/無麥麩化蛹的比例差異已經不是隨機產生的!!
再由信賴區間上限得知, 有/無麥麩的化蛹比例差距高達20.41%
這數據可不得了了,
100隻的差距就是20隻化蛹, 20隻如果有一半是母的
一隻母的生一斤的大麥蟲, 那差距就是20斤大麥蟲 不可謂不大啊!!
再來, 我們關心的是
*利用存活分析中的無母數方法計算,隨著時間過去,有/無麥麩隔離盒的大麥蟲化蛹比例變化
講白話就是, 我們要關心的事是:
到底加入麥麩是否可以讓大麥蟲早點進入蛹期呢?
我們在收集數據時, 除了記錄了化蛹時間外
我們也發現, 被隔離的大麥蟲還太嫩, 所以可能會有蟲蛻皮,
或者是沒蛻皮, 可是死亡了!!
發生以上兩件事, 就表示他們化蛹時間是在這兩件事之後才會發生
我們已經無從觀察了, 也就是右設限(Right-Censoed)了
所以, 要回答這個問題, 我們可以使用存活分析(Survival Analysis)來回答
先比較隨大麥蟲在隔離盒的時間過去, 已經完成化蛹的比例是否不同
> survdiff(Surv(Time,Status)~Recode.Food)
Call:
survdiff(formula = Surv(Time, Status) ~ Recode.Food)
N Observed Expected (O-E)^2/E (O-E)^2/V
Recode.Food=無麥麩 30 18 31.8 5.99 24.6
Recode.Food=有麥麩 27 26 12.2 15.60 24.6
Chisq= 24.6 on 1 degrees of freedom, p= 6.93e-07
這一個方法叫做對數秩檢定(Log-Rank test)
是用來比較完成化蛹比例在不同時間是否不同的檢定方法
P值, p=6.93e-07, 告訴我們, 有無麥麩對話蛹的影響非常顯著!!
一定有一個方法可以讓大麥蟲早日進入蛹期!!
所以我們可以再進一步劃出完成化蛹比例圖形
請看....
紅色的線是有加麥麩的化蛹表現
黑色是沒加麥麩的化蛹表現
結果已經很明顯:
有加麥麩的隔離盒, 7-8天就有蛹出現了
反觀沒加麥麩的, 超過10天才有蛹出現
各只看紅色與黑色的第二條線
完全沒有交錯, 這表示: 有加麥麩的化蛹表現一直遙遙領先....
最後, 利用R程式得到50%的大麥蟲化蛹時間的比較
> KM.estimate
Call: survfit(formula = Surv(Time, Status) ~ Recode.Food, conf.type = "log-log")
records n.max n.start events median 0.95LCL 0.95UCL
Recode.Food=無麥麩 30 30 30 18 15.21 13.04 18.4
Recode.Food=有麥麩 27 27 27 26 9.71 9.04 11.1
看到沒?
無麥麩, 50%大麥蟲化蛹時間最長是18.4天 最短13.04天
有麥麩, 50%大麥蟲化蛹時間最長是11.1 最短是9.04天
有加麥麩真的是遠遠大勝啊!
有加麥麩, 只花了11.1天, 就有一半的大麥蟲化蛹了........
沒加麥麩, 最長需要18.4天才有一半化蛹.....
所以結論就是, 加麥麩吧....
後記:
*要挑選在飼料表面遊走的大麥蟲, 化蛹率才高
*加麥麩可以讓被誤抓的大麥蟲有東西吃, 免得餓死
*我們後續有實驗過, 最終齡的大麥蟲在有麥麩的隔離盒, 第三天就會蜷曲, 進入化蛹狀態
沒有麥麩的, 三天過去, 所有大麥蟲還是不安分地動來動去
所有有賣膚真的可以讓他們放心地化蛹
上圖是進入化蛹狀態的大麥蟲, 他們要化蛹前都是捲起來的
上圖是已經排出黑便的大麥蟲, 排黑便是化蛹前必經的
(注意看大麥蟲屁股後面的小小黑塊...)
化蛹後的狀態......
剛羽化的大麥蟲成蟲
羽化後就集中放在飼育盆, 開始交配生出一大堆小大麥蟲了....
留言列表