Presentation is loading. Please wait.

Presentation is loading. Please wait.

利用SAS统计软件定位 家畜群体的QTL 报告人: 殷 宗 俊 安徽农业大学动物科技学院.

Similar presentations


Presentation on theme: "利用SAS统计软件定位 家畜群体的QTL 报告人: 殷 宗 俊 安徽农业大学动物科技学院."— Presentation transcript:

1 利用SAS统计软件定位 家畜群体的QTL 报告人: 殷 宗 俊 安徽农业大学动物科技学院

2 报告内容 用于QTL检测定位的资源群体类型. 家畜群体QTL检测定位的方法和手段. 利用SAS统计分析软件进行QTL检测定位. 实例分析.

3 一、用于QTL定位的资源群体 资源群体构建(试验设计)→数据收集→构建标记连锁图谱→标记-QTL连锁分析或连锁不平衡分析及定位。

4 资源群体类型 近交系杂交 分离群体 BC 设计 TC设计 F2 设计 重组近交 FS 设计 HS 设计 孙女设计 动物模型

5 二、家畜群体QTL定位的方法和手段 现状 目前,用于分离群体QTL区间定位的统计方法主要有最大似然法和贝叶斯方法等,不同的方法也都形成了很多相应的分析软件,但从现有软件的使用情况来看,在使用上存在许多不便: 专用性太强; 对分子与表型数据的录入要求十分严格; 输出的结果不易被完全理解,使用起来有一定的难度。 )

6 SAS统计软件的强大功能 SAS统计分析软件具有强大的计算分析功能,并具有易编程和直观易懂等优点。 Proc mixed;
Proc nlin; Proc ilm; SAS Macro

7 三、回交群体的QTL分析 1、SAS程序(最小二乘分析) Proc Nlin;
R=0.30; /* Recombination fraction from the linkage map */ parms a=1 r1=.15 mu=0; bounds 0 <= r1 <= R; r2=(R-r1)/(1-2*r1); if M=1 and N=1 then x=r1*r2/(1-R); if M=2 and N=2 then x=1-r1*r2/(1-R); if M=1 and N=2 then x=r1*(1-r2)/R; if M=2 and N=1 then x=r2*(1-r1)/R; model y=mu-a*x; run;

8 运行结果 Gauss-Newton Iterative phase Iter r1 a Squares 0 0 1.0000 2069.6
… … … …

9 标准差与置信区间 Approx Parameter Estimate Std Error CI(95%)

10 四、半同胞群体的QTL分析 分析模型: 通用SAS程序编制

11 混合模型过程的最大似然求解 proc mixed data=sire method=ml; class sire;
model y=sire; run; proc mixed data=sirel method=ml; by l1; random z1*sire;

12 运行结果 Obs L0 interval L1 LA LRT 1 24062.46 1 0.0001 24028.12 34.3348

13 LRT剖面图

14 五、应用总结 SAS软件能进行一般群体的QTL定位研究。 Proc iml 和SAS宏能解决极复杂计算过程。
可以在统计基因组学的多个方面扩展应用。

15 谢谢大家!


Download ppt "利用SAS统计软件定位 家畜群体的QTL 报告人: 殷 宗 俊 安徽农业大学动物科技学院."

Similar presentations


Ads by Google