零代码不求人,GWAS在线分析网站,建议收藏

公众号:生信漫谈,获取最新科研信息!

零代码不求人,GWAS在线分析网站,建议收藏icon-default.png?t=N7T8https://mp.weixin.qq.com/s?__biz=MzkwNjQyNTUwMw==&mid=2247485814&idx=1&sn=d79745ce0254b32662b6566f2a6fc798&chksm=c0e9ec86f79e6590bfceef1835dabf657fcbee1eafc81adac8dcc6fb9ba859270871856be228&token=1502972688&lang=zh_CN#rd

最新杂交水稻资源

最近得益于韩斌院士发表NG提供的最新杂交水稻资源构建的在线分析网站:http://ricehybridresource.cemps.ac.cn/#/,功能还是非常强大的(主页me哪里显示怪怪的,但是不影响分析)。

添加图片注释,不超过 140 字(可选)

前期相关推文:给你一个国家杰青资助金额,你会进行韩斌院士水稻杂种优势的研究工作吗

1变异位点查询

其实这个网站的功能还是非常强大,可以查看杂交品种的变异位点,如下所示:

添加图片注释,不超过 140 字(可选)

结果如下所示:

添加图片注释,不超过 140 字(可选)

添加图片注释,不超过 140 字(可选)

2样品信息查询

具体的样品代号信息可以通过这个查询:

添加图片注释,不超过 140 字(可选)

添加图片注释,不超过 140 字(可选)

重点是这个GWAS部分

3 GWAS部分

通过输入物种信息,具体表型信息就可以得到结果:

添加图片注释,不超过 140 字(可选)

我们以分蘖角度做关联分析:

添加图片注释,不超过 140 字(可选)

可以看到具体的表型数据,以下是具体的结果,可以看到在9号染色体这个区段有一个非常显著的QTL位点,效应非常强:

添加图片注释,不超过 140 字(可选)

表型数据还是非常新的,稻米品质类数据都有:

添加图片注释,不超过 140 字(可选)

产量数据:

添加图片注释,不超过 140 字(可选)

最后用生育期数据做一个GWAS分析图:

添加图片注释,不超过 140 字(可选)

具体的GWAS结果也可以下载:

添加图片注释,不超过 140 字(可选)

ricevarmap

相较于韩斌院士的杂交水稻GWAS分析网站,华中农业大学的这个网站,应该大部分比较熟悉,但还是有很多人不会使用。

网站地址:http://ricevarmap.ncpgr.cn/

添加图片注释,不超过 140 字(可选)

网站的功能也是非常强大的,我们今天专门看下GWAS部分:

添加图片注释,不超过 140 字(可选)

选择所有种质资源的生育期表型进行全基因组关联分析:

添加图片注释,不超过 140 字(可选)

结果如下,提供了表型分布图:

添加图片注释,不超过 140 字(可选)

表型的具体信息:

添加图片注释,不超过 140 字(可选)

添加图片注释,不超过 140 字(可选)

CottonMD

华中农业大学构建的棉花数据库,数据库内容非常丰富,功能齐全。

网站地址:https://yanglab.hzau.edu.cn/CottonMD

添加图片注释,不超过 140 字(可选)

功能模块:

添加图片注释,不超过 140 字(可选)

GWAS、eQTL等在线分析:

添加图片注释,不超过 140 字(可选)

添加图片注释,不超过 140 字(可选)

结果分析图:

添加图片注释,不超过 140 字(可选)

 

棉花数据库还是非常齐全,分析图片也非常精美,希望对大家有用!

回顾往期内容

1、生物博士这些研究网站建议收藏

2、生物博士互作蛋白研究收藏网站汇总

3、Cell文章教你怎么找互作蛋白

4、超简单的比较基因组学Ka/Ks进化速率图绘制

5、超详细非模式物种GO数据库全新代码构建

6、生物研究的一些常用在线生信工具建议收藏

7、植物基因的时空表达及胁迫处理可视化网站

8、生信小白送你超好用分子克隆工具SnapGene

9、三个团队三个故事,一个基因连续三年发表三篇高水平文章

10、从西湖大学柴大教授看当代博士研究生的自我突围

11、如果浙江大学李兰娟院士在MDPI杂志发表的SCI论文结果不能复现,会怎么样?

12、为什么浙江大学李兰娟院士选择在MDPI杂志发表SCI论文

13、施一公为什么没有推荐颜宁/柴继杰成院士候选人,李兰娟推荐了谁?

14、顶尖学者纷纷加入西湖大学,颜宁什么时候加入施一公团队

15、由万建民院士Cell论文引起的一点思考和疑惑

16、清华大学施一公团队为什么选择《International Journal of Molecular Sciences》发表SCI

17、西湖大学施一公团队在国际著名杂志IJMS发表重要论文给予我们的启示

18、可以在施一公团队投稿的《International Journal of Molecular Sciences》上发表文章吗?

通过以上步骤是否学会了呢?关注我们,学习更多生信干货,原创不易,请各位多多点赞加转发朋友圈进行分享!

今天就先给大家介绍到这里,希望大家的科研能有所帮助!祝您科研顺利快乐!

生信漫谈

生信漫谈,认识生信,学习生信,跨越生信入门路上的障碍,从而利用生信技术解决科研学习路上的绊脚石!
 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
GWAS meta分析是一种将多个基因组关联研究结果进行综合分析方法,以下是一个Python脚本的例子,用于进行GWAS meta分析: ```python import pandas as pd import numpy as np import statsmodels.api as sm # 读取所有研究的GWAS结果文件 study1 = pd.read_csv("study1.csv") study2 = pd.read_csv("study2.csv") study3 = pd.read_csv("study3.csv") # 将每个研究的p值进行变换,转化为z值 study1["z"] = np.sqrt(2) * sm.stats.proportion.proportions_ztest(study1["n_cases"], study1["n_total"], value=study1["OR"])[0] study2["z"] = np.sqrt(2) * sm.stats.proportion.proportions_ztest(study2["n_cases"], study2["n_total"], value=study2["OR"])[0] study3["z"] = np.sqrt(2) * sm.stats.proportion.proportions_ztest(study3["n_cases"], study3["n_total"], value=study3["OR"])[0] # 将每个研究的z值和样本量进行合并 meta_data = pd.concat([study1[["SNP", "z", "n_cases", "n_total"]], study2[["SNP", "z", "n_cases", "n_total"]], study3[["SNP", "z", "n_cases", "n_total"]]]) # 计算每个SNP的z值和加权样本量 meta_data["wz"] = meta_data["z"] * np.sqrt(meta_data["n_cases"] + meta_data["n_controls"]) meta_data["w"] = np.sqrt(meta_data["n_cases"] + meta_data["n_controls"]) # 计算meta分析的z值和p值 meta_z = meta_data["wz"].sum() / meta_data["w"].sum() meta_p = 2 * (1 - sm.stats.norm.cdf(abs(meta_z))) # 输出meta分析结果 print("Meta-analysis result:") print("z value: ", meta_z) print("p value: ", meta_p) ``` 这个脚本首先读取每个研究的GWAS结果文件,将每个研究的p值转化为z值,然后将每个研究的z值和样本量进行合并,计算每个SNP的加权z值和加权样本量,最后计算meta分析的z值和p值,输出结果。需要注意的是,不同研究的样本量、OR值等参数可能存在差异,需要根据具体情况进行调整。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值