Waddington D, Burt D W
AFRC Institute of Animal Physiology and Genetics Research, Edinburgh Research Centre, Roslin, Midlothian, UK.
Comput Appl Biosci. 1992 Dec;8(6):539-48. doi: 10.1093/bioinformatics/8.6.539.
A generalized linear model with Gamma errors is used to estimate the coordinates of a restriction map when the site order is known. This can be conveniently programmed in a wide range of statistical packages (e.g. Genstat 5, Minitab, SAS), and gives maximum likelihood estimates with their associated optimal properties. Regression diagnostics allow the checking of assumptions and help to identify mis-specified, influential or discordant fragment lengths. A specific diagnostic for identifying fragment lengths causing reversal of restriction site order is derived. Exact 'fragment' lengths from DNA sequencing can be conveniently included in an approximate manner by giving them a larger weight than observed restriction fragment lengths. Two examples and the Genstat 5 codes used in their analysis are presented.
当位点顺序已知时，使用具有伽马误差的广义线性模型来估计限制酶切图谱的坐标。这可以方便地在多种统计软件包（如Genstat 5、Minitab、SAS）中进行编程，并给出具有相关最优性质的最大似然估计。回归诊断允许检查假设，并有助于识别指定错误、有影响或不一致的片段长度。推导了一种用于识别导致限制酶切位点顺序颠倒的片段长度的特定诊断方法。通过给予来自DNA测序的精确“片段”长度比观察到的限制酶切片段长度更大的权重，可以方便地以近似方式将其纳入分析。给出了两个例子及其分析中使用的Genstat 5代码。
