有限差分近似法在数量性状基因座定位上最大概似估值变异矩阵估算.PDFVIP

研究報告 有限差分近似法在數量性狀基因座定位之應用 201 有限差分近似法在數量性狀基因座定位上 最大概似估值變異矩陣估算上之應用 江欣容、劉 清* 國立臺灣大學農藝系 摘要 資料分布為何，利用概似函數的二次微分式 與有限差分近似法的近似二次微分式，兩者 在數量性狀基因座的定位與分析上，混 所算出的漸近變異矩陣幾乎完全相同。因此 合 模 式 的 最 大 概 似 估 值 可 以 EM 法 在數量性狀基因座的定位與分析上，當各類 (expectation maximization, Dempster et al. 數學模式之概似函數無已知的二次微分式或 1977) 、ECM 法 (expectation-conditional 不易以解析方法導出時，建議可先以 EM 法、 maximization, Meng and Rubin 1993) 、IRLS ECM 法、IRLS 法等方法來計算模式之最大 法(iteratively reweighted least squares)等 概似估值。一旦求得最大概似估值的解析 方法求出，但最大概似估值的漸近變異矩陣 解，則可利用數值方法上的有限差分近似法 因為需要利用概似函數的二次微分式求算， 來計算概似函數的近似二次微分式及最大概 又由於一般混合模式之概似函數的二次微分 似估值的漸近變異矩陣。 式相當複雜不易導出，因此如何計算混合模 關鍵詞︰數量性狀基因座、遺傳標識、簡單 式之最大概似估值的漸近變異矩陣為一相當 區間定位法、綜合區間定位法、有 重要的課題。因為僅有估值卻不知估值的變 限差分近似法、最大概似估值、漸 異，將無法評估根據此估值所做統計推論的 近變異矩陣。 可靠性。本研究提出以數值方法上的有限差 分近似法 (finite difference approximation Application of Finite Difference method)來計算概似函數的近似二次微分式 Approximation on the Estimation of 及最大概似估值的漸近變異矩陣。為證實利 Dispersion Matrix of Maximum 用有限差分近似法所算之漸近變異矩陣結果 Likelihood Estimates for Mapping 的正確性，特將模擬之常態分布、二項分布 and Analysis of Quantitative Trait 與卜瓦松分布的 F2 子代的數量性狀資料，分 Loci (QTL) 別 以 概 似 函 數 的 二 次 微 分 式 (