㈠ 全基因組選擇之模型篇
在介紹GS模型之前，我們有必要先來了解一下混合線性模型（Mixed Linear Model，MLM）。混合線性模型是一種方差分量模型，既然是線性模型，意味著各量之間的關系是線性的，可以應用疊加原理，即幾個不同的輸入量同時作用於系統的響應，等於幾個輸入量單獨作用的響應之和（公式1）。
𝒚= Xβ + e =𝒃 𝟎 +𝒃 1 𝒙 1 +𝒃 2 𝒙 2 +⋯+𝒃 𝒌 𝒙 𝒌 +𝒆 （公式1）
式中𝒚表示響應變數的測量值向量， X 為固定效應自變數的設計矩陣， β 是與 X 對應的固定效應參數向量；𝒃 𝟎 、𝒃 1 、⋯、𝒃 𝒌 是未知參數；𝒙 𝟎 、𝒙 1 、⋯、𝒙 𝒌 是影響各因素的觀察值；𝒆是殘差。同時需要滿足條件： E(y)=Xβ，Var(y)=σ 2 I， y 服從正態分布。
既然是混合效應模型，則既含有固定效應，又含有隨機效應。所謂固定效應是指所有可能出現的等級或水平是已知且能觀察的，如性別、年齡、品種等。所謂隨機效應是指隨機從總體中抽取樣本時可能出現的水平，是不確定的，如個體加性效應、母體效應等（公式2）。
y = Xβ + Zμ +𝒆 （公式2）
式中 y 為觀測值向量； β 為固定效應向量； μ 為隨機效應向量，服從均值向量為0、方差協方差矩陣為G的正態分布 μ ~ N(0,G) ； X 為固定效應的關聯矩陣； Z 為隨機效應的關聯矩陣；𝒆為隨機誤差向量，其元素不必為獨立同分布，即 𝒆 ~ N(0,R) 。同時假定 Cov(G,R)=0 ，即G與R間無相關關系， y 的方差協方差矩陣變為 Var(y)=ZGZ+R 。若 Zμ 不存在，則為固定效應模型。若 Xβ 不存在，則為隨機效應模型。
在傳統的線性模型中，除線性關系外，響應變數還有正態性、獨立性和方差齊性的假定。混合線性模型既保留了傳統線性模型中的正態性假定條件，又對獨立性和方差齊性不作要求，從而擴大了適用范圍，目前已廣泛應用於基因組選擇。
很早以前C.R.Henderson就在理論上提出了最佳線性無偏預測（Best Linear Unbiased Prediction，BLUP）的統計方法，但由於計算技術滯後限制了應用。直到上世紀70年代中期，計算機技術的發展為BLUP在育種中的應用提供了可能。BLUP結合了最小二乘法的優點，在協方差矩陣已知的情況下，BLUP是分析動植物育種目標性狀理想的方法，其名稱含義如下：
在混合線性模型中，BLUP是對隨機效應中隨機因子的預測，BLUE（Best Linear Unbiased Estimation）則是對固定效應中的固定因子的估算。在同一個方程組中既能對固定效應進行估計，又能對隨機遺傳效應進行預測。
BLUP方法最初應用在動物育種上。傳統的動物模型是基於系譜信息構建的親緣關系矩陣（又稱A矩陣）來求解混合模型方程組（Mixed Model Equations，MME）的，因此稱之ABLUP。Henderson提出的MME如下所示：
式中X為固定效應矩陣，Z為隨機效應矩陣，Y為觀測值矩陣。其中R和G：
其中A為親緣關系矩陣，因此可轉化公式為：
進一步可轉化為：
通過求解方程組，計算殘差和加性方差的方差組分，即可得到固定因子效應值 (BLUE)和隨機因子效應值 (BLUP)。
作為傳統BLUP方法，ABLUP完全基於系譜信息來構建親緣關系矩陣，進而求得育種值，此方法在早期動物育種中應用較多，現在已基本不單獨使用。
VanRaden於2008年提出了基於G矩陣的GBLUP（Genomic Best Linear unbiased prediction）方法，G矩陣由所有SNP標記構建，公式如下：
GBLUP通過構建基因組關系矩陣G代替基於系譜信息構建的親緣關系矩陣A，進而直接估算個體育種值。
GBLUP求解過程同傳統BLUP方法，僅僅在G矩陣構建不同。除了VanRaden的基因組關系構建G矩陣外，還有其他G矩陣構建方法，但應用最多的還是VanRaden提出的方法。如Yang等提出的按權重計算G矩陣：
Goddard等提出的基於系譜A矩陣計算G矩陣：
目前GBLUP已經廣泛應用於動植物育種中，並且因為它的高效、穩健等優點，現在仍飽受青睞。GBLUP假設所有標記對G矩陣具有相同的效應，而在實際基因組范圍中只有少量標記具有主效應，大部分標記效應較小，因此GBLUP仍有很大的改進空間。
在動物育種中，由於各種各樣的原因導致大量具有系譜記錄和表型信息的個體沒有基因型，單步法GBLUP（single-step GBLUP，ssGBLUP）就是解決育種群體中無基因型個體和有基因型個體的基因組育種值估計問題。
ssGBLUP將傳統BLUP和GBLUP結合起來，即把基於系譜信息的親緣關系矩陣A和基因組關系矩陣G進行整合，建立新的關系矩陣H，達到同時估計有基因型和無基因型個體的育種值。
H矩陣構建方法：
式中w為加權因子，即多基因遺傳效應所佔比例。
構建H矩陣後，其求解MME過程也是與傳統BLUP一樣：
ssBLUP由於基因分型個體同時含有系譜記錄和表型數據，相對於GBLUP往往具有更高的准確性。該方法已成為當前動物育種中最常用的動物模型之一。在植物育種中，往往缺乏較全面的系譜信息，群體中個體的基因型也容易被測定，因此沒有推廣開來。
如果把GBLUP中構建協變數的個體親緣關系矩陣換成SNP標記構成的關系矩陣，構建模型，然後對個體進行預測，這就是RRBLUP（Ridge Regression Best Linear Unbiased Prediction）的思路。
為什麼不直接用最小二乘法？最小二乘法將標記效應假定為 固定效應 ，分段對所有SNP進行回歸，然後將每段中顯著的SNP效應相加得到個體基因組育種值。該方法只考慮了少數顯著SNP的效應，很容易導致多重共線性和過擬合。
RRBLUP是一種改良的最小二乘法，它能估計出所有SNP的效應值。該方法將標記效應假定為 隨機效應 且服從正態分布，利用線性混合模型估算每個標記的效應值，然後將每個標記效應相加即得到個體估計育種值。
一般而言，基因型數據中標記數目遠大於樣本數（p>>n）。RRBLUP因為是以標記為單位進行計算的，其運行時間相比GBLUP更長，准確性相當。
GBLUP是直接法的代表，它把個體作為隨機效應，參考群體和預測群體遺傳信息構建的親緣關系矩陣作為方差協方差矩陣，通過迭代法估計方差組分，進而求解混合模型獲取待預測個體的估計育種值。RRBLUP是間接法的代表，它首先計算每個標記效應值，再對效應值進行累加，進而求得育種值。下圖比較了兩類方法的異同：
直接法估計 ，間接法估計標記效應之和 M 。當K=M』M且標記效應g服從獨立正態分布（如上圖所示）時，兩種方法估計的育種值是一樣的，即 = M 。
基於BLUP理論的基因組選擇方法假定所有標記都具有相同的遺傳方差，而實際上在全基因組范圍內只有少數SNP有效應，且與影響性狀的QTL連鎖，大多數SNP是無效應的。當我們將標記效應的方差假定為某種先驗分布時，模型變成了貝葉斯方法。常見的貝葉斯方法也是Meuwissen提出來的（就是提出GS的那個人），主要有BayesA、BayesB、BayesC、Bayesian Lasso等。
BayesA假設每個SNP都有效應且服從正態分布，效應方差服從尺度逆卡方分布。BayesA方法事先假定了兩個與遺傳相關的參數，自由度v和尺度參數S。它將Gibbs抽樣引入到馬爾科夫鏈蒙特卡洛理論（MCMC）中來計算標記效應。
BayesB假設少數SNP有效應，且效應方差服從服從逆卡方分布，大多數SNP無效應（符合全基因組實際情況）。BayesB方法的標記效應方差的先驗分布使用混合分布，難以構建標記效應和方差各自的完全條件後驗分布，因此BayesB使用Gibbs和MH（Metropolis-Hastings）抽樣對標記效應和方差進行聯合抽樣。
BayesB方法在運算過程中引入一個參數π。假定標記效應方差為0的概率為π，服從逆卡方分布的概率為1-π，當π為1時，所有SNP都有效應，即和BayesA等價。當遺傳變異受少數具有較大影響的QTL控制時，BayesB方法准確性較高。
BayesB中的參數π是人為設定的，會對結果帶來主觀影響。BayesC、BayesCπ、BayesDπ等方法對BayesB進行了優化。BayesC方法將π作為未知參數，假定其服從U(0,1)的均勻分布，並假設有效應的SNP的效應方差不同。BayesCπ方法在BayesC的基礎上假設SNP效應方差相同，並用Gibbs抽樣進行求解。BayesDπ方法對未知參數π和尺度參數S進行計算，假設S的先驗分布和後驗分布均服從(1,1)分布，可直接從後驗分布中進行抽樣。
下圖較為形象地說明了不同方法的標記效應方差分布：
Bayesian Lasso（Least absolute shrinkage and selection operator）假設標記效應方差服從指數分布的正態分布，即拉普拉斯（Laplace）分布。其與BayesA的區別在於標記效應服從的分布不同，BayesA假設標記效應服從正態分布。Laplace分布可允許極大值或極小值以更大概率出現。
從以上各類貝葉斯方法可看出，貝葉斯方法的重點和難點在於如何對超參的先驗分布進行合理的假設。
Bayes模型相比於BLUP方法往往具有更多的待估參數，在提高預測准確度的同時帶來了更大的計算量。MCMC需要數萬次的迭代，每一次迭代需要重估所有標記效應值，該過程連續且不可並行，需消耗大量的計算時間，限制了其在時效性需求較強的動植物育種實踐中的應用。
為提高運算速度和准確度，很多學者對Bayes方法中的先驗假設和參數進行優化，提出了fastBayesA、BayesSSVS、fBayesB、emBayesR、EBL、BayesRS、BayesTA等。但目前最常用的Bayes類方法還是上述的幾種。
各種模型的預測准確度較大程度的取決於其模型假設是否適合所預測表型的遺傳構建。一般而言，調參後貝葉斯方法的准確性比BLUP類方法要略高，但運算速度和魯棒性不如BLUP。因此，我們應根據自身需求權衡利弊進行合理選擇。
除了基於BLUP和Bayes理論的參數求解方法外，基因組選擇還有半參數（如RKHS，見下篇）和非參數，如機器學習（Machine Learning, ML）等方法。機器學習是人工智慧的一個分支，其重點是通過將高度靈活的演算法應用於觀察到的個體（ 標記的數據 ）的已知屬性（ 特徵 ）和結果來預測未觀察到的個體（ 未標記的數據 ）的結果。結果可以是連續的，分類的或二元的。在動植物育種中， 標記的數據 對應於具有基因型和表型的訓練群體，而 未標記的數據 對應於測試群體，用於預測的 特徵 是SNP基因型。
相比於傳統統計方法，機器學習方法具有諸多優點：
支持向量機（Support Vector Machine，SVM）是典型的非參數方法，屬於監督學習方法。它既可解決分類問題，又可用於回歸分析。SVM基於結構風險最小化原則，兼顧了模型擬合和訓練樣本的復雜性，尤其是當我們對自己的群體數據不夠了解時，SVM或許是基因組預測的備選方法。
SVM的基本思想是求解能夠正確劃分訓練數據集並且幾何間隔最大的分離超平面。在支持向量回歸（Support Vector Regression，SVR）中，通常使用近似誤差來代替像SVM中那樣的最佳分離超平面和支持向量之間的餘量。假設ε為不敏感區域的線性損失函數，當測量值和預測值小於ε時，誤差等於零。SVR的目標就是同時最小化經驗風險和權重的平方范數。也就是說，通過最小化經驗風險來估計超平面。
下圖1比較了SVM中回歸（圖A）和分類（圖B）的差別。式中ξ和ξ*為鬆弛變數，C為用戶定義的常數，W為權重向量范數，ϕ表示特徵空間映射。
當SVM用於預測分析時，高維度的大型數據集會給計算帶來極大的復雜性，核函數的應用能大大簡化內積，從而解決維數災難。因此，核函數的選擇（需要考慮訓練樣本的分布特點）是SVM預測的關鍵。目前最常用的核函數有：線性核函數、高斯核函數（RBF）和多項式核函數等。其中， RBF具有廣泛的適應性，能夠應用於訓練樣本（具有適當寬度參數）的任何分布。盡管有時會導致過擬合問題，但它仍是使用最廣泛的核函數。
集成學習（Ensemble Learning）也是機器學習中最常見的演算法之一。它通過一系列學習器進行學習，並使用某種規則把各個學習結果進行整合，從而獲得比單個學習器更好的效果。通俗地說，就是一堆弱學習器組合成一個強學習器。在GS領域，隨機森林（Random Forest，RF）和梯度提升機（Gradient Boosting Machine，GBM）是應用較多的兩種集成學習演算法。
RF是一種基於決策樹的集成方法，也就是包含了多個決策樹的分類器。在基因組預測中，RF同SVM一樣，既可用做分類模型，也可用做回歸模型。用於分類時，注意需要事先將群體中個體按表型值的高低進行劃分。RF演算法可分為以下幾個步驟：
最後，RF會結合分類樹或回歸樹的輸出進行預測。在分類中，通過計算投票數（通常使用每個決策樹一票）並分配投票數最高的類別來預測未觀察到的類別。在回歸中，通過對ntree輸出進行求平均。
有兩個影響RF模型結果的重要因素：一是每個節點隨機取樣的協變數數量（mtry，即SNP數目）。構建回歸樹時，mtry默認為p/3（p是構建樹的預測數量），構建分類樹時，mtry為[圖片上傳失敗...(image-10f518-1612450396027)] ；二是決策樹的數量。很多研究表明樹並非越多越好，而且構樹也是非常耗時的。在GS應用於植物育種中，通常將RF的ntree設置在500-1000之間。
當GBM基於決策樹時，就是梯度提升決策樹（Gradient Boosting Decision Tree，GBDT），和RF一樣，也是包含了多個決策樹。但兩者又有很多不同，最大的區別在於RF是基於bagging演算法，也就是說它將多個結果進行投票或簡單計算均值選出最終結果。而GBDT是基於boosting演算法，它通過迭代的每一步構建弱學習器來彌補原模型的不足。GBM通過設置不同的損失函數來處理各類學習任務。
雖然已經有不少研究嘗試了將多種經典機器學習演算法應用於基因組預測中，但提升的准確性仍然有限，而且比較耗時。在無數的機器學習演算法中，沒有一種方法能夠普遍地提高預測性，不同的應用程序及其最優方法和參數是不同的。相比於經典的機器學習演算法，深度學習（Deep Learning，DL）或許是未來應用於基因組預測更好的選擇。
傳統的機器學習演算法如SVM，一般是淺層模型。而深度學習除了輸入和輸出層，還含有多個隱藏層，模型結構的深度說明了它名字的含義。DL的實質是通過構建具有很多隱藏層的機器學習模型和海量的訓練數據，來學習更有用的特徵，從而最終提升分類或預測的准確性。DL演算法的建模過程可簡單分為以下三步：
在GS領域，研究較多的DL演算法，包括多層感知器（Multi-layer Perceptron，MPL）、卷積神經網路（Convolutional neural network，CNN）和循環神經網路（Recurrent Neural Networks，RNN）等。
MLP是一種前饋人工神經網路（Artificial Neural Network，ANN）模型，它將輸入的多個數據集映射到單一的輸出數據集上。MLP包括至少一個隱藏層，如下圖2中所示，除了一個輸入層和一個輸出層以外，還包括了4個隱藏層，每一層都與前一層的節點相連，並賦予不同權重（w），最後通過激活函數轉化，將輸入映射到輸出端。
CNN是一類包含卷積計算且具有深度結構的前饋神經網路，通常具有表徵學習能力，能夠按其階層結構對輸入信息進行平移不變分類。CNN的隱藏層中包含卷積層（Convolutional layer）、池化層（Pooling layer）和全連接層（Fully-connected layer）三類，每一類都有不同的功能，比如卷積層的功能主要是對輸入數據進行特徵提取，池化層對卷積層特徵提取後輸出的特徵圖進行特徵選擇和信息過濾，而全連接層類似於ANN中的隱藏層，一般位於CNN隱藏層的最末端，並且只向全連接層傳遞信號。CNN結構如下圖3所示。
需要注意的是，深度學習不是萬能的。使用DL的前提是必須具有足夠大和質量好的訓練數據集，而且根據GS在動植物方面的研究表明，一些DL演算法和傳統的基因組預測方法相比，並沒有明顯的優勢。不過有一致的證據表明， DL演算法能更有效地捕獲非線性模式。因此，DL能夠根據不同來源的數據通過集成GS傳統模型來進行輔助育種。總之，面對將來海量的育種數據，DL的應用將顯得越來越重要。
以上是GS中常見的預測模型，不同分類方式可能會有所區別。這里再簡單介紹一下上述未提及到但比較重要的方法，其中一些是上述三類方法的拓展。
再生核希爾伯特空間（Reprocing Kernel Hilbert Space，RKHS）是一種典型的半參數方法。它使用高斯核函數來擬合以下模型：
RKHS模型可採用貝葉斯框架的Gibbs抽樣器，或者混合線性模型來求解。
GBLUP仍然是動植物育種中廣泛應用的方法，它假定所有標記都具有相同的效應。但在實際情況中，任何與目標性狀無關的標記用來估計親緣關系矩陣都會稀釋QTL的作用。很多研究對其進行改進，主要有幾種思路：
沿用以上的思路，sBLUP(Settlement of Kinship Under Progressively Exclusive Relationship BLUP, SUPER BLUP)方法將TABLUP進一步細化為少數基因控制的性狀，這樣基因型關系矩陣的構建僅僅使用了與性狀關聯的標記。
如果要在親緣關系矩陣中考慮群體結構帶來的影響，可根據個體遺傳關系的相似性將其分組，然後將壓縮後的組別當做協變數，替換掉原來的個體，而組內個體的親緣關系都是一樣的。因此在構建基因組關系矩陣時，可用組別的遺傳效應值來代替個體的值，用個體對應的組來進行預測，這就是cBLUP（Compressed BLUP）。
以上思路都提到了將已驗證和新發現的位點整合到模型中，這些位點從何而來？最常見來源自然是全基因組關聯分析（Genome Wide Association Study, GWAS）。GS和GWAS有著天然的聯系，將GWAS的顯著關聯位點考慮進GS中，直接的好處是能維持多世代的預測能力，間接的好處是能增加已驗證突變的數量。
下圖比較了GWAS輔助基因組預測的各類方法比較。a表示分子標記輔助選擇方法（MAS），只利用了少數幾個主效位點；b表示經典GS方法，利用了全部標記，且標記效應相同；c對標記按權重分配；d將顯著關聯標記視為固定效應；e將顯著關聯標記視為另一個隨機效應（有其自身的kernel derived）；f將染色體劃分為片段，每個片段構建的G矩陣分配為不同的隨機效應。
GWAS輔助基因組預測的結果會比較復雜，單純地考慮將關聯信號納入模型不一定都能提高准確性，具體表現應該和性狀的遺傳構建有關。
GS對遺傳效應的估計有兩種不同的策略。一是關注估計育種值，將加性效應從父母傳遞給子代。而非加性效應（如顯性和上位性效應）與特定基因型相關，不能直接遺傳。當估計方差組分時，非加性效應通常和隨機的環境效應一起被當成噪音處理。另一種策略同時關注加性和非加性效應，通常用於雜種優勢的探索。雜交優勢一般認為是顯性和上位性效應的結果，因此，如果非加性效應很明顯，而你恰好將它們忽略了，遺傳估計將會產生偏差。
雜種優勢利用是植物育種，尤其是水稻、玉米等主糧作物的重要研究課題。將非加性遺傳效應考慮進GS模型進行雜交種預測，也是當前基因組預測在作物育種中研究的熱點之一。
當然，雜種優勢效應的組成也是隨性狀而變化的，不同性狀的基因組預測需要與鑒定雜優QTL位點結合起來。由於一般配合力GCA（加性效應的反映）和特殊配合力SCA（非加性效應的反映）可能來自不同遺傳效應，所以預測雜交種F 1 應該分別考慮GCA和SCA。GCA模型可以基於GBLUP，重點在基因型親緣關系矩陣構建。SCA模型有兩種方法：一是將雜優SNP位點的Panel作為固定效應整合進GBLUP模型中；二是使用非線性模型，如貝葉斯和機器學習方法。據報道，對於加性模型的中低遺傳力性狀，機器學習和一般統計模型比較一致。但在非加性模型中，機器學習方法表現更優。
傳統的GS模型往往只針對單個環境中的單個表型性狀，忽略了實際情況中多性狀間或多環境間的相互關系。一些研究通過對多個性狀或多個環境同時進行建模，也能提高基因組預測的准確性。以多性狀（Multi-trait，MT）模型為例，多變數模型（Multivariate model，MV）可用如下公式表示：
多性狀選擇一般用於性狀間共有某種程度的遺傳構建，即在遺傳上是相關的。尤其適用於對低遺傳力性狀（伴隨高遺傳力性狀相關）或者難以測量的性狀。
農作物的環境條件不如動物容易控制，而且大部分性狀都是數量性狀，很容易受到環境影響。多環境（Multi-environment，ME）試驗發揮了重要作用，基因型與環境互作（Genotype by E nvironment，G × E）效應也是當前基因組選擇關注的焦點。
除了GBLUP，多變數模型也可基於貝葉斯框架的線性回歸，或者基於非線性的機器學習方法。
我們知道，基因經過轉錄翻譯以及一系列調控後才能最終體現在表型特徵上，它只能在一定程度上反映表型事件發生的潛力。隨著多組學技術的發展，整合多組學數據用於基因組預測也是目前GS研究的一個重要方向。
在植物育種中，除基因組外，轉錄組學和代謝組學是當前GS研究相對較多的兩個組學。轉錄組將基因表達量與性狀進行關聯預測，代謝組則將調控表型的小分子含量與性狀進行關聯預測，對於某些特定的性狀而言，可能會提高預測能力。最好的方法是將各個組學的數據共同整合進模型，但這樣會大大增加模型的復雜度。
表型測定的准確性直接影響模型的構建。對於一些復雜性狀，單憑肉眼觀察記錄顯然已不可取，而且表型調查費時費力，成本很高。因此，高通量表型組也是GS發展的重要方向。表型的范疇非常之廣，當個體性狀不可簡單測量時，我們也可採用多組學數據，如蛋白組、代謝組等數據來替代。
考慮到成本效益問題，多組學技術在動植物育種中仍處於研究階段，但代表了未來的應用方向。
㈡ 什麼是BLUP法
BLUP法即最優線性無偏估計法，是1973年由美國提出的一種評定種公牛育種值的方法。它的基本出發點是從女兒的表型值（產奶量）中將公牛育種值剖分出來；也可將牛群效因或來源效因剖分出來，這樣所得的公牛育種值（公牛效因）消除了牛群差異的影響，其估測精確度高，誤差最小，可用線性函數表示。
BLUP法將所有重要的系統環境影響和遺傳分組的固定效應都考慮時，可以通過混合模型方程組得到最准確而又可靠的個體育種值的預測值。由於方程組的解與BLUP估值等價，所以在家畜實際育種中，混合模型方程組法已成為BLUP法的同義詞。當一個混合模型中隨機遺傳效應為該個體本身的加性遺傳值，即育種值時，該模型為個體動物模型，也稱為加性遺傳值模型。
BLUP法的動物模型公式為：
Y＝Xb＋Za＋e式中，Y——某動物的表型觀察值；X——與固定效應有關的個體數矩陣；b——固定效應的估計值（包括場、年、季、胎次等）；Z——與加性遺傳效應有關的個體數矩陣；a——需估計的某動物的育種值（加性遺傳效應）；e——隨機誤差。
BLUP動物模型有以下優點：
（1）親屬資料的最佳利用，將父母、本身、旁系及後代資料有機地結合。
（2）可校正選配造成的配偶誤差。
（3）當使用多泌乳期記錄時，可將由淘汰造成的偏差校正到最小。
（4）最重要的是可校正牛群、年份、季節、性別、胎次、地區等固定環境效應並進行育種值最佳無偏估計。
動物模型BLUP法計算過程比較復雜，必須用電腦才能完成繁瑣的計算。所以，雖然它具有上述優點，但目前在我國仍未真正用於奶牛育種工作。
㈢ 恐龍的重量是怎樣測出來的
怎樣估算恐龍身體重量？ 由美國古脊椎動物學家、恐龍專家科爾伯特改進的用復原的模型測量體重的方法，是目前世界上公認的行之有效的方法，原理非常簡單，即體重等於體積乘以比重，恐龍的體積是利用恐龍的模型求出來，首先依據恐龍的骨架做一個縮小的模型，再把模型放入一個箱子里,把箱子用沙子添滿，然後量出箱子和沙子的體積，二者之差就是恐龍模型的體積。把恐龍模型的體積的數字乘以原大與復原模型的比例的三次方，就得出了恐龍生前的體積。用恐龍的體積乘以現代與恐龍親緣關系較近的爬行動物的比重，最後就得出了恐龍的大致體重。
日前，英國曼徹斯特大學的生物學家開發了一種新技術，能藉助激光掃描更精確地測量恐龍的體重和大小，並發現它們並未像之前預想的那樣重。科研人員同時表示，這一技術可應用於所有恐龍體重的測量。相關研究論文發表在近日出版的《生物學快報》雜志上。
論文的主要作者比爾·塞勒斯博士說：「對於古生物學家而言最重要事情之一，就是需要了解已經變為化石的動物有多重。這項工作困難得出乎意料，因此我們測試了這種新途徑。」科學家利用激光掃描了馴鹿、北極熊、長頸鹿和大象等14種大型哺乳動物的骨架，並計算出了包裹它們骨架主要部分的最少皮膚量。
他們發現，這些動物的實際重量比上述的計算量大概多出了21%。隨後，他們將這一方法應用在了德國柏林自然博物館龐大的腕龍科恐龍——長頸巨龍骨架上，通過激光掃描計算出這種恐龍的皮膚和骨架包裹量，再增加21%，從而計算出它的實際重量。此前對於長頸巨龍重量的估測十分多樣，最高可達80噸，但研究小組此次的計算可將這一數字降低至23噸。雖然從現有狀況較好的骨架來看，最高和最重的恐龍仍可能是長頸巨龍，但這種龐然大物卻比之前預想的輕了許多……
㈣ 曹沖稱象還有別的辦法可以稱到象的重量嗎
可以找個游泳池，水高過大象，量下原始水位，將大象踢入池中，重新計算水位。
可以算出體積。
再想辦法估出密度，重量即可得到。
拉開冰箱門
把大象塞冰箱里
關上冰箱門
然後稱冰箱的重量
減去冰箱原來的重量即可(冰箱的自重可以咨詢廠家獲得)
把稱拿來
把大象放在稱上，稱之
把稱收起來
將大象從高樓上扔下，測量其落地的時間。
用牛頓第二運動定律，易得其質量m. (很明顯這位兄弟學的是鴨梨世鐸得的物理學…)
題目要求比曹沖先進的方法，你們都沒有
突出「先進」二字。 我的辦法是，用太空梭把大象送到一個
小行星上用普通稱稱，利用小行星軌道用銀河巨型計算機算出小行星質量，再用廣義相對
論算出大象重量。
挖個足夠大的坑，把大象趕進去，之後填土
最後剩下多少土方乘以大象的密度就可以了
讓大象以某速度沖撞某儀器，測量其沖量，除以速度就是質量.
稍微改進一點，不用石頭，還得搬來搬去的，讓幾個人跳上船去，再讓他們自己報體重。
(嗯，當年一休哥就是這么稱大象的，結果差點被宰了)
切一小塊象肉求其密度
把大象放進浴缸中，漫出的水的重量就是其體積
當年阿基米德 就是這樣算王冠重量的
把大象煮煮吃了，看能味飽多少個人，計算方式就是：
5兩×男生人數＋2兩×女生人數＝大象體重. (大象里邊又不都是肉，還有屎啥的,象牙，
象骨，象毛去了才好算象肉的密度)
重新定義質量單位，
以一隻大象的重量定為一噸，
所以所有的大象質量都是一噸了
恩
思考良久，覺得杠桿原理比較合適。
嗯，大家知道了吧，我用的是桿秤。
切了塊稱，一斤一塊!
用一氫氣球充氣後使其恰好讓大象懸浮,餘下就慢慢算了
把大象趕到地秤（就是警察查汽車超載的那種）上，不就可以直接得到其體重嗎？小曹那
時候沒有這么先進的東東，所以只好用「巧妙」的笨方法了：）
來一頭母象，誘惑它說
找個冰面放上木板看看用多大的力拉動，算出摩擦系數，再趕上大象再用一個力拉動可能
算出質量
帶大象去滑雪
把大象無限微分，處處可導
讓大象去pielove徵友 一般要說身高體重...
用西遊記裡面隨便找個法寶收了大象，然後稱前後重量變化. 對了，說起這個我想起來
那個一叫名字就收人的葫蘆是哪個神仙的來著？
不知道他走在路上怎麼和人打招呼。。。
找一個承重N噸的橋
趕n只大象上去
橋塌
大象重量=N/n-1
嗯 橋載重越大 最後得出的數值越准
打電話或者發電子郵件問問曹沖，或者乾脆google一把
嚴肅點，嚴肅點，我們這兒稱象呢！！
說個比較正規的吧，找一輛卡車，托著大象，行駛一段距離，看用了多少油（質量
和用油的模型應該有吧，沒有就建一個）。
在曹沖的船上綁上細的鐵絲，將象趕上去，測量電線的電阻的變化，根據電阻和長度、截
面面積的關系，計算出長度的變化，然後做有限元分析，推算出重量，然後轉換為質量
從大象出生起，每日記錄大象的吸收與排泄，加上皮膚散發，呼吸帶出的蒸氣，這樣就可
得到大象的精確質量。
估計完了之後就是一個動物生理學家。
你找個蹺蹺板，把大象放在一段
然後從一百層樓上跳下，看大象能彈多高
然後就行了
雙規大象，看它招不招
哼哼，把大象用油燒成渣渣。收集排放的CO2和H2O，用油量以及O2量。
渣渣－用油量－O2重量＋H2O質量＋CO2質量＝大象重量
騙大象說它太瘦了，他只有1克，他氣得不行了。
就說：「很豐滿的，體重是****kg」
用錄音機錄下來，然後找人做錄音整理，就可以得出結論了。
㈤ 迷你世界如何製作3000血的生物
製作微縮模型，生物模型方塊修改，替換生物模型。
製作微縮模型後保存，通過背包工具裡面將生物模型方塊找到對其進行編輯，然後點擊+號添加微縮模型，選擇好對應的身體部件進行替換，最後進行調整就好了。
在編輯模型人偶時，微縮模型的大小，方向旋轉，位置都可以設置，微縮模型是無法直接替換生物模型，不能直接就去點替換，要去資源庫裡面操作，修改新建生物必須要載入對應的新生物插件。
㈥ 5 - 混合模型(MME)的多種變型
本文從最簡單的單隨機變數到多個隨機變數的復雜模型。
最簡單的基礎模型(每個動物有一個記錄)：
y = Xb + Za + e
其中E(y) = Xb; E(a) = E(e) = 0
ZGZ'是遺傳方差，Z連接個體G與記錄
Var(y) = V = var(Za + e) = ZGZ' + R; cov(y, e) = ZG; cov(a, y ) = GZ'；
再次寫下，經過Henderson轉為MME的一般式子
Mixed model equations(MME)
這個模型常常被稱為動物模型（animal model）:一般為單性狀，每個動物最多一個記錄，對每個動物直接估計育種值。但是需要注意的是，動物模型不是第一種混合模型。
根據上述的式子，求解時，需要求G或C的逆矩陣， 因計算量過大， 1960年以前，不可能能實現。
所以人們就先簡化了MME
所以的公牛都沒有親緣關系
但是其I表示的不合理， 應有關系矩陣
加上公牛之間的關系矩陣A s
問題：計算3個牛場的牛體高的平均值
結果：
但我們知道這些奶牛的是來自兩頭公牛的後代：
這樣OLS的前提假設（記錄之間不相關）就不成立了。
因為：具來自同一個sir的奶牛會傾向有相似的表型，其均有公牛的一部分相關的基因。
所以需要在各記錄之間加入方差-協方差
(co)variances是描述記錄值之間的相關
variance => 所有的變異來源
covariance => 動物之間的共同差異（如有相同的父親）
如有共同父親，則相同父親傳遞的遺傳力基本相同
但是這里也就需要知道這些新的參數：
假設：
但是公牛的貢獻在哪？將通過以下實例（顏色）進行展示：
生成公牛的（協）方差遺傳矩陣G *
增強（總，來自遺傳方差的遺漏部分）環境（協）方差Var(e * ) = R*
所以總（表型）（協）方差矩陣變為V（SI中稱為P）
當OLS考慮（協）方差矩陣時，其變為了GLS
我們引入Z(10 x 2)矩陣，使其關聯公牛（無親緣關系）通過其女兒的記錄
G* = ZVar(s)Z' = Z(I x 25)Z『 = ZZ『 x 25(前面我們已經假設了遺傳方差為100)
前面我們對R和G進行假設：
再在SI中測試：
則結果：
當具有相同父親或者母親的協方差為25.
當具有相同父親，與母親的協方差為50.
Note: 以下R與G都是指的父親模型
sir model 變為animal model演算法：
(1) 方差-協方差矩陣不同；
(2) 非對角線元素在非近交群體和非克隆體時，最大為50
(3) 對角線元素在非近交群體,當所有遺傳方差被解釋時為100
是對sir model的升級，增加了母親的遺傳效應。但是直接加入母親效應，模型當時沒法計算，就加入了外祖父效應，間接反映母親效應：
（1）直接加入目前效應m:
模型為：
但隨著近30年的演算法和計算機的發展，可以直接計算A及其逆矩陣
所以我們現在普遍使用動物模型，但有時也需要根據數據調整。
最基本的動物model：
1 增加「環境」效應（random）， 如： 觀察值之間具有協方差（豬的窩， 相同母親（非遺傳））， 重復觀察值（奶牛的泌乳量）， 這些都需要使用(co)variances
其中加入永久環境隨機效應： y = Xt + Zg + Zp + e， 這樣會使e = e initial -Zp
2 有時，會受到共同遺傳因素的影響，如：直接母性作用（calving ease, 斷奶前體重）, 繁殖性能（fertility）
3 random regressions 與（協）方差結構相關的回歸效應, eg:年齡， DIM， 通常用於分析軌跡的「縱向」數據的模型，使用的模擬軌跡的函數最常用的有LP, spline.
4 具有遺傳分組的動物模型
系譜中有些個體買的親本未知，通常把這些個體成分基礎群體（base population animals）。動物模型中，親緣關系矩陣是假設所有的個體從一個總體中抽樣獲得。但是基礎群體的遺傳方差都是不同的，需要對這些未知父母分組。目前也是一個應用研究熱點。