24小时热门版块排行榜    

查看: 1432  |  回复: 1

jiejing

银虫 (初入文坛)

[交流] [求助]Phase2.1 已有1人参与

求助原因: 学习需要
软件名称: Phase2.1
版本号: 2.1
软件公司或主页:
http://www.stat.washington.edu/stephens/home.html
软件功能说明:
Phase 是一个以人口数据为基础构建单倍型频率的软件
官方网站网页:
http://www.stat.washington.edu/stephens/home.html
求助者email: chenjiezs@163.com
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cheungsober

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
http://www.stat.washington.edu/stephens/phase/download.html
下载即可,
PHAS采用Bayesian统计方法,利用群体的基因型数据来构建单倍型。软件可以处理SNP、微卫星及其他多等位形式的位点 (如三态的SNP和HLA等位基因)。
较之先前的版本,它提供了:重组热点模型,并且重组率计算的偏差小
2.        运行条件
虚拟dos环境下运行,基本命令:./PHASE test.inp test.out
3.        默认的输入文件格式:
个体数目(正整数)
位点数目(正整数)
P 位点(1) 位点(2) ``````位点(n) )(位点的位置,默认单位为bp,正整数)
位点基因型(1) 位点基因型(2) ... 位点基因型(n)(以M或S,分别代表二态或三态的等位形式)
ID(1)
基因型(1)(二态位点可以碱基或数字表示,数据缺失以?代替;多态位点必须以正整数表示,缺失以-1代替)
ID(2)
基因型(2)
.
.
.
ID(n)
基因型(n)

4.        输出文件
1)        Freqs:给出单倍型频率及标准偏差
2)        Pairs:给出每个个体最可能的单倍型及P值
3)        Recom:给出一定模式下整个区域重组参数的估计
4)        Hotspot:给出使用热点选项条件下的重组参数的估计
5)        Signif:给出病例-对照空值假设的置换检验的P值
6)        Monitor:估计某一模型下单倍型统计的适用性
7)        Hbg:空文件,忽略

5.        模型的选择
1)        第一种是最初的模型(original model),即非重组模型。由Stephens et al. (2001), and Stephens and Donnelly (2003)提出, 它忽略重组的存在。
第二种是新的模型,明确地允许重组,即重组模型(recombination model)。
它们都需要在-M选项下运行:-MS或-MR。
此外还有混合模型(hybrid model),先用较快的初始模型进行初步的运算,再在用重组模型进行最终运算,运行的选项为-MQ。这种方法在构建单倍型上基本与-MR一样准确,运行起来却更快。其他的条件的比较表明:-MR选项的性能通常更好,因此被作为默认值。
2)        如果不使用-M选项,程序将新的重组模型作为运行的默认选项。这样使计算更为精确些,但是运算量也加大了。
3)        如果数据信息相对较少(小于20个个体,10个位点),采用-MR法能得到更可靠的LD的估计。
2楼2011-03-25 15:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiejing 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[硕博家园] T^T +6 m584r2z4j3 2026-06-11 7/350 2026-06-13 23:10 by 3eaav2pea1
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[公派出国] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:33 by lkmi35jdzt
[考研] T^T +6 m584r2z4j3 2026-06-11 6/300 2026-06-13 18:29 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 14:45 by 80j16ujcu8
[考研] O(∩_∩)O +3 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 14:42 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[考博] →_→ +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 08:49 by t4hz8376wy
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +6 yexuqing 2026-06-07 7/350 2026-06-12 16:46 by wwwkkk83
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见